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SYNOP SI S 


"Studies in Stochastic Modelling and Simulation 
of Mixed Vehicular Traffic"- a thesis submitted in partial 
fulfilment of the requirements for the Degree of Doctor of 
Philosophy by Baldev Raj Marwah to the Department of Civil 
Engineering/ Indian Institute of Technology/ Kanpur, 

July, 1976. 

Characterisation and estimation of traffic are two 
important aspects in the planning and design of a highway 
transportation system. The capacity of a highway is generally 
expressed in terms of passenger car volume per hour. In 
developing countries like India, the traffic is of mixed nature 
consisting of both slow and fast moving vehicles between 
which there are wide variations in speed. Ministry of 
Transport (MOT) and Central Road Research Institute(CRRI) in 
India specify the Passenger Car Equivalents (pCEs) for 
different categories of vehicles. These PCES are considered 
constants and they do not account for the variations in 
traffic composition and volume. 

This study considers the nonlinear interaction 
between the different vehicles in mixed traffic flow situa- 
tions and attempts to characterise the variations of pCEs 
in terms of traffic corrposition and Volume, Traffic flow is 
stochastic in nature, the free speeds are probabilistic and 
the process is affected by logical decisions concerning 
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acceleration, retardation, overtaking etc. Analytic 
techniques of traffic flow theory like hydrodynamic analogy, 
car following theory, queuing theory etc,, are suitable 
for the analysis of simple homogeneous vehicular traffic and. 
it does not seem feasible to use exact analytic procedures 
for the analysis of complex mixed traffic flow. Computer 
simulation is adopted in this study to analyse the mixed 
traffic flow and infer the PCES under different volumes 
and compositions of traffic. 

Traffic volume on a highway has hourly, daily and 
monthly variations and may have secular trends and random 
fluctuations. An understanding of the stochastic charac- 
teristics of traffic volume is needed to determine the 
design volume. Using historical monthly and daily traffic 
data of seven years on five main highways approaching Kanpur, 
stochastic models have been developed for the monthly and 
daily traffic data for three categories of goods carriers 
viz,, incoming -motor vehicles, transit motor vehicles and 
incoming slow moving vehicles. Two types of stochastic 
models were developed for the time series of monthly traffic 
data. They are: 

(i) A general time series model of the form 

Z-£-= f^(t) + f£(t:) + £3 (^t—i/ Zrfc„2 '*»*»f Zt-k ) + e t 
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whore Z^= traffic flow at txme t? f-^(t) = trend corponent? 

f 2 (t) = cyclic component, persistence component; and 

©t = random corponent. 

The secular trend in the monthly traffic series was 
indentified by moving average method and estimated by 
regression analysis. Trend was eliminated from the data 
series and the trend free series were subjected to corre- 
logram and spectral analyses. It was found that 12 monthly 
cycle was significant in the truck (both incoming and transit) 
traffic and in addition 6 monthly cycle was also significant 
in the bullock cart traffic. The trend and cycle free data 
can be represented by a simple autoregressive model of first 
order. The parameters were estimated by standard procedures 
and this representation was found to be satisfactory for all 
the three categories of traffic on the five roads. 

(ii) a general autoregressive moving average (ARMa) model 
can be r ep resented by the equation 0p(B) ( Z^. -Z ) = 0 q (B) e^ , 

where Z^= value of the variable at time t with a mean value 
Z ; e t = random corponent; B= backward shift operator; and 
0 p (B) , e q (B) are polynominals in B of degrees p and q, 
representing respectively an au toregressive component of 
order p add moving average corponent of order q. 

In case the time scries is non stationary, it is 
possible that the dth difference of (Z^-Z ) may satisfy 
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the ARMA process. This is referred, to as autoregressive 
integrated moving average (ARIMA) process. When seasonal 
effects are inportant# it may be possible to use a general 
multiplicative seasonal ARIMA model of order (p, d, q) x 
(P, D / Q) s given by 

- D d _ s 

$ p (Bb) 0p 7 s V (2 t - Z ) = 0Q (B )e q (B) a t 

where in addition D= degree of seasonal differencing; 

s = period of season; V = difference operator; 4^ (B s ) # 

(B s ) = polynomials in B s of degrees P and q respectively; 

and at is the white noise process. This model has the 

advantage that it requires fewer parameters than other models 

of comparable complexity. 

For the monthly traffic series# multiplicative 
seasonal ARIMA model of order ( 1,0, 0) x ( 0#1# 1) 12 was 
found to fit the series well. The value of the various 
parameters of the model for three categories of vehicles on 
the five roads were determined by standard procedures. Tests 
for randomness were also conducted. 

Daily fluctuations of traffic volume were also 
analysed for three rep resen tative months on two roads. 

Weekly differencing was found suitable for the daily series. 
The differenced series were found to follow a simple autore- 
gressive process of first order. The stochastic models of 
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monthly ind daily traffic can be used in estimation and 
forecasting of traffic demands. 

Detailed data were needed for simulation and so 
field studies were conducted for ten days round the clock 
on a 2.25 km stretch of Grand Trunk Road. From the obser- 
vations, peak hourly volumes; variations of volumes within 
the day of six vehicular types; and interarrival time gap 
distributions were estimated. The free speeds/ delay time 
and operating speeds of different categories of vehicles 
were also estimated at varying volume levels and traffic 
conposi tions. 

A conputer simulation model of mixed traffic flow 
on a two lane highway was formulated. Vehicles move from 
either direction of the roadway section abd scanning is done 
at one second intervals. Any vehicle in the roadway moved 
at its free speed where it had enough headway to travel 
unimpeded. In case it cannot move unimpeded and it has a 
higher speed than the vehicle ahead of it, then it may try 
to overtake. Overtaking is possible only if there is no 
conflict with the opposing traffic stream during overtaking 
operation. Otherwise the vehicle is forced to reduce its 
speed to that of the vehicle ahead until sometime later 
overtaking is possible. The above logic is applied to each 
of the vehicles in the section and those entering the section 
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for each time interval. As a vehicle leaves the section/ 
time of leaving is noted ahd its characteristics like travel 
time, delay time and overall running speed etc., are determined. 

Initially historical data of arrival times/ traffic 
corrposition and system parameters were used to identify the 
system and validate the system model. The roadway is 
divided into one metre long sections and vehicles are moved 
from either direction at one second intervals. Vehicles 
move at their free speed when there is enough headway to 
travel unimpeded for two seconds. The change of speed of 
vehicles is assumed to be instantaneous and a looking for 
gap time of two seconds is provided for overtaking operations. 
The overtaking vehicle accelerates only if the speed 
difference wlch the overtaken vehicle is less than 16 krrph. 
Otherwise overtaking takes place at normal speed. For 
vehicles being overtaken, minimum spacing related to speed 
and length of vehicle, is specified. The results from the 
sirrulated model were consistent with the observed character- 
istics and the model is thus validated. Simulation is 
hence found to be a versatile tool for understanding of 
the flow process and identifying the system model. 

Traffic on a highway is highly nonstationary with 
hourly trends/ and cycles, and persistence. However, for 
design considerations a stationary peak rate corresponding 
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to a specified level of risk may be assumed, so this 
study deals with simulation of stationary traffic flow. 

It may be noted that when appropriate/ more complicated 
nonstationary p rocesses can be simulated to derive compara- 
ble results, Conerated data were used in simulation of 
stationary mixed traffic flow at different volumes and 
compositions. When generated data are used in computer 
simulation, it is necessary to initialise the system. In 
this study initially the stretch of the road was assumed to 
be empty. The process was run until a steady state was 
reached. This initial time period varied generally with 
the volume level, being more for lower traffic volumes. 
Simulation was carried out further for a sufficiently long 
time for estimation of the characteristics of the process. 
Simulation runs were performed for tine following 

casts* 

(i) Homogeneous traffic at different volumes for passenger 
cars and trucks taken separately, 

(ii) Mixed traffic at varying volumes and traffic compo- 
sitions for the following combinations; 

(a) five cases consisting of two types of vehicles only,, 
including passenger car and one of the other vehicles 
viz. /truck, motorbike, horse driven vehicle, bulloch- 
cart , or bicycle; 
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(b) four cases consisting of combination of three 
vehicles including car, truck and one of the 
remaining four categories of vehicles? and 

(c) all the six types of vehicles at the peak hourly 
compositions only. 

On an average the processing time for one simulation 
run on IBM 7044/ i4oi system was more than half an hour, in 
all more than loo hours of computer time was consumed. 

When traffic is homogenous, the vehicles move 
essentially at their free speed at low volume levels. As 
the volume increases/ the headways are reduced and there is 
a significant reduction in the operating speed of the 
vehicles till at a certain volume level, the density of the 
roadway section continues to increase thereby jamming the 
traffic. There is a nonlinear relationship between the 
volume level and the average operating speed of the vehicles. 
The capacity of a highway can thus be described for homo- 
geneous traffic at any level of setvice(i. e. , cp era ting speed) 
by knowing the free speed distribution parameters of the 
vehicles. 

In the case of mixed traffic flow, the effect on 
the operating speed of passenger cars was studied at varying 
traffic volumes and composition. The operating speed of 
the car in the mixed flow reduces even at low volume levels 
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due to restrained overtaking operations, The operating 
speed was also found to vary with the composition of slow 
moving vehicles in the mixed flow. A generalised conditional 
rrul tiplicative model was developed to estimate the operating 
speed of cars in mixed traffic and the equivalent volumes 
of the homogeneous traffic at the same level of service 
were determined. Other vehicles affect the operating 
speed of cars and their PCES were estimated for combination 
of two or three types of vehicles at different volumes 
and compositions. The PCEs so determined were not constants, 
but wore fcund to vary nonlin early with the traffic volume 
and composition, curves indicating the relationships 
based on the generalised conditional rail tiplicative model 
are plotted and they may be used for design. 

This study demonstrates that the daily and monthly 
traffic volume can be represented by stochastic models 
which can be used for estimation and forecasting of traffic 
demand. It also demonstrates that the corrpl ex p robl em of 
mixed traffic flow on a two lane highway can be mathematica- 
lly modelled and analysed using corrputer simulation techniques. 
The results of this study may be used for design in the 
neighbour- hood of Kanpur and other comparable regions. For 
regions with significantly different traffic characteristics, 

the methodology demonstrated in this study may be used with 

1 

appropriate modifications. The demand and capacity models can 
also be used in combination for the design of highway systems. 



1. INTRODUCTION 


1.1 G enc ral 

Knowledge of capacity and performance characteris- 
tics of highways is fundamental to traffic engineering. 
Highway capacity concerns itself both with the ultimate 
carrying ability of various facilities and also with the 
relative service characteristics of the facilities opera- 
ting at some fraction of the capacity volume. The study 
of highway capacity is both qualitative and quantitative,, 
which permits evaluation of both the adequacy and the quality 
of the vehicle service being provided by the facility 
(Pignataro, 1973), The deficiencies of an existing highway 
system can be evaluated by coirparing demand volumes to 
the capacity of the existing facilities. The ef f ectivoness 
of changes in the geometries of the highway system can 
also be evaluated with respect to capacity. Design of a 
new facility is based on the capacity analyses coupled 
with the projected traffic demand. Thus traffic demand and 
capacity analyses are two major aspects of highway trans- 
portation system planning. 



1.1.1 Highway Capacity 


Highway capacity is defined as the maximum number 
of vehicles that can pass over a given section of lane 
or roadway during a given time period (mostly an hour) , 
under prevailing roadway and traffic conditions (Pignataro, 
1973) . It is a common practice to express capacity in 
terms of passenger car units (PCUS). Level of service 
is associated with different operating conditions that 
occur on a facility when it accomodates various traffic 
volumes. It is a quantitative measure of the effect of 
a number of factors that include speed and travel time, 
traffic interruptions, freedom to manoeuvre , drivers com- 
fort and convenience etc. six levels of service, designated 
'A to F 1 , have been defined in the Highway Capacity 
Manual (HRB, 19 65) . The capacity of various highways 
have also been defined for each of the six levels of 
service. 

In India, highway capacity has often been 
expressed as the maximum daily traffic volume that can 
operate on the highway. The daily capacities of single 
and double lane highways were considered to be 1000 and 
4000 pcu s respectively (sehgal, 1967) . Very limited field 
studies on highway capacity seem to have been done in India. 



3 


Eastern Region Transport survey Group (1966) has recommended 
a daily capacity of 3000 pCUS for single lane and 10000 PCUS 
for two lane roads. A study team of National Council of 
Applied Economic Research (NCAER) has suggested that single 
and double lane roads could carry 3000 and 7500 PCUS 
respectively per day ( NCAER, 1965) . _The factors for 
reducing capacity to compensate for pavement width and 
lateral obstructions below ideal conditions have also been 
specified by Eastern Region Transport Survey Group (1966) , 

1.1.2 Traffic Demand 

Estimation of future traffic demand is necessary 
for designing the transportation facilities. Future traffic 
demand is generally estimated by a three step procedure: 

(i) Trip generation forecasting; (ii) Distribution of these 
trips between zones, and (iii) Assignment of these trips to 
a future network (Martin et al,, 19 65) , Trips generated 
in a region are dependent upon a number of variables like 
population, income of people, number of vehicles, agricultu- 
ral and industrial production etc. (Roberts, 19 66) . 

Based on present and past traffic data, the structure of trip 
generation model is established and used to obtain an 
estimate of future traffic for some target periods. Traffic 
demand levels are not deterministic, but are probabilistic. 
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Furthermore, the annual forecasted traffic for a design 
year is not constant throughout the year, but there exist 
seasonal cyclic variations within the year. The major 
variations may bo expressed by monthly, weekly and daily 
time patterns of the traffic flow { HRB, 1965) . peaking 
characteristics are observed even within the peak hour. 
The entire phenomenon of traffic flow volume is thus time 
dependent and probabilistic. 

Traffic facilities are usually designed on the 
basis of Design Hourly Volume (DHV) . EXie to the extreme 
variations in traffic flow during the day and throughout 
the year, it was noted in the capacity manual (USDC, 19 5o) 
that s 


.... If a roadway facility is to be so designed 
that traffic will be properly served, consi- 
deration must be given to the brief but frequen- 
tly repeated rush-hour periods. It is neither 
wise nor economical, however, to provide for the 
extreme hourly volumes of traffic that may 
occur but a few times during a year. The law 
of diminishing returns must be applied to fix 
the highest hourly volumes which will justify 
the necessary expenditure of funds to provide 
the added capacity. 

The hourly volumus of traffic may be expressed 
as a percentage of Annual Average Daily Traffic (AADT) . 

The cumulative frequency distribution of hourly traffic 
indicates that at about the thirtieth highest hour, the 
slope of curve changes rapidly and it may not be economical 
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to design for a volume higher than this volume, where 
data are available, or for major locations, DHV may be 
derived from an actual analysis of data, where detailed 
data are not available, data for similar highways may be 
used to determine the DHV by expressing it as a percentage 
of A&DT (wohl and Martin, 1967) . The DHV for the highway 
under consideration is obtained by multiplying this 
percentage by AADT, 

1,2 Statement of the problem 

Traffic flow on a highway is a stochastic process 
which may consist of deterministic components including 
trends and cycles, and probabilistic components. Stochastic 
modelling deals with representing these components in 
analytic terms. Time series analysis considers the sequence 
of values of the variable in time (Box and Jenkins, 19 70) , 
Several time series models are available and can be used 
for representing the temporal variations of traffic flow. 
Different categories of vehicles generally have different 
temporal variations. It is proposed to develop time series 
models to represent monthly, daily and hourly variations 
of different categories of vehicles in the traffic flow. 
These models may be used for forecasting future traffic 


demand 



s 


Traffic generally is not homogeneous/ but consists 
of different categories of vehicles. In developing 
countries like India, it includes slow moving (man driven 
and animal driven) and fast moving vehicles among which 
there are wide variations in speed. It however seems 
necessary to represent the characteristics of each type of 
vehicle in terms of a common unit preferably the passenger 
car. They are referred to as Passenger Car Equivalents (PCES) . 
The mixed traffic flow process is complicated because of 
some of the following factors: (i) It is a complex stochastic 
process; (ii) Free speeds of vehicles are not constant 
but are probabilistic; and (iii) There are interactions 
between vehicles moving in the same direction and also with 
vehicles moving in the opposite direction. These include 
acceleration/ deceleration/ overtaking etc. 

In India, PCES are generally considered to be 
constants and to be independent of composition and volume 
level of traffic. Because of the complex nature of the 
process/ it seems reasonable to assume that interactions are 
nonlinear and hence PCES are functions of volume and compo- 
sition of traffic. It is necessary to investigate the 
interaction between slow and fast moving vehicles in the 
mixed traffic flow. It seems difficult, if not irrpossible. 
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to formulate and analytically solve for the interactions 
particularly because of the complex logical and probabilistic 
characteristics of the process. It is however possible to 
simulate the process in a digital computer, study the 
process and infer about the complex interactions. 

The problem may hence be stated as follows: 

(i) To represent monthly, daily and hourly 
Variations of traffic using time series models; and 

(ii) To study the interactions between slow and fast 
moving vehicles in the mixed traffic flow using computer 
sirrulation and hence to derive PCES of different categories 
of vehicles as a function of traffic composition and volume 
1 cvel . 

1.3 Objectives of the study 

The following are the main objectives of the study: 

(i) To develop stochastic models for the monthly 
and daily goods traffic on five highways touching Kanpur and 
to use these models for traffic forecasts; 

(ii) Data concerning characterisation of traffic 
including probabilistic distribution of flow, free speed, 
operating speed and delays etc., are not available. Field 
studies were conducted to formulate appropriate component 
model s; 


r 
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(iii) To formulate and validate a simulation modul 
for the mixed traffic flow on a two lane highway; 

(iv) To sinulate the mixed traffic process under 
Varying volume levels and compositions to study the interacti- 
ons; and 

(v) To derive the PCES for each category of vehicle 
as a function of traffic volume and composition. 

1.4 Scope of the Study 

The subject matter under study is very broad. Further- 
more, there are limitations of availability of data and 
computer time. So the present study was restricted to the 
f ol lowing: 

(i) Univariate time series models are developed for 
the monthly traffic data over the years and daily traffic 
data over the months. They are limited to the considerations 
of three categories of goods carriers approaching Kanpur on 
five different roads. 

(ii) Simulation of the mixed traffic flow was 
carried out only for a two lane highway. 

(iii) simulation was done only for straight sections 
and effect of various geometries like grade, curves and 


lateral obstructions/ etc.. 


were not studied 
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(iv) Free speed distribution parameters of 
different categories of vehicles used in this study were as 
estimated from field observations on one road only. 

(v) Traffic volume was taken to bu equal in 
each direction. The effect of different directional distri- 
butions on speed volume relationships could not be studied. 

1.5 Organisation of the Report 

The study is reported in the following sequence: 

(i) problems of mixed traffic phenomenon are 
reviewed. Traffic field studies and analysis of data 
for characterisation of mixed traffic flow are presented 
(Chapter 2) . 

(ii) Time series models are presented and used to 
represent monthly, daily and hourly traffic variations. 

Time series models of monthly flows are also used for 
forecasting (Chapter 3) . 

(iii) problems in mathematical modelling of mixed 
traffic flow are considered. Because of limitations of 
analytic models, the necessity of computer simulation of the 
process is considered (Chapter 4) . 

(iv) Computer simulation model for the process is 
formulated and validated using available field data 
(Chapter 5) . 
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(v) Simulation experiment s for different conpo- 
sitions and volume levels are planned. The sirailation 
results are analysed to evaluate the interactions between 
the vd^icles. The complex interactions are represented 
in terms of PCEs for any catagory of vehicle as a function 
of traffic volume and conposition . A nonlinear conditional 
multiplicative model is derived to represent the interactions 
( Chapter 6) . 

(vi) Study is summarised, conclusions are drawn 
and suggestions are made for future investigations (Chapter 7) . 
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2, CHARACTERISATION of mixed traffic flow 


2.1 General 

The ability of a roadway to accommodate traffic 
depends to a large extent on the physical features of the 
roadway. However, there are other factors not directly 
related to the roadway features which are also of major 
importance in determining the capacity. Most of these 
factors relate to the variations in traffic demand and 
the interaction of vehicles in the traffic stream (Pignataro, 
19 73) . 


2.1.1 Stream Characteristics 

Three basic measures of stream characteristics 
are volume, speed and density as they together describe 
the quality of service. Speed-volume relationship for 
uniterrupted flovj is shown in Fig. 2.1. As volume increases, 
the space mean speed of traffic decreases until a critical 
density is reached. Thereafter, because of jamming, both 
volume and speed decrease (Wohl and Martin, 19 67). Fig, 2,2 
shows that density increases with volume upto the point 
of critical density. Thereafter volume decreases as density 
continues to increase to a maximum value known as jam 
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density. Speed- density relationship in the upper range 
is similar to speed-volume relationship in that speed 
decreases with increasing density. Beyond the point of 
critical density, density continues to increase whereas the 
volume decreases, 

2,1,2 Variations of Traffic Flow 

Traffic flow is a stochastic variable characterised 
by the peak hourly rate. Traffic volume variations within 
the peak hour may seriously affect the stream characteristics. 
Capacity is normally defined in terms of vehicles per hour. 

If all the vehicles are equally spaced, determination of 
operating characteristics is simple. However, the vehicles 
do not move at uniform headways; rather they tend to form 
groups. Thus the rate of traffic flow for intervals of 
less than an hour can substantially exceed the peak hour 
rate, a peak hour factor that relates the peak interval 
and the peak hour volume, may be defined empirically 
(HRB, 19 65) . When demand exceeds capacity, congestion 
starts bail ding up within the peak interval or extend over 
a longer period depending upon the traffic flow subsequent 
to the peak rate. Hence it is necessary either to provide 
excess capacity over the full hour to accomodate peak intervals 
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of flow or to permit limited congestion . For a rational 
analysis of the stream characteristics/ the actual operating 
headway distribution is needed. 

2,1.3 Passenger Gar Equivalent s(PCES) 

PCEs of vehicles, as specified by Central Road Research 
Institute (CRRI) and Ministry of Transport (MOT) of India, 
are given m Table 2.1 ( Sehgal , 1967) . They show a wide 
difference in values. Furthermore, these PCEs are constant 
and are considered to be independent of traffic volume and 
composition. Highway traffic at low volumes does not genera- 
lly influence the speeds, and vehicles continue to move 
essentially at their free speeds. But as the volume 
increases, it is no longer possible for the vehicles to move 
at free speeds. When a fast moving vehicle is behind a 
slow moving vehicle immediately ahead of it, the fast moving 
vehicle desires to overtake. For a two lane road, the 
overtaking vehicle has to travel for some distance in the 
wrong lane (right lane in India) till the overtaking 
operation is conpleted. Overtaking is possible only if Ro 
vehicle coming from opposite direction comes in conflict 
with the overtaking vehicle. If overtaking cannot be carried 
out, then the overtaking vehicle has to reduce its speed, 
resulting in considerable delay to the fast moving vehicle. 
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TABLE 2.1 passenger car equivalents (PCES) for STRAIGHT 
SECTIONS OF ROAD 


Type of Vehicle Passenger Car Equivalents (PCES) 

Central Road Research Ministry of 

Instituted CRRI) Tran sport (MOT) 


Car, Jeep, Van 

1.0 

1.0 

Bus 

3.6 

3.0 

Truck 

2.7 

3.0 

4 Seater Rickshaw 

0.7 

- 

2 Seater Rickshaw 

0 . 6 


Motor Cycle 

0.3 

1.0 

Scooter, AU tocycle 

0.2 

1.0 

Bicycle 

0.4 

1.0 

Tonga 

2.6 

6.0 

Bullock Cart 

10.7 

6.0 

Hand Cart 

2.7 

6.0 

Cycle Rickshaw 

1.4 

1.0 
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With increases in volume level and proportion of slow 
moving vehicles in the mix, the number of gaps available 
in the opposite direction to allow overtaking decreases, 
resulting in considerable restriction on the overtaking 
operations. The operating speeds of the fast moving 
vehicles thus considerably get reduced. The flow charac- 
teristics of the mixed traffic flow are quite conplex due 
to the interactions between the vehicles. Because of these 
complex interactions, PCEs of different categories of vehicles 
are considered to vary with traffic volume and composition, 

2.2 Traffic Field Studies 

2.2.1 Necessity 

For characterising the mixed traffic flow, it is 
necessary to identify the probabilistic characteristics 
of the flow process and the effect of interactions between 
the vehicles. The following information is needed for 
system modelling ( Naylor et al 19 66) : 

(i) Character! sties of the input like traffic 
volume, composition, gap distribution ^ free speeds/ etc., 

(ii) Characteristics of the system like physical 
features of the roadway and other control measures, if 
any; and 
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(iii) Characteristics of the output like density, 
travel time, delay time , operating speed etc,, of 
various vehicles under different input conditions. 

The accuracy of the system model depends upon 
the data available. Information concerning the above 
characten sties was not available from records and so it 
was necessary to conduct field studies. 

2.2.2 Site Selection 

Grand Trunk (G.T,) Road -NH 2 , is a major link 
between the north and east of India with a considerable 
amount of traffic, a 2.25 km long section of this road, 
between Kalyanpur and Rawatpur, was selected for field 
studies. This section lies within Kanpur Metropolitan 
limits and thus has considerable amount of local fast and 
slow moving traffic alongwith the inter-city traffic. It 
had the following advantages with respect to the objectives 
of the study: 

(i) The selected section has 6.0 metre wide 
cement concrete pavement with about 1.5 metre wide unpaved 
shoulders on either side. The pavement is in good condition 
for free flow; 

(ii) The section lies on straight reach with no 
restrictions on sight distance. Thus the effect of aurves. 
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sight distance etc. on operating speeds does not come into 
play. 

(iii) There is no intersection on the selected 
section and so there is no cross movement affecting the 
flow of vehicles on the road; 

(iv) The road profile is almost flat with the 
maximum gradient being even less than one percent. Running 
speeds of vehicles are thus not affected by gradients; and 

(v) There is a railway track on one side of road 
and very little man made construction on other side. There 
is/ therefore, little possibility of pedestrians crossing 
the road and delaying the vehicular movement. 

The section is hence considered ideal for estimation 
of interaction effects and hence chosen for field studies. 

2.2.3 siting of posts 

Five observation posts were set up for data 
collection as shown in Fig. 2.3. Traffic volume of 
different categories of vehicles was recorded at 5 minute 
intervals at posts 1 and 5 respectively for incoming and 
outgoing traffic. Data for free speed of vehicles were 
recorded at posts 1, 2 and 3. Information concerning travel 
time, delay time , gap distributions/ etc., were 
recorded at posts 1, 3, 4 and 5. 
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2.2.4 Duration of Survey 

Traffic survey was conducted round the clock 
for ten days in June, 1973, Twenty five engineering 
students were trained for proper collection of data* The 
study was conducted in four shifts per day* Traffic 
volume was recorded for all the four shifts. Free speed, 
gap distribution and travel time data were recorded for 
one shift on each day so as to cover all the four shifts 
twice during the period of study. On the whole, infor- 
mation for about 7o,000 vehicles was recorded. 
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2.2,5 Data Collection 

Stop watches were synchronised for each shift 
before starting the data collection, procedures concerning 
tho recording of data are briefly discussed. 

(i) Traffic Volume Data: Incoming and outgoing 
volumes of different categories of goods and passenger 
carriers were recorded at 5 minute intervals. Volume of 
bicycle traffic was very high and this was recorded at 
30 second intervals. For goods carriers# the volume of 
loaded and empty vehicles were separately recorded, 

(ii) Free Speed Data:- Three posts namely 1# 2 and 
3 were set up with the distance between successive posts 
as 0.25 km. ( Fig, 2.3) , Free speed information was recorded 
by noting the time of crossing of a mark on the road by 
each vehicle at the three posts. In order to identify the 
vehicle, additional information including its category, vi z. , 
car# truck# bullock cart, bicycle , etc. ? make , colour , 
etc., were also noted. The positions of the observers were 
such that they could easily see the movement of any vehicle 
over a 0.25 km section on either side of the posts. So the 
observers also noted whether the vehicle is not accelerating/ 
decelerating or overtaking and hence is moving freely during 
the entire section of 0.5 km* If it is observed that there is 
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any restriction to free flow of a vehicle, information 
concerning that vehicle was crossed out from the records. 

If information about a vehicle was crossed even by one 
observer, then that vehicle was not considered for the free 
speed analysis. Scanning of vehicles was done from the 
type and identification recorded at each of the three posts. 
Travel time for free flow over 0.5 km distance was confuted 
and free speed for each vehicle estimated. 

(iii) Travel Time and Operating Speed Data: When a 
vehicle moving from either direction crossed a post ( 1, 3, 4 or 
5 ) , information concerning the type of vehicle, its make, 

colour , time of crossing, etc. , was recorded, separate 
observers were stationed for the two directions. Knowing 
the time of entry into-the section at post 1 ( or 5), travel 
times were computed from observation records at posts 3, 4 
and exit post 5 ( or 1) , The data were also used for 
determination of the following:- (i) operating speed of 
vehicles? (ii) delay time of vehicles? (iii) density of 
different categories of vehicles in the section at any 
instant? and (iv) interarrival time gaps from either 


di recti on 
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2,3 Frequency Analysis of Field Data 

2.3.1 General 

Field data were analysed to obtain the necessary 
inputs and outputs of the system model. Analysis was done 
for the following characteristics: 

(i) Traffic volume distributions of various categories 
of vehicles at different time periods. This is discussed 
in sec, 3,6, 

(ii) Interarrival time gap distributions of the vehicles 
on the road. 

(iii) Free speed distributions of different categories 
of vehicles. 

(iv) Travel time, delay time and operating speed 
distributions of various vehicles at different volume levels 
atid traffic composition. These deal with the output of the 
system and are useful for validating the system model 
(Sec. 5.10) , 

2.3.2 Interarrival Time Gsps 

Vehicles arrive at the section from either direction, 
apparently without any regularity. The flow may be considered 
as random, i,e,, the number of vehicles arriving in any 
interval of time follow a probability distribution and is 
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independent of the number of vehicles that arrived during 
any previous time interval* From the arrival time data/ 
interarrival time gaps were confuted. The frequency dis- 
tribution of arrival gaps is shown in Fig. 2. 4 for one shift 
of a day. Fig, 2,5 shows the cumulative distribution of 
arrival gaps. Assuming that the vehicle arrivals are 
randomly distributed, Poisson distribution may be used for 
the number of vehicles arriving in unit time and in turn 
exponential distribution may be used for the interarrival 
time gaps ( Gerlough, 1964 ). Vehicles moving on a highway 
in the same direction require at least some physical 
spacing and hence a minimum time gap. Fig, 2.4 and 2.5 
show that the minimum observed headway, h m ^ n Was 1 second. 

By shifting the exponential curve from the origin along the 
x-axis by an amount h m j_ n = l second, the minimum headway, 
the following shifted exponential distribution -was fitted 
to the intcrarrival time gaps, 

~ h min ^ 

_ V 

prob ( h < t) =[ 1-e h min ] far t^/h^n (2,1) 

^ ^rrvin ^ 4 

Freq ( h< t) h^ f ] for t^h^ (2.2) 

where h = interarrival time gap in seconds; C = T/V = average 
interarrival time gap in seconds; V = total number of vehicles 
arriving during time T; and h min = i second. 
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The theoretical frequencies of arrival gaps as 
obtained from shifted exponential distribution are also 
shown in Pigs. 2.4 and 2.5 alongvnth the observed frequencies. 
Goodness of fit was tested by Chi-square test. Results 
confirm that shifted exponential distribution fits well 
with the interarrival time data in other cases also. 

2.3.3 Free Speed Distributions 

The free speed of vehicles were computed from the 
times to travel 0.5 km, under free flow conditions. The 
vehicles were categorised into the following six types, 
each having similar characteristics; 

Ci) Heavy motor vehicles like trucks and buses; 

(ii) Medium motor vehicles like passenger cars, 
jeeps, tempos and taxis etc.; 

(iii) Light motor vehicles like scooters and 
motorcycles? 

(iv) Horse driven vehicles like tongas or 
kharkhara? 

(v) Bullock carts; and 

(vi) Bicycles and cycle rickshaws* 

Free speed data were recorded for about 2300 vehicles. 


They include all motor vehicles, tongas and bullock carts 
moving at free speed during the period of study and only 
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around 4oo bicycles whose volume was much higher. The 
free speed of vehicles of a particular category were found 
to vary. The variations may depend upon the make of vehicle, 
ics age , driver's characteristics, etc. The frequency 
of free speed distributions of 480 heavy motor vehicles is 
shown in Fig, 2,6, It has a peak in the middle with low 
values at very high and very low speeds. The observed mean 
was 44,6 kirph with a standard deviation of 6,65 knph, A 
normal distribution was fitted to the field data and parameters 
were estimated. 

For estimating the parameters, it is necessary to 
determine the p robability density function f(x) or the 
cumulative distribution function F(x) from the observed 
data, F(x) was determined by the plotting position formula 
(Benjamin and Cornell, 1970) , 

p(x) = — — (2,3) 

N + 1 

where p(x) is the exceedence probability P(x) = 1 - F(x) ; 
m = rank of the observation with the highest value having 
m a 1/ and N = total number of observations. 

Plotting positions P(x) vs, x were plotted on 
normal probability paper and a straight line was fitted 
by eye sight. The fitted straight line represents a 
normal distribution with a mean of 44,3 knph and a standard 
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deviation of 6.8 krrph, Chi-square test -indicated the 
goodness of fit at 10 percent significance level. Normal 
distribution was fitted to free speed of other categories 
of vehicles also and the estimated parameters are shown 
in Table 2.2. The maximum deviation of speeds from mean 
was found to be not more than thrice the standard 
deviation ( sd) and so the normal distributions were censored 
at ( x ± 3 sd ) limits. 

The results of frequency analysis are used ip 
the validation of simulation model and in the simulation 
of mixed vehicular traffic. 
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table 2.2 parameters of normal distributions fitted to 

FREE SPEED OF VEHICLES 


S.No. 

Type of Vehicle 

Sample 

Size 

Mean 

Free 

Speed 

kmph 

Standard 
Deviation of 
Free Speed 
kmph 

1. 

Heavy Motor Vehicles 
ie. Trucks and Buses 

480 

44.30 

6.8 

2. 

Medium Motor Vehicles 
ie. Cars , Jeeps/ Tempos etc* 

585 

49.95 

8.05 

3 • 

Light Motor Vehicles 
ie, scooters, Motocycles 
etc. 

298 

39.80 

7.6 

4. 

Horse Driven Vehicles 
i e. Tonga, Kharkhara 

354 

12.8 

1.5 

5. 

Bui lock Carts 

82 

5.2 

1.1 

6. 

E£cycles and Cycle 
Rickshaws 

512 

11.9 

1.8 
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3. STOCHASTIC ANALYSIS OF TRAFFIC FLOW 

3.1 introduction 

Traffic volume varies within the year with seasonal 
weekly and. daily variations and they are to be considered 
in traffic forecast. The system should be designed for the 
peak traffic in a suitable period, say an hour, perhaps 
with a certain probability of meeting the demand. The 
design should be based on stochastic modelling and analysis 
of historical data for a sufficiently long period. 

3.2 Stochastic Processes 
3.2,1 General 

When the future variations cah be predicted uniquely on 
the basis of the present and past values of the characteristics/ 
a process may be considered to be a deterministic process. 
Actually the variations of traffic volume are erratic and 
future volumes cannot be uniquely predicted, predictions are 
then to be made only by using probabilistic models and in 
terms of the chance of occurrence of an event, Furthermore, 
the sequence of occurrence of events is important in the 
analysis of traffic flow. Mathematical models may be used 
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to describe the process and to analyse its characteristics. 
The theory of stochastic processes deals in analytic terms 
with phenomena which develop in time in accordance with 
laws of probability ( Box and Jenkins, 19 70) . 

3.2.2 Time Series 

A time series is an ordered set of realisations 
of the process. Let Z (t) be a characteristic defined in 
the time domain t 8 t ie,, [ z(t) ; t &T ] . The chara- 

cteristic may be defined as a continuous function of time 
and then it is referred to as a continuous time series 
( Ramaseshan, 19 72) . However, traffic volume data are 
given for discrete values of time as sampled or quantised 
data. The time series are then referred to as discrete 
time series. Traffic volume observations may be made at 
constant intervals of it; say ( t 0 + At ) , ( b 0 + 2 At ) , 

( t Q 3 At) ,•••■( ( t Q t n At) . Then the N successive 
observations Z^, Z 2 r Z 3 ,..., Zj^ respectively constitute 
a discrete time series. 

3.2.3 Components of a Time Series 

A historical time series record may be considered 
to consist of ( Kisiel, 19 69) : 
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(i) Deterministic components like trend and 
cycles which are purely functions of time; 

(ii) Other nonrandom components like autoregressive (ah) 
and moving average (MA) components that represent 

the internal and external correlation respectively; 

(iii) Pure random component; and 

(iv) Catastrophic events, which are extremely rare 
events that have occured during the limited period of 
record. 

The time series may thus be represented as; 

Z(t) = f^(t) + f 2 (t) + e t + m6 (t - t c ) (3.1) 

where f-j_ = deterministic components like trend , cycles 
etc.; fg= persistence conponent; e = pure random corrponent; 
and in = magnitude of catastrophic event occuring at time 
t c . Catastrophic events are generally ignored in the repre- 
sentation of time series models. 

A time series is said to be stationary if the 
generating mechanism of the process is independent of time 
(Yevjevich , 19 7 2) i.e. the deterministic conponent and 
the parameters of the random conponent (and its distribution) 
are all independent of time. A stochastic process, z(t), 
is said to be stationary in a wide sense if it possesses 
finite second moment and covariance of Z(t) and z(t + h) 
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depends only on h for all t g T. Nonstationary time 
series thus need to be transformed to stationary series 
for estimating the persistence components. 

3.2.4 Stages in the Selection of a Model 

The selection of a stochastic model from historical 
records requires an iterative approach and the steps 
generally used are indicated in Fig. 3.1 (Box and Jenkins, 19 7o) , 
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(i) A useful general class of models is considered 
from the int© raction of theory and practice. 

(ii) Rough methods for identifying subclass of the 
models are developed. The identification process can be 
used to yield rough preliminary estimates of the parameters 
in the model. 

(iii) The tentative model is fitted to the data 
and its parameters are estimated. The rough estimates 
obtained daring the identification stage may be used as 
starting values in an iterative method for estimating the 
parameters, 

(iv) Diagnostic checks are applied with the object 
of uncovering possible lack of fit and diagnosing the 
cause, if any inadequacy is found, the iterative cycle 
of identification, estimation and diagnostic checking is 
repeated till a suitable representation is found, 

3,3 Traffic Studies 

Estimation of traffic demand on highways 
needs the following variations to be evaluated: 

(i) .Variations within the year and the month 
of the volume of different categories of vehicles using 
the road; and 



( ii) Variations within the day to account for 


hourly volume characteristics. 

For case (i) goods traffic data were availabje for 
five main roads approaching Kanpur (Subsec. 3,3.2) whereas 
for case (ii) field studies were conducted for data 
collection (Sec. 2 . 2 ) . 

3,3.1 Choice of Region 

Kanpur is the largest industrial city in the 
state of Uttar Pradesh (U.P.) . The population of Kanpur 
Metropolis/ as per 1971 census/ is 12.65 lakhs and it 
enjoys the eighth position among the metropolitan cities 
of India (KDA, 1975) , It is served by two National 
Highways/ G.T. Road (NH 2 ) and Lucknow-Kalpi (NH 25 ) Road. 

Thus it is connected to Delhi, Allahabad, Lucknow and Jhansi 
and thence to all the major cities of India, Hamirpur 
Road, a State Highway connects Kanpur to the districts 
of southern U.P. It is also connected to all the major 
cities of India by the broad gauge Central and Northern 
Railway lines and metre gauge North-Eastern Railway .lines 
(KTC, 19 75). Fig. 3.2 shows the transportation network in 
Kanpur Metropolis. 
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3,3,2 Classification of Data 

There is considerable road traffic between Kanpur 
and other adjoining cities due to various industries and 
commercial activities located at Kanpur. The agricultural 
produce of adjoining areas is also brought to Kanpur mostly 
by bullock carts. Further/ the direct traffic between Delhi 
and Calcutta, and between Western U.P. and Bombay side, also 
pass through Kahpur on the two National Highways. 

Delhi Road, Allahabad Road, Lucknow Road, Kalpi Road 
and Hamirpur Road are the five main road sections connecting 
Kanpur, Historical traffic data on these five roads were 
taken from the octroi records of Kanpur Municipal Corporation. 
The data were available only for goods traffic approaching 
Kanpur, which has to pay octroi. Daily data, recorded in shifts 
of eight hour duration, were available for seven years 
(April, 1966-March, 1973) for motor vehicular traffic. 

Bullock cart traffic data were available for three years. 

(April, 19 7o-March, 19 73) . The available data were classified 
into three categories, viz,, (i) transit truck traffic 
i.e., commercial vehicular traffic touching Kanpur, but 
destined to some other place beyond Kanpur; (ii) incoming 
truck traffic l.e.,, commercial vehicular traffic destined 
for Kanpur; and (iii) incoming bullock cart traffic that 
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comes from adjoining villages to Kanpur, 

Traffic volume studies were conducted on G.T. Road 
(Sec, 2, 2 ) so as to analyse the hourly variations of 
all categories of vehicles in both directions, 

3,3,3 .Analysis of Data 

Historical daily data of goods traffic and hourly 
data from field studies were used to develop and validate 
appropriate stochastic models. They include! 

(i) stochastic models for monthly traffic flow 
for each of the three categories of vehicles on all the 
five roads? 

(ii) stochastic models for daily traffic within 
a month in two roads for three representative months in a 
year; and 

(iv) distribution of hourly traffic over a day 
for different categories of vehicles on G.T, Road, 

3,4 Analysis of Monthly Traffic Flows 

3.4,1 Introduction 

Monthly traffic volumes# represented by the mean 


daily traffic of the month, were plotted for all the series, 
viz., incoming truck and bullock cart traffic on Hamirpur 
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Road and transit traffic also in case of other four roads. 

Data for transit truck traffic on Delhi Road (19 66-19 73) 
are shown in Pig, 3,3. The series were analysed using 
two different stochastic models/ viz., the general time 
series model(Eq. 3.1), and a seasonal multiplicative 
Autoregressive Integrated Moving Average ( ARIMA ) model to 
be discussed later in Subsec. 3.4,3. 

3,4,2 General Time Series Model 

Trend Analysis: It deals with the identification 
and evaluation of long term trends, if any, in the time 
series. Fig. 3.3 shows that transit traffic on Delhi road 
increases generally with time, indicating the presence of 
periodic and random fluctuations. In some of the other 
series, annual traffic has more than doubled over a period 
of seven years, indicating a significant trend component. 

Due to the presence of periodic and random fluctuations, 
it is difficult to identify the actual nature of trend and 
its parameters. Method of moving averages (Kendall and 
Stuart, 1966) was used for fitting the trend line. Let the 
weighted average , 


bj Z ( t + j ) 


Z‘ (t) 


j =m/2 
Z 

j=-m/2 


( 3.2 ) 
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with j 


m 



m 


2 


+ 1 


1 , 0 , 1 , 2 


m , m 

— — ** ,1, — - ; and rhe 

2 2 

1 

weights b s = - 

J m + 1 

For the monthly traffic flows, 12 monthly moving 
average was used as it eliminates the annual cycle or its 
subharmonics , if any* The nature of z' (t) series was 
mostly exponential though linear trend was also observed 
for two series. The relationships were of the following 
types; 


Z 1 (t) 

- A . 

1 


z' (t) 

= At B 

1 

(3.3) 

z‘ (t) 

— A + Bt 1 



The regression coefficients A and B were estimated 
by the method of least squares for various time series and 
are shown in Table 3.1. The results indicate that generally 
bullock- cart traffic increases at a faster rate than the 
motor vehicular traffic. The trend free series Z^(t) were 
obtained after subtracting the trend component from the 
original values. 


Periodicity Analysis: Seasonal fluctuations of 

traffic volumes are indicated by the plots of the time 
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TABLE 3.1 REGRESSION ANALYSIS OF TREND FREE SERIES 


Road 

Category 

Trend 

Regression 

Coef fxci ents 


of Traffic 

Rel a ti on ship 

A 

B 

Delhi 

Road 

Transit Truck 

Z' (t) 

= Ae Bt 
= At B 

60.0 

0.012 


Incoming Truck 

z 1 (t) 

88.0 

0,170 


Bullock Cart 

Z' (t) 

= A + Bt 

3o«0 

2.500 

Allahabad 

Road 

Transit Truck 

Z'(t) 

Bt 

= AC 

61.0 

0.0 20 


Incoming Truck 

Z' (t) 

= A + Bt 

75.0 

0.210 


Bullock Cart 

Z' (t) 

= At* 

36.0 

0.350 

La cknow 

Transit Truck 

z l (t) 

= At 5 

42.6 

0.0593 

Road 

Incoming Truck 

z' (t) 

a At B 

138.0 

0.0782 


Bullock Cart 

z' (t) 

= 7ie Bt 

149.0 

0.01543 

Kalpi 

Transit Truck 

z' (t) 

= Ae Bt 

25.5 

0.0207 

Road 

incoming Truck 

z' (t) 

*= Ae Bt 

68.0 

0 .017 


Bullock Cark 

z' (t) 

= Ab B 

34.0 

0. 430 

Hami rpu r 

Incoming Truck 

No significant 

0.0 

0,0 

Road 


Trend 





Bullock Cark 

z' (t) 

= At B 

30.0 

0.460 
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series. Hence the monthly mean of trend free traffic 
volumes wore calculated and are shown in Figs. 3.4 and 3.5 
respectively for incoming truck and bullock cart traffic. 
They indicate the seasonal variations of the means with 
a peak volume around January and the minimum around August- 
September for truck traffic. For bullock cart traffic, 
there are generally two peaks, the larger one around 
January and the other around June. Hence seasonal cycles 
seem to be present. It may be possible to relate the 
seasonal variations to agricultural operations and monsoon 
season. 

The trend free monthly traffic flow series were 
found to have a skewed distribution, A log normal distri- 
bution was fitted to the monthly data. Let X^(t) denote 
the natural logarithm of the trend free series Z^Ct) ie., 
Xi<t) = Log e Zi(t) , The periodic component of a time 
series can be represented by a sum of a number of harmonic 


components ( Kisiel, 1969) ie.. 


Xi(t) = 


2 ft jt 

Cj Co s ( — ^ + 


(3.4) 



m 

2 A.j Cos 
j=l J 


2 71 jt 


m 

T. B. sin 

1=1 


2 71 j t 


or 


Xi(t) 


2 


T 


T 

(3.5) 
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3o 

where — — is the mean of X-j, (t) ; 
period; fundamental frequency 
the Fourier coof ficients? Cj and 
phase angle respectively of the j 


T is the fundamental 

1 , r-v 

= ; Aj and Bj are 

a are the amplitude and 
th , 

harmonic. 



2 N 2 uj t 

— S X l< fc > ^ —y 


B 


j 



N 

S X. ( t) sin 
t=l 1 


2 itj t 
T 


2 


G j 



+ 



(3.6) 

(3.7) 

(3.8) 


and a . 

J 


arc tan 


( - 



) 


(3.9) 


The periodic and persistence components may be 
identified through (i) correlogram (Box and Jenkins, 1970) 
which is a plot of serial correlation coefficients of the 
time senes for different lags; (ii) partial autocorrelation 
coefficients; and (ni) power spectra which indicates the 
distribution of variance of the process among the different 
frequencies (Jenkins and Watts, 19 68) . Standard procedures 
were used in the estimation of the above functions. 

Hie correlogram and spectra were computed for all 
the trend free traffic series, some of the results are 
shown in Figs* 3,6 and 3,7 alongwith the 95 percent 
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confidence level for a pure random series. Correlogram 
and pov/er spectra analyses indicate the presence of a 
12 monthly cycle for all categories of traffic ( vide 
spikes at f = 0.0833 in Figs. 3.6 and 3.7 ) . For the 
bullock cart traffic, 6 monthly cycle is also indicated 
(vide spike at f = 0.1667 in Fig. 3.7 ) . The harmonic 
coefficients and ct ^ for the 12 monthly cycle 

were estimated for each of the series and were eliminated 
from the X^(t) series. The residual series were again 
subjected to correlogram and spectral analyses. Six 
monthly cycles were found to be significant for the 
bullock cart traffic only. The harmonic coefficients 
B g/ Cg and a 2 for the 6 monthly cycle were also estimated. 
Fig. 3,8 shows the X^(t) series with the 12 monthly cycle 
superimposed; the residual series ( after elimination of 
12 monthly cycle ) with 6 monthly cycle superimposed, 
and the residual series after elimination of 6 monthly 
cycle as well. The parameters of the harmonic coefficients 
of various series are shown in Table 3,2. 

persistence Components: Persistence refers to the 
linkage or relationship existing between the value of 
a series at a given time with the earlier values/ viz./ 



6 


12 


24 


30 


36 



18 

MONTHS 

HARMONIC COMPONENTS AND RESIDUALS OF 
TRANSFORMED SERIES 


51 


TABLE 3,2 HARMON I C COEFFICIENTS OF PERIODIC CYCLES 


Road 

Category of 

Annual 

Cycle 

Six Monthly Cycle 


Traffic 

*1 

B 

1 

*2 

B 2 

Delhi 

Transit Truck 

0.164 

-0.034 


- Ate 

Road 

Incoming Truck 

0. 273 

-0.136 

- 

- 


Bullock Cart 

0.462 

-0. 207 

-0. 219 

-0.0475 

Allahabad 

Transit Truck 

0. 133 

-0.027 



Road 

I n comi ng Tru ck 

0.098 

-O ■ 0 44 

- 

- 


Bullock Cart 

0.275 

-0. 121 

-0.145 

-0 *03 2 

La cknow 

Transit Truck 

0. 145 

-0.036 


_ 

Road 

Incoming Truck 

0.174 

-0.094 

- 

- 


Bullock Cart 

0.448 

-0. 214 

-0.271 

“0 *0 53 

Kalpi 

Transit Truck 

0.09 2 

-0.00 27 


.. 

Road 

Incoming Truck 

0.079 

-0.032 

- 

- 


Bullock Cart 

0.403 

-O. 124 

-0.224 

-0.083 

Hami rpu r 

Incoming Truck 

0.251 

-O. 113 



Road 

Bullock Cart 

0. 264 

-0.107 

-0.125 

-0.034 


I « - i . ,/K 

CENIRaL HORARY 


Acc* No. 


A 
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Z(t) = f' 2 fz(t-i) , Z(t-2), Z(t-k) , t ] 

+ £"[y ( t-1) , Y ( t~ 2) , . . . . Y( t-m) , t J (3.10) 

where f^ is due to the internal dependence and ££ is due 
to external dependence. 

Due to the presence of seasonal cycles, the time 
series are nonstationary. It is necessary to transform 
them into stationary series before fitting the persistence 
models ( Thomas and Fiering, 1962? Mckerchar and Dell tur, 1972) . 
The following transformation was used in this study; 

X (t) - X,(j) 

X 9 (t) = — — (3.11) 

2 S (j) 

where X-^(j) and s(j) are the mean and standard deviation 
of the logarithms of the trend free scries Z-j_(t) in month j, 
j = 1, 2 , 3, 12, X 1 ( j) and Stj^eru estimated by 

maximum liklihood method, The transformed X£(t) scries 
having zero mean and unit variance is stationary in mean 
and variance. The standardised scries X 2 (t)were analysed 
for persistence. 

Model Identification; Before a time series model may 
be fitted to the data, it is necessary to postulate or 
identify theoretical forms of the model. Oorrelogram, partial 
autocorrelation function (pacf) and power spectra are used 
for this purpose. For any given model, its ccrrelogram, pacf 
and/or power spectra have specific characteristics. These 
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may be compared with the empirical values estimated from 
available data to identify a possible theoretical model 
that can be used for representing the time series. Persistence 
may be represented by .Autoregressive (AR) and/or Moving 
Average (MA) models. Only AR model is described in this 
section. MA model is described later in Subsec. 3.4,3. 

Ah AR process of order p ( Box and Jenkins, 19 70) 
may be represented by; 

X 2 (t) = X 2 (t-1) + 0 2 x 2 (t~2) + .... 

+ 0 p X 2 (t-p) + a t (3.12) 

a t is the random component and 
the other terms represent the persistence 

component, i.e,, the present value of X 2 (t) depends upon 
the proceeding p values* 

Let B = backward shift operator i.c. ( , B X 2 (t) = X 2 (t-'l) 

and ( 1-B) X 2 (t) = X 2 (t) -* X 2 (t,-i) , The AR operator of 


order p is defined by; 

0(B) = 1- B - 0 2 b 2 - - 0 p Ep (3.13) 

So the AR model may be written as? 

0 (B) X 2 (t) = a t (3.14) 


The model contains ( p + 2) unknown parameters, ie % , 

^1 , ^2/ ...... i P n / ft, and 0 (variance of the random 

series a t ) . 
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The autocorrelation function (acf) and pacf -were 
computed for each of the standardised scries X 2 (t) P They 
are shown in Pig. 3,9 for incoming truck traffic on Lucknow 
road alongwith the 95 percent confidence limits. The 
decaying exponential behaviour of the estimated correlogram 
as well as significant partial autocorrelation at lag one 
identifies a first order aR model to represent the time 
series X 2 (t) i.e,, 

X 2 (t) = f> x X 2 (t-l> + a t (3.15) 


Estimation of Parameters: It is necessary to estimate 

the values of the parameters of the identified model from 
available data. For an AR process, p coefficients can be 
initially estimated by Yule-Walker equations (Mckerchar 
and Delleur, 1972). Let r^_, r 2/|1 ,, etc, be the acf for 
lags 1, 2 etc. 


Then, for p = l , 

and for p = 2 , 



(3.16) 

(3.17) 

(3.18) 


Initial estimate of the residual variance cT e 
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the AR process may be obtained as follows (Box and Jenkins, 
19 70) : 
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cj a 2 - C Q ( 1 - 0 1 r-L - 0 2 r 2 - - 0 p r p ) 

(3. 19) 

x N 2 

where C = f (x 2 (i) - % ) (3.20) 

For the first order model identified for data 
series, parameter estimates are obtained as follov/s: 

H = 

° a 2 " c o ( 1 - 0i ^ ) = q, - 0 1 C x (3.22) 

Box and Jenkins ( 1965, 197o) suggest a nonlinear 
iterative algorithm that leads to the maximum livelihood 
estimation of the model parameters by minimising the variance 
of the residuals, viz., the residual sum of squares la^ « 

This procedure was adopted, with the estimates from Yule-Walker 
relations as initial values, and the parameters were 
determined. The initial values and successive estimates 
for bullock cart traffic on Delhi road are given in 
Table 3.3. Convergence is reasonably fast; and the initial 
values estimated are not far from the final values in this 
case. The final estimates for all the cases are listed in 
Table 3,4 , Residual variances lie between 0.498 and 0.816 
indicating that between 18,4 to 5o.2 percent of variance is 
explained by persistence. 
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TABLE 3,3 ITERATIVE ESTIMATION OF PARAMETERS 


(Delhi Road - Bullock cart traffic ) 


Iteration 

01 

Sum of Squares 

Starting value 

0 . 2200 

29.39 

1 

0. 1952 

24.01 

2 

0.1913 

23.98 

3 

0. 1887 

23.97 

Final value 

0.1887 

23,97 


Diagnost-id Checking: It rufL4rs"1rf> t -€he shocking 

and validation of the fitted model and if it is not 
adequate, the determination of inadequacies to indicate a 
better model. The techniques normally used are (i) overfit- 
ting and (ii) examining the residuals. 

Overfitting means fitting a more elaborate model 
than identified and comparing the fits in terms of the 
residual sum of squares of the more elaborate and less elabo- 
rates models thus giving an indication as to whether the model 
identified may be irrp roved upon or not* Examining the 
resi dials gives a more general check on the adequacy of 
the model. If the model is adequate and estimated parameters 
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T7.BLE 3,4 FlWZiL ESTIMATES OF PARAMETERS FOR FIRST ORDER 
7iR MODEL 


Road 

Category 
of Traffic 

0 2 
a 

*1 

Q 

(23 DOF ) 

Delhi 

Transit 

0.743 

0,3713 

25.82 

Road 


Incoming Truck 

0,624 

0.640 4 

30.73 


Bullock Cart 

0.778 

0.1887 

18.55 

7»11 ah ab ad 

Transit 

0,785 

0.28 76 

24.98 

Road 


Incoming Truck 

0.602 

0.4732 

8.86 


Bullock Cart 

0.516 

0.3179 

25.27 

Lucknow 

Transit 

0 . 730 

0. 1917 

35.21* 

Road 


Incoming Truck 

0. 59 2 

0.4276 

30.30 


Bullock Cart 

0. 758 

0.3938 

21.87 

Kalpi 

Transit 

0.728 

0.38 54 

11.57 

Road 


Incoming Truck 

0. 498 

0.5579 

22.85 


Bullock Cart 

0.617 

0.4445 

20.73 

Hamirpur 

Incoming Truck 

0. 528 

0 . 48 5 4 

33.38* 

Road 


Bullock Cart 

0.816 

0,5728 

20.26 

^23 ( 10 %) 

= 32.0 

^4 (5 c/o) 

= 35. 

20 


* Significant at 10 % level tout not at 5 % level 


59 


are close to the true values, the residuals should be a 
purely random series. If a-j- series are not purely random 
but shov; some serial dependence, the tentatively identified 
model is inadequate and it is to be modified. The nature 
of the dependence gives some indication of the modification 
that need be made. In the present study, diagnostic 
checking was done by examining the residuals. The two 
checks that were errployed are (i) autocorrelation of the 
residuals and (ii) ouiTUlative periodogram of the residuals. 

The rcsidial autocorrelation , , phould be 

independent and normally distributed about ?ero with 

variance l/N . For lower lags and for short 'series, the 

estimates of r L and r 2 are biased and it is not proper 

to use the standard error as a basis for determining 

\yhether the corrclogram represents a random series. Instead 

k 2 

a statistic Q = N 2 r, can be calculated for the first 

i=l 1 

k values of autocorrelations of the residuals, to indicate 
the randomness or otherwise of the residual series (Eurfoin, 
1969; Box and pierce, 1970) . For the truly randan series, 

Q is distributed as ^ with (N - p - q ) deg roes of 

freedom where p and q are the orders of AR and hA components 

respectively. 
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For seasonal time senes/ acf of residuals may 

riot be a sensicivu indicator of th>- departures from randomness, 

because periodic effects will typically dilute themselves 

among several autocorrelations. The periodogram can bo 

used for detecting the periodic pattern in the background 

of white noise ( Barlett, 1955) , The power spectra p (f) 

2 

for white noise has a constant value 2 cf over the 

cl 

frequency domain 0-0,5 cycles. The cumulative spectrum 

f 

for white noise p ( f) = / p (g) d g plotted against f 

° 2 
is a straight line running from ( o , 0 ) to ( 0.5, ) 

„ 2 

ie., P ( £ ) / °a is a straight lino running from (0,0) 

to ( o . 5 , 1 ) , J 

£ 1 < £ i> 

G (f.) = (3.23) 

j n (sd) 2 

2 2 

provides an estimate of P(fj) / a where (sd) is an 

2 J a 

estimate of o' and I (f .) an estimate of P (f .) . The 
a J J 

cuitulativo periodogram is the sample estimate of the 
integrated spectrum. On the basis of Kolmogorov- smimov 
test ( Hald, 1952), the confidence limits of the cumulative 
periodogram are given by j. K e /Yq where K e is the 
coefficient for calculating probability limits and 
q = ( N - 2) / 2 for N even and ( N - 1 ) / 2 for N odd. 
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The; autocorrelations for the residual a t series 

were computed. They vjere not significant at 95 percent 

confidence levels/ thereby indicating serial independence. 

Value's of statistic Q are shovm in Table 3.4 for all the 

series. These values were generally not significant compared 
2 

withjl^ values at lo percent level of significance. Only 
in two cases/ Q values were significant at 10 percent level 
but not at 5 percent level. The cumulative periodogratn 
for the residual a t scries is shown in Fig, 3.io alongwith 
5 and 25 percent probability limits. They also confirm that 
the residual series do not have serial dependence. The 
fit of the models is hence good and they are accepted, 

3.4,3 Seasonal ARIMA Model 

The autoregressive model has already been presented 
in Subsec, 3.4.2, The deviation of a process from mean, 
say, w t = ( z(t) - z ) / may be represented by a moving 
average process (Box and Jenkins, 1970 ) of order q, via,, 

w t = Z(t) - Z = a t - Q 1 a t _ tl - e 2 a t _ 2 - 

- 9 a. (3.24) 

q t-q 

where 0 , 0,, , 9 are a finite get of 

parameters and a t is a random genes with zero mean and 
2 

variance cs , Ihc process can also be represented in 
a 



ALLAHABAD ROAD 
INCOMING TRUCK TRAFFIC 



FREQUENCY, fj 


3-10 CUMULATIVE PERIODOGRAM OF 
RESIDUAL AUTOCORRELATIONS 
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t^rms of backward shift operator, B by 

2 a 

w, = Z(t) - Z = ( 1 - 0 B - 9 B - -0 B ) a. 

^ 12 q t 

= 0 ( B ) a fc (3.25) 

where 0 (B) is a polynorranal in B of order q . 

A more general model is the Autoregressive Moving 
Average ( ARMa ) model of order ( p , q ) with both AR 
and Ma components, viz., 

w t = Z(t) - Z - w t _ 1 + 0 2 w t _ 2 + .... + 0 P w t _ p 

+ a fc 0 i " 0 2 2 “* * ■ » • — 0 q 

q (3.26) 

or 0 (B) ( Z (t) - Z ) = 9 (B) a t (3.27) 

The series of traffic volume generally exhibit 

a trend along time, representing growth, and seasonal 

variations. When trend is exhibited, it may be possible 

to represent nonstationanty in terms of differences rather 

than in terms of the original time series. ARIMA process 

assumes that the dth degree of differencing of the series 

can be represented by an ARMA process of order ( p , q ) , 
d v d 

namely 7 = ( 1 - B ) and 

0p(B) (z(t) - Z ) = 6^ (B) a t (3,28) 

is referred to as an ARINA process ( Mckorchor and Dell cur, 1972) 
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of order ( p , d , q ) , 

When seasonal effects alone are present, they 

may be eliminated by usincj seasonal differencing; Let 

s = period of differencing ( s = 12 for monthly series data 

D s D 

with an annual cycle) ; V = { 1 - B ) ; where 

s 

D = degree of seasonal differencing. The time scries 
may thus be represented by ARIMA model ( Bacon, 19 65) of 
order ( P , D , Q ) g 

(3 s ) V° (Z(t) -Z) = o (B S ) e. (3.29) 

p s Q 


where 

4> p (B S ) = ( i _ * B*- $ 2 B 2S - .... - $ p B PS ) , 

a seasonal autoregressive operator of order P, and 

O ( B S ) = ( 1 - 0 B S - 0 B 23 - ~ 0 B QS > ' 

U 12 Q 

a seasonal moving average operator of order Q, 

When both trend and seasonal effects are present, 
the multiplicative seasonal ARIHA model of order 
( p , d, q) x ( P , D, Q ) _ can be represented by 

D a 

* p (B S ) 0 p (B) V g V ( Z (t) - Z) » 0 q (B ) 6 (B) a t 

(3.30) 

ARIMA models have the advantage of representing 
complex time series with only a relatively small number 
of parameters. 
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Model Identifications -Ltigarithnis. (: Tratilrai ) -of 
monthly traffic flows tend to be normally distributed and 
were used for farther analysis. Autocorrelations were 
estimated for the following cases 5 (i) undifferenced series/ 
i.e ,, original log transformed series; < i i) first differencing 
of the series w.r.t months only, i.e., d = 1 ; (iii) first 
seasonal differencing of the series,i,e„ D = 1 v;ith s = 12 
to take into account the annual cycle; and (Lv) first 
seasonal and serial differencing of the series d = 1 , 

D = 1. 

The estimated autocorrelations of differenced series 
for transit traffic on Lucknow road arc shown in Table 3.5, 

It is seen that for undifferenced senes as well as for 
simple differencing (d = l) and for serial and seasonal 
differencing ( d = 1 , Dal), serial correlations are 
not negligible for a number of lags and they fail to die 
out at higher lags. By contrast only seasonal differencing 
(D « 1 ) alone gives an improved result. The autocorrel ac- 
tions of this scries (9 = 1) are statistically significant 
2 

a t -j- confidence limits only for lag 12. The 

' N 

insignificant values between lags 2-11 and large values 
at lags 1, 12, 13, 14, 24, 25 and 26 indicate a multipli- 

cative seasonal ARIM/* process of order ( 1, 0 , 0 ) x ( 0 , 1 , 1)22 
i.e,„ p=Q = l,P=q=0, d =o,and D=l„ it is of the form 
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TABLE 3.5 ESTIMATED AD TO CORRELATION JUNCTIONS FOR DIFFERED CED 
SERIES ( TRANSIT TRAFFIC - LUCKNOW ROAD ) 


d 

D 

n 

Lags 

Estimated Autocorrelation 
Functions 


Confiden ce 
Limits 

2 

1 

2 

3 

4 

5 

6 

3 

0 

34 

1-6 

k 

. 25 1 

.12 

.13 

-.10 

-.12 

.02 

+ 0.213 




7-12 

-*11- 

*17 

.13 

• 16 

’ 15 * 

.38* 





13-18 

.27* 

.05 

,05 

-. io 

-.23* 

-.02 





19-24 

-.11 - 

.15 

-.09 

.04 

,03 

.25* 





25-30 

.04 - 

. 10 

,06 

-.12 

-.15 

-.15 





31-36 

-.13 - 

. 19 

-.01 

,01 

-.07 

. 16 


1 

0 

83 

1-6 

-.40*- 

.09 

. 12 

-.12 

10 

.20 

± 0.21R 




7-12 

— . 0 6 — 

. 23 

.18 

.01 

-.16* 

.24* 





13-18 

.06 - 

.13 

.10 

-.02 

-.22 

.23* 





19-24 

-.07 - 

.05 

-.0 6 

,08 

-.17 

.32* 





25-30 

-.06 - 

. 17 

.23* 

-.11 

-.02 

-.02 





31-36 

.05 - 

. 15 

. 10 

.05 

20 

. 20 


0 

1 

72 

1-6 

. 12 

.0 2 

-.00 

,05 

.06 

.03 

O • 23 5 




7-12 

.01 - 

.03 

. 10 

.07 

-.03 

-.35* 

i 




13- 18 

.U 

.15 

-.01 

.02 

-.11 

. 12 





19-24 

.01 

.01 

-.12 

-.06 

.03 

.13 





25-30 

-.17 - 

. 19 

.08 

.0 5 

.07 

-.08 





31-36 

.0 1 - 

.03 

.06 

.0 5 

-.14 

,03 


1 

i 

71 

1-6 

, , * 

-.44 - 

.04 

-.08 

-.00 

.13 

-.07 

+ 0.235 




7-12 

.03 - 

.19 

.28 

- .06 

.13 

-.45* 





13-18 

. 24* 

. 11 

-.10 

.07 

-.20 

.21 





19-24 

- .08 

.14 

-.24* 

.05 

.0 5 

,13 





2 5-30 

-.12 - 

. 16 

.17 

-.02 

.10 

-.15 





31-36 

.08 - 

.07 

.06 

.10 

-.21 

.08 



* 


Significant at 95 % confidence level 
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( 1 - fS ± B ) V 12 w t = ( 1 - 0 1 B 12 ) a t 

where 

w t = X(t) - X (t) with X(t) = Log e Z(t) (3.31) 

Various other series also gave similar results. 

Estimation for ( 1,0/0) x ( 0, 1 , 1 ) Model:- The 
model parameters were estimated by a nonlinear iterative algo- 
rithm ( Box and Jenkins, 19 7o) that minimises the sum of 
squared errors. Values of the parameters were estimated as 
followss 

Sam of square functions s( 0 , 0 ) were computed 

for each time senes over a range of ^ and , values. 

A plot for transit traffic on Lucknow road is shown in 
Fig. 3.11. A minimum is indicated at 0-^ = 0.3 , 0^ = 0,7, 

These values can bo used as starting values for estimation 
of parameters by the iterative procedure. Table 3.6 shows 
that the iterative process converges quickly. The final 
estimates are close to the initial estimates. 

The final estimates of the parameters for other 
series were also obtained by the same procedure and are 
given in Table 3,7, 
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TA3L.E 3,6 ITERATIVE ESTIMATION OF PARAMETERS FOR ARIMA 
MODEL ( TRANSIT TRUCK TRAFFIC - LUCKNOW ROAD ) 


Iteration 

Sum of Squares 

*1 


Starting Value 

1, 40017 

0 . 3000 

0, 7000 

1 

1.35495 

O. 2291 

0,7774 

2 

1,32785 

0.2325 

0.8363 

3 

1.31971 

0.2304 

0,8652 

4 

1,31868 

0 , 2247 

0.8747 

5 

1.31861 

0,2215 

0.8 778 

Final Value 

1.31861 

0. 2215 

0,8778 


Diagnostic Checkings The adequacy of the fit of the 
model was checked by examining the residuals from the 
fitted model. Autocorrelations for any series were not 
significant at 95 percent confidence level. The q statistic 
is given for all the series in Table 3.7, Compared with 
the critical value of at 10 percent significance level 
and 22 Degrees of freedom ( = 3o.8 ) , only one Q value is 
marginally significant though it is also not significant at 
5 percent significance level. The residual series may thus 
be considered to be pure random. Hence the multiplicative 
seasonal ARIMA model of order ( 1 , 0, O ) x ( 0, l, 1) ^2 
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TABLE 3.7 PINAL ESTIMATES OF PARAMETERS FOR MULTIPLICATIVE 



seasons arima 

MODEL ( 

1/ 0, 0 ) X 

( 0, 1, 

12 

Road 

Category of 

-2 

a 

0, 

0 1 

Q 


Traffic 

1 

L 

(22 DOF ) 

Delhi 

Transit Truck 

0,886 

0.4163 

0.7533 

21.95 

Road 

Incoming Truck 

0.38 7 

0.68 22 

0.9117 

25.72 


Bullock Cart 

0 . 670 

0. 1524 

0.7664 

10.8 7 

Allahabad 

Transit Truck 

0.730 

0.3007 

0.7914 

14. U 

Road 

Incoming Truck 

0 * 546 

0.4591 

0.8326 

16.28 


Bullock Cart 

0.43 5 

0.3488 

0.6255 

33.74* 

Lucknov? 

Transit Truck 

0. 508 

0,2215 

0.8778 

27.83 

Road 

Incoming Truck 

0 .406 

0.3837 

0 , 7o6o 

28.54 


Bullock Cart 

0.178 

0.4016 

0.5884 

18.06 

Kalpi 

Transit Truck 

0,7 28 

0.4349 

0 . 79 69 

15.82 

Road 

Incoming Truck 

0,917 

0.5110 

0.908 2 

27,54 


Bullock Cart 

0.453 

0.4253 

0.8193 

23.92 

Kami rpu r 

Incoming Truck 

0.803 

0 . 469 1 

0.8375 

21.35 

Road 

Bullock Cart 

0.956 

0.5134 

0,758 6 

24.66 

2 ( 10 %) 

= 30.8 3d, 2 

( 5 % ) 

= 33,92 




* Significant at 10 % level but not at 5 % level 
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appears to fit the log transformed traffic flow series well. 


3,4,4 Traffic Flow Forecasting 


Stochastic models for monthly traffic volumes can 
be used to forecast future values of the time series. The 
approach to forecasting is described ahd the traffic 
forecasts using seasonal ARIMA model on Delhi road, upto 
April 19 76 are also presented. 

The objective is to make a forecast of Z ( t + L ) 

( or its transform w t ^ ) / at time t for L 1 , Lis 

the lead time and corresponding forecast is known as 

L - step ahead forecast denoted respectively by Z^( L ) or 

w t ( L ) , Box and Jenkins ( 19 7o) show that the minimum mean 

square error forecast w t ( L ) is the conditional expectation 

of w . . at time t , i • e* / 

^ + L 

w t < L ) = E t [ « t ] <3 - 32) 

EUrthcr, the error of estimate of the forecast is 

Error t ( L ) = w t ^ - (l) (3.33) 

The above error function can also be written as (Bari ett, 1955) 


Error t ( L ) = a t+L + ^a t+L ^ 1 + .. 
With L = 1 / Error t (l) = a t+1 = w fc+]> - w t (l) 


“ * ^L-l a t+l 
(3.34) 

that is/a^ 
t 
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series are one step ahead forecast errors. 

For the seasonal ARIMa model with p =1/ d = q = P = 0, 
D = Q = 1/ 

12 2 

( 1 ** 0 ^ B ) ( 1 - B )( ^ B +■ ^ B + ».,,) 

= 1 - 01 B 12 (3.35) 

or 

( 1-0 1 B - B 12 + S^B 13 ) ( * + ^ B + i(i 2 B 2 + ) 

= 1 - OjB 12 (3.36) 

Thus equating the coefficients of B 

*0 = 1 

= 0i 

*2 -01 *1 


*11=01 \ 0 ( 3 ‘ 3 0 

*12” 0i *u + l - ©! 

* 13 = 01 *12 + *1 “ **1 

= 01 *12 

*j = 01 *j-l + *j-12 " 01 *j- 13 for 3 V 13 
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The variance of the forecast error can be written as 

2 2 

V [l] = Var [Error t ( L ) ] =* ( 1+ * + ^ + . . . . 

+ £ , ) a 2 (3.38) 

L- 1 a 

The general ARIMA model can be written m three 

alternative forms; 

(i) as a dif for ence equation; 

(ii) as an infinite su.ro of the currant and previous 

values of shocks a, ? and 

t 

(iii) as an infinite sum of previous observations 

plus the current value of a. . 

t 

Conditional expectation of any of these forms 
gives a forecasting function . In this study, forecasts 
were obtained using the difference equation form. 

The developed model is of the form ; 


U-0 L B ) ( w t+L - W t+L _ 12 ) = ( 1 - S l B 12 ) a t 

(3.39) 

Talcing expectations / 

[ w t+ L ] ~ 01 [ w t+ L- 1 ] + [ w t+ L-12]“ ^I[ w t+L-13 ] 

+ [ a t+ h ] " <9 £ [ a t+L-l2 ] (3.40) 

Thus, for 

L = 1 ! C^U) = 0 1 w t + w t . 12 - 0jat.ll 
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L ;= 2,3,. . . 12 : w^Ch) - 0^ (L- 1) + + L- 12 1 W t+ L - 13 

" 0 1 a t+ L - 12 

Ij = 13 s w t (l3) w t (12) +w S t (l) - fS x 


b = 1 4, 15 , . , jw^.(L) — 0^ w t (L-l) + w^_(L/— 12) — 0 X vJ^(L— 13) 

The standard error of w t (L) is (3.4l) 

22 2 1/7 

t\7 CL> = ( i + * + * + + * T1 ) v * 

1 2 b - 1 a 

(3 . 42) 


The function w t (L) gives the minircum mean square error 
forecasts of X|_(L) - 3^. where x^L) is the logarcthm 

of Z ( t + L ) ♦ ^(L) is a conditional mean which is expected 

to he normally distributed about the actual flow logarithms 
with standard deviation S^(L) . Using the method of moments 
(Mekercher and Dellur, 1972) 

Z^(L) = exp [ X t (L) + (L)/ 2 ] (3,43) 

A r 2 .1/2 

S t (L) = Z t (L) L exp ( Sq (l) ) - 1 J (3.44) 

In this study, traffic forecasts of three different 
categories of goods carriers were made only for Delhi road 
(Fig. 3.12) using seasonal ARIMA model. This model involves 
much less number of parameters compared to the general time 
series model. The forecasts were made for 36 months ahead 
of the last available value (March 1973) of the series. The 
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actual braffic volume data beyond March 19 73 could not 
be obtained due to changed pattern of octroi collection. 

Thus the adequacy of forecasts could not be checked though 
their general characteristics are similar to historical 
data. 

3.5 Analysis of Daily Traffic Flows Within a Month 

Time senes models of monthly traffic flows consider 
the average daily traffic flows of various months. It is 
also desireable to study the daily traffic volumes within 
a month so as to estimate the daily traffic demand. Stochastic 
models have been developed m this study for a few represen- 
tative months. Analysis was carried out for the daily 
traffic data on two roads, namely, Delhi road and Lucknow 
road , during the peak traffic months of December, 19 72, 

January, 19 73 and February, 19 73. 

Figs, 3,13 and 3.14 show the daily volumes of 
three categories of goods carriers during January, 1973, These 
figures indicate that all categories of goods traffic 

i 

fluctuate throughout the month. Fig. 3.13 shows that 
incoming truck traffic volume was minirrum on Sundays in 
both roads. The peak occurs generally on Monday in Lucknow 
road and sometime mid-week in Delhi road • For the transit 
truck traffic ( Fig, 3.13 ) on Lucknow road, the daily 
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variations are smaller wirb the minircum volume occuring 
generally on Mondays. On the Delhi road, transit truck 
traffic also fluctuates throughout the week with no marked 
peak or minimum on a particular day. Bullock cart traffic 
volume shown in Fig, 3,14 indicates minimum traffic flow 
on Sundays similar to incoming truck traffic. The above 
results indicate that there are marked daily variations 
within th<- week for each category of goods traffic. 

Daily traffic volume data were normalised by log 
transformation to X (t) series. Autocorrelations and 
partial autocorrelations were estimated for the following 
three cases: 

(i) undifferenced log transformed series X{t) ; 

(ii) serial differencing of the x(t) senes w.r.t 
days onlyi.e.,d = 1; and 

(lii) weekly differencing of the series ie., D = 1 
and s = 7 days . 

It was observed that for case (iii) , estimated 
autocorrelations are smaller than those of cases (i) and 
(ii) , The log transformed series thus need to be weekly 
differenced before identifying the nature of the model. 

The estimated autocorrelations of all the series 
are shown in Table 3,8. For one series, acf and pacf are 
also shown in Fig, 3.15. The decaying exponential behaviour 
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Road 


Delhi 

Road 


TABLE 3.8 ESTIMATED AUTOCORRELATION FUNCTIONS FOR 
WEEKLY DIFFERED ced data series 


Category 

of 

Traffi c 

Month 

Lags 

Estimated Autocorrelation 
Function s 

12 3 4 

5 

Tran sr t 

12/7 2 

1-5 

» 28 

. 16 

.19 

.15 

.13 

Tm ck 


6-10 

. 11 

.09 

.05 

.06 

.01 


1/73 

1-5 

.31 

.25 

.19 

.13 

.17 



6-10 

. 10 

. 10 

.07 

-.04 

,00 


2/73 

1-5 

.27 

. 12 

.16 

.11 

.07 



6-10 

- .03 

.04 

.01 

-.00 

-.02 

1 n comi ng 

12/72 

1-5 

.35 

. 25 

.21 

.18 

.13 

Truck 


6-10 

.15 

. 10 

- .03 

-.05 

.0 2 


1/73 

1-5 

.32 

.29 

.22 

.24 

.20 



6-10 

.14 

. 13 

.Q6 

.07 

.05 


2/73 

1-5 

.34 

. 25 

.21 

.17 

. 19 



6- 10 

.13 

. 11 

.09 

— .06 

.0 2 

Bullock 

12/72 

1-5 

.45 

. 29 

. 25 

.22 

. 20 

Cart 


6-10 

. 16 

. 11 

.07 

.03 

.01 


1/73 

1-5 

.42 

.30 

.27 

.21 

.19 



6-10 

.13 

. 10 

.10 

.08 

-.04 


2/73 

1-5 

.44 

.34 

.34 

.24 

. 20 



6-10 

.15 

. 11 

.05 

-.07 

.02 


Contd 
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TABLE 3.8 CONTD. ... 


Road 


Lucknow 

Road 


Category Month Lags Estimated Autocorrelation 
of Functions 


Traffic 



1 

2 

3 

4 

5 

Tran sit 

12/72 

1-5 

.12 

.11 

.09 

. 10 

.07 

Truck 


6-10 

.02 

-.04 

.03 

.01 

.00 


1/73 

1-5 

. 10 

.07 

.04 

.02 

.01 



6-10 ■ 

-.01 

.00 

.0 2 

.01 

-.03 


2/73 

1-5 

.11 

.09 

.10 

.06 

.03 



6-10 

.01 

-.03 

.0 2 

-.01 

,00 

Incoming 

12/72 

1-5 

.39 

.31 

.21 

.26 

.20 

Truck 


6-10 

.14 

.09 

.04 

.01 

-.03 


1/73 

1-5 

.40 

.32 

.25 

.21 

.16 



6-10 

.10 

.05 

.0 4 

.05 

,02 


2/73 

1-5 

.36 

.29 

.23 

. 20 

.21 



6-10 

. 16 

.13 

.09 

-.02 

.04 

Bullock 

12/72 

1-5 

.45 

. 25 

.25 

. 18 

. 19 

Cart 


6-10 

.12 

.11 

.09 

1 

* 

o 

.00 


1/73 

1-5 

.48 

.33 

.25 

. 19 

.13 


6-10 

.18 

.19 

. 11 

.06 

.01 


2/73 

1-5 

.48 

.30 

.24 

. 22 

.12 


6-10 

.16 

.14 

.08 

.07 

-.02 
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DELHI ROAD 



a) ESTIMATED AUTOCORRELATIONS 


4 U/y 




33 


0jL acf and substantial pacf at lag one suggest that the 

first order AR model of the form x(t) = 0 X ( t-1) + a*. 

2 

may be suitable. Preliminary estimates of 0 , and a were 

JL 3 

made using Yule-Walker relationship ( Eqs. 3,16 and 3.19) • 

Final estimates were obtained by minimizing the residual 

surn of squares ( Subsec. 3.4.2 ) and are shown in Table 3.9. 

The fined, estimates are close to the preliminary estimates. 

The autocorrelations of the residual seri es were calculated 

and very fev? values are significant cortpared with the 

confidence limits of ± . Values of Q statistic 

fH 

for various series are also shown in Table 3.9, None of 
them is significant at 5 percent level and only one is 
significant at 10 percent level# 

These results indicate that log transformed 
daily traffic series, after weekly differencing, can be 
adequately e>p»ressed by an AR process of first order. For 
peak months 0^ values are very close to each other ( Table 3.9) 
for the same traffic volume series and hence an average 
value of 0^ may be used in analysis. 



TABLE 3.9 FINAL ESTIMATES OF PARAMETERS FOR THE 
FIRST ORDER AR MODEL 


Road 

Category of 

Month 

01 

2 

a 

Q 


Traffic 


X 

a 

( 10 DOF ) 

Delhi 

Transit 

12/7 2 

0. 29 

0.774 

10.9 2 

Road 

Tru ck 

1/73 

0. 28 

0.714 

11.85 



2/73 

0. 29 

0,627 

3.27 


Incoming 

12/72 

0.32 

0.58 5 

5.98 


Truck 

1/73 

0.36 

0.392 

9.86 



2/73 

0.31 

0.547 

8.45 


Bullock 

12/72 

0.42 

0.736 

13.92 


Cart 

1/73 

0.39 

0.952 

3.45 



2/73 

0.46 

1.293 

17.82* 

Lucknow 

Transit 

12/7 2 

0.11 

0.607 

5.95 

Road 

Truck 

1/73 

0 .09 

0.742 

U.67 



2/73 

0,12 

0.738 

12.54 


Incoming 

12/72 

0.37 

0 .506 

12.02 


Truck 

1/73 

0.41 

0.607 

10,05 



2/73 

0.38 

0.325 

9.34 


Bullock 

12/7 2 

0. 43 

0.985 

11.35 


Cart 

1/73 

0.46 

O .886 

16.79 



2/73 

0.47 

0,891 

5.44 

( io * ) 

= 15.99 

x 2 ( 
10 

5 % ) = 

18.31 



* Significant at 10 % level hut not at 5 % level 
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3.6 Analysis of Traffic Flows Within the Day 
3,6,1 Hourly Variations 

traffic volume was recorded at 5 minute intervals 
for 10 days for all categories of incoming and outgoing 
vehicles on G.T, Hoad ( Sec. 2.2) , It v/as found that 5 
minute volume data fluctuate randomly within an hour and 
follow a Poisson distribution for all categories of 
vehicles. Because of daily fluctuations within the week, 
the average daily traffic (ADT) for each category v/as cal- 
culated from data for one week only. The mean hourly 
volumes of different categories in terms of ADT were 
computed from field data for a week. They indicate the 
hourly variations within the day. They are shown for 
both incoming and outgoing directions for different 
categories of vehicles in Figs. 3.16 to3.l8. 

Analysis of daily data for car traffic in a 
week indicated that generally minimum volume occurs on 
Sunday for both incoming and outgoing directions and 
maximum occurs on Monday for incoming and on Saturday for ■ 
outgoing directions. For cars ( Fig. 3,16), outgoing 
traffic is maximum during 1000 to 1100 hours/ whereas 
for incoming traffic, there are twd peaks/ one in the 
morning ( 1000 to 1100 hours) and the other in -the evening 
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( 1700 to 1800 hours) . The peak hourly volume is between 8*0 
to 8.5 percent of ADT. It was also observed that 90 percent 
of ADT travels during 060 O to 2200 hours ( 16 hours ) and 
the volume is insignificant during the remaining eight 
hours . 

Truck traffic ( Fig. 3 . 16 ) fluctuates throughout 
the day with hourly volumes varying generally between 2.2 
to 6.5 percent of ADT, Maximum volume is generally in the 
afternoon between l40o to 1600 hours for outgoing traffic and 
between 130o to 1500 hours for incoming traffic. It was 
also observed that loaded trucks vary between 78 to 87 
percent of the daily truck traffic with a mean of 8 3 percent 
and standard deviation of 4 percent. Thus on an average 
17 percent of the trucks arc errpty. His traffic ( Fig. 3.17) 
is mostly during the day with 9o percent of ADT between 0600 
to 2100 hours. Incoming bus traffic is maximum during 0900 to 
11O0 hours with maximum hourly volume being 9 percent of 
ADT. peaks for the outgoing traffic are around 7 percent 
of aDT and they occur between 0800 to 1000 hours and 1700 to 
1900 hours. 

The daily variations of scooters and motorcycles 
are similar to those of cars , Hourly variations, shown in 
Fig. 3,17, are quite significant. Maximum hourly volume is 
as high as 14 percent of aDT and the peaks occur between 
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1000 to 1100 hours for outgoing traffic and 1700 to 1800 
hours for incoming traffic. Also 95 percent of daily traffic 
moves between 0700 to 2200 hours, Bicycle traffic is minimum 
on Sunday. Hourly variations ( Fig. 3.18) indicate that 
incoming traffic is maxirrum botv/cen 0700 to 1000 hours and 
outgoing traffic between 1700 to 1900 hours, peak hourly 
volume is about 12 percent of ADT and more than 90 percent 
of daily traffic travels between 060 O to 2100 hours. 

Hourly variations of animal driven vehicles 
( tongas and bullock carts), shown in Fig. 3. 18 , indicate 
that the maxirrum outgoing traffic moves between 1500 to 1800 
hours and is ll,5 percent of ADT; whereas peak for incoming 
traffic is between lOOo to 1200 hours and is only 10 percent 
of adt. 

3.6.2 Peak Hourly Conposition of Mixed Traffic 

Hourly variations of different categories of traffic 
indicate that the peaks occur at different times for the 
different categories of incoming and outgoing traffic. For 
the outgoing traffic maximum volume for ca^s and scooters 
is generally in the morning between lOoo to 1200 hours 
with a smaller peak in tho evening. For other categories 
of vehicles , the peak is generally in the evening. Thus 
the critical period for outgoing traffic seems to be between 
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16oo to 1800 hours in the evening . 

For the incoming traffic/ peaks for buses/ trucks/ 
bicycles and slow moving vehicles occur generally in the 
morning; whereas the maximum volume for cans and scooters 
occur in the evening with a smaller peak in the morning. 

Hence the critical period for incoming traffic seems to be 
in the morning during 1000 to 1200 hours. The peak hour volumes 
of different categories of vehicles and the corresponding 
compositions for the tv/o critical periods are given in 
Table 3.io. 

table 3.10 peak hojr volumes and compositions 


Outcoming Incoming 


Category 
of Vehi- 
cle 

ADT 

Morning 
( 1000 - 
1200 ) 

%of APT 

Evening * 

( 1600 - 
1800 ) 

°/ 0 of ADT 

ADT Morning* 

( looo- 
1200 ) 

%o£ ADT 

Evening 
( I 600 -I 8 OO) 
% of ADT 

Cars 

484 

8.5 

6.6 

431 

7.8 

8.1 

Tracks 

377 

4.5 

5.5 

380 

6.0 

4.0 

Bu sos 

128 

7.0 

6.8 

121 

9 .O 

6. 4 

Scooters 

352 

14.0 

G.O 

3 43 

7.0 

12.5 

Bicycl os 

2341 

5.5 

12.0 

2418 

7.0 

7.0 

Tongas 

63 



6f| 





7. 0 

w 

• 

LT 


10 .0 

7.7 

Bull ock 

47 



51 i 



Carts 








* l^parantly critical corvposition 
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Subsequent studies ( Chapter 6 ) indicate? that 
the PCUS for mixed traffic for both directions are 
nearly equal during the peak hours. 

3,6,3 Estimation of peak Hourly Volume 

Field data and their analyses indicate that the 
traffic volume varies from month to month within and over 
the years; from day to day within a month; and from hour to 
hour within the day. Stochastic models for monthly and 
daily goods traffic on five roads approaching Kanpur and 
hourly variations within the day for one road, have been 
determined. The critical peak hours within the day and 
compositions of traffic for the critical periods have also 
been determined, it seems possible to determine the DHV 
by considering the different variations, e.g., from the 

monthly and daily variation of traffic, it may be possible 
to estimate the ADT for different categories of vehicles 
with a given recurrence interval. From the composition 
of vehicles at critical periods, the volume of each category 
of vehicle in the critical period can be determined. 

When PC E s for different categories of vehicles for given 
volume and composition are known, the equivalent pcUs for 
mixed traffic in critical periods can be determined. The 
highways can then be appropriately designed. 
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4. MATHEMATICAL MODELLING 

4,1 introduction 

A dynamic traffic system is characterised toy the 
occurrence of changes in the system variables. A mathematical 
model consists in representing analytically the system 
relationships in terms of the parameters of the system and 
hence representing the changes in the system in analytic 
terms. The model should represent the system as closely as 
possible so that the estimates are meaningful, ?,.ccuracy of 
prediction can generally be improved only by making the 
model more corrplicated and hence less convenient to use, 
particularly when a sinple model may not be satisfactory. 

There is hence a necessity for a compromise between accuracy 
and simplicity ( Hall, 19 62; Chestnut, 1967) . 

The advent of digital corrputers and developments 
in system analysis and optimisation have radically changed the 
approach to mathematical modelling of physical systems. 

While lack of solution techniques was a prime hinderance 
in the use of sophisticated and more realistic system models, 
at present the level of sinplicity is governed very often 
by lack of knowledge or understanding of the system. 
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4.2 Mathematical Models for Traffic Flow 

4.2.1 General 

The objective of traffic flow modelling is to derive 
theoretical or empirical relationships between the variables 
so as to determine the characteristics of the traffic 
stream. There are three variables of particular interest/ 
viz., flow , density and speed; and they together define 
the quality of service experienced by the drivers ( Wohl 
and Martin, 1967; Haight, 1963) , A complete behavioural 
understanding of the traffic stream phenomenon is currently 
not available. Most of the earlier studies were designed to 
develop workable approximate relationships. In recent 
years, some of the models based on operations research 
techniques have also been developed. Some of the deterministic 
and probabilistic models that describe the stream chara- 
cteristics of highways are briefly described. 

4.2.2 Empirical Methods 

The deterministic approach began with a suitable 

algebraic equation to explain the flow concentration 

relationships. Green shields (193 4) gave a relationship 

of the form; __ 

Uf 

tr B = tr f - ( ) d ( 4. i) 

D 

j 
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whore U s = average space mean speed; U £ = mean free speed; 

D = density; and Dj = jam density. 

Greenberg (1958) assumed that high density traffic 
flow is analogous to continuous fluid flow and based on 
field studies fitted an e:xponential function to speed density 
observations, viz., 

D = C e b1T s (4.2) 

where C and b are constants, 

A number of empirical relationships based on travel 
time and volume data have been suggested by Guerin (196l) , 
Normann (1942) , etc. They are of the following general 
type. 

V = C + b U s (4.3) 

where Vis the volume; and b and C are parameters. 

4.2.3 Deterministic Models 

Hydrodynamic Analogies? Analogies have been drawn between 
the flow of fluids and the movement of vehicular traffic. 
Principal contributions have been made by Greenberg (1958 ) , 
Lighthill and Whitham ( 19 64 ) and Richards (19 56) , The basic 
assumption is that high density traffic behaves like a 
continuous fluid. The fundamental equation is; 

d U s c 2 3D 

9x 


d t ” 


D 


(4.4) 
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where U s = fluid velocity or space mean speed in krrph; D = 
density in vehicles per km; x = distance in kms; t = time to 


travel distance x and c = roadway parameter. 

The solution of Eq. 4.3 yields the speed xn 
terms of density as follows: 


U 


s 


c Log e 



( 4.5 ) 


Lighthill and Whitham ( 19 64 ) have shown that 


the speed of waves carrying continuous changes of volume 


through a vehicular flow is given by 

d V dl D d U , 


u w = 


a d 


d D 


U s + D 


a d 


(4.6) 


The interaction between the vehicles is very 
small at low densities and wave velocity equals the vehicle 
speed. At densities above the point of maximum volume, the 
waves will move backward relative to the road and at maximum 
volume, the wave is stationary relative to the road. Fig. 4.1 
shows the speed volume density relationships. Due to the 
changes in wave velocity with vehicular density, it is 
possible to have different waves traveling through a 
traffic stream, when traffic moves from a low density 
road section to a high density section, the waves in the 
low density traffic will travel forward at a higher speed 
than the waves in the high density flow. The new wave 
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formed on the meeting of the two waves is termed as a 
shock wave ( Fig, 4,1), some of the examples dealing 
with different conditions occuring in vehicular flov; have 
also been analysed by hydrodynamic analogies ( Lighthill and 
Vjhitham, 19 64) , 

Car Following Models; A driver probably reacts depending 
on the distance he is away from the vehicle in front and 
its rate of acceleration or deceleration, A car following 
model relates the movement of a single vehicle to the 
vehicle it follows. Herman and his associates ( 1959, I960 ) 
at the General Motors Research Laboratories have studied 
this problem extensively. They assumed that a driver will 
try to keep the relative speed between his vehicle and that 
ahead, as small as possible. The mathematical representation 
is as follows: 

a n (t) =K[U n _ 1 ( t - T s ) - U n ( t- T s )] (4,7) 

whore a^t) = acceleration of nth vehicle at time t; 

U n (t) = velocity of nth vehicle at time t; T s = time lag 
or stimulus - response time of driver - car system; and 
K = driver - sensitivity coefficient. 

The driver's sensitivity to the other vehicles 
increases as the vehicles get closer to one another , 




Ej<pcriments have also shown that K varies with spacing 
between the vehicles. Herman (I9 60 ) modified Eq. 4.7 as 
fol low ss 
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a^t) = [ u n _ x ( t - T s > _ u n < t- Tg)] 

mi-1 _x n 

I (4.8) 

where K 0 = a constant; 3 ^^ and are positions of ( n - 1) th 
and n th vehicles. 

The general solution of the above equation indicates 
that the spacing between successive vehicles has oscillatory 
characteristics that depend upon the nature of the constant 
C , defined as 

C » K c T s ( 4.9 ) 

Edie ( i 960 ) has suggested that driver - sensitivity 
coefficient ( K ) be taken as a function of both vehicle 
speed and vehicle spacing, i.e., 


= 


K 2 V * - ) 

< *n-l “ *h * 


f~ [ Vl ( - T s>- U n <t-T ) 3 


S' 

(4. IO) 


where K 2 is a constant. It leads to the following 
rel ation ship . 


U. 

D = D . Log 

J e U 


s 


where Dj, the jam density = — 

J K t 


( 4.11) 
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4 , 2. 4 p robabil i sti c Model s 


Due to the probabilistic nature of headway gaps, 
free speeds etc. ( Sec, 2.3 ) , traffic flow is probabilistic. 
Based on queuing theory, flow travel time relationships 
have been developed by Davidson ( 19 66 ) . He defines the 
service time as the travel time of a vehicle when no other 
vehicle is using the road. Queuing delay is taken as the 
increase in travel time due to the presence of other vehicles 
on the road. The delay time for a vehicle in queue for 

random arrival and random service is given by ( Cleveland, 19 64) . 

2 


a = 


(4.12) 


u(u-a) a ( 1 - c ) 
where d = the delay per vehicle; a = arrival rate; u = service 
rate; and c = a/u . 

The terms of Eq, 4,12 may be translated into traffic 
stream variables as follows? a = V = flow, u = Sjj = the 
saturation flow and c = v/Sf. The ratio of the delay(d) in 
the queue to service time is given by 


d 

1/u 


2 

c u 


( 4.13) 


a(l-c) 1-c 
Let t Q = 1/u be the travel time at no flow. Then 
from Eq, 4.13 , the delay time is given by 

c 


d = h 


1-c 


( 4.14 ) 
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Traffic flow is not truly a single continuous 
queuing situation. Delay is caused by a succession of 
queuing situations such that a varying amount of the total 
service time is subject to queuing delays. Davidson ( 1966) 
modified Eq, 4.14 to 

d = t J —1— ( 4. 15 ) 

° 1 - c 

where J is a factor that varies with road type and frequency 
of delay producing situations along its length. The total 
travel time t can be expressed as 

t = t Q ( 1 + J ) (4.16) 

° 1 - c 

4,2.5 Simulation Models 

When suitable conceptual models cannot be formulated 
for the several components of the process and the complexity 
is due to the stochastic nature of the process and coherent 
system components and perheps the presence of multiple 
alternatives and constraints/ computer simulation of the 
process may be adopted. Actually simulation analysis does 
not solve any mathematical problem. It consists in setting 
up of a mathematical model of the system, scanning the 
process through the system and study the characteristics 
of the process as it evolves in time ( Mize and Cox, 19 68) . 
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In simulation, it is possible to vary the parameters of 
the system, to include stochastic inputs or system responses, 
to investigate the independence or interdep on donee of the 
elements and also, if necessary, to adept sanpling techniques 
to determine an optimal or near optimal design of the system 
( Wohl and Martin, 19 67 ) , 

Analog computers are preferred when the model 
consists largely of differential equations; if conpontents 
are to be tested separately? frequency of repetition is 
large? and for 1 real-time 1 and for 'man-machine' system 
simulation. Digital corrputers are preferable if the model 
consists of a large number of arithmatic or logical operations 
and difference equations and for large size problems 
{ Drew , 19 68) . 

Simulation is hence very suitable for analysis 
of conplex stochastic processes with interacting components. 

It has been used to a limited extent to study freeway merging 
operations, interchange design, signalised intersections 
etc. ( Gerlough, 19 64) , and are generally limited to the 
consideration of homogeneous traffic. 

4,3 Mathematical Modelling of Mixed Traffic Flow 

4.3.1 Corrplexity of the problem 


Mathematical models described in Sec, 4.2 describe 
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the flow characteristics of motor vehicles only. The 
problem of analysing mixed traffic flow is quite complex 
because of the following factors: 

(1) The traffic flow is stochastic in nature and 
the free speeds of vehicles are also probabilistic, 

(ii) There are large variations in the speeds of 
the fast and slow moving vehicles; thereby increasing 
the interaction, 

(iii) Slow moving vehicles like bullock carts are 
unaffected by other stream characteristics like volume 
level and composition and they continue to move at their 
respective free speeds. 

’ (iv) Vehicles like trucks, cars and bullock carts 
etc, occupy one full lane, whereas two wheelers like 
bicycles, scooters etc, occupy very little width and two 
of such vehicles can very easily move parallel in one lane. 
Thus the flow logics for two wheelers are significantly 
different from big vehicles. 

(v) The overtaking operations of the fast 
moving vehicles are restrained in single and double lane 
roads* due to nonavailability of sufficient gap in the 
opposite traffic stream, 

- (vi) The slow moving vehicles have a tendency to 


move in platoons 


A fast moving vehicle has thus to 
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overtake a full platoon in one operation, thereby needing 
a lot of time and free gap . 

4.3.2 Limications of Analytic Techniques 

Analysis of the mixed vehicular traffic flow by 
analytic techniques is difficult due to their limitations 
wh i ch i n clu de : 

(i) Mathematical equations developed from 
hydrodynamic analogies hold good only for high traffic 
densities as they rely on the assumption that high density 
streams of vehicles behave like an incompressible fluid. 

The character of intervehicle relationship has not been 
taken into account in the analogy. 

(ii) Car following models are generally restricted 
to the flow of homogeneous vehicles in a single lane with no 
overtaking. In mixed traffic flow, overtaking is quite important. 

(iii) Flow-travel time relationships take into account 
random distribution of traffic flow. But they are limited to 
the consideration of homogeneous traffic with constant free 
speed. 

4.3.3 Need for simulation 


Mixed vehicular traffic flow is very complex and 
it does not seem feasible to use analytic techniques for 



studying the nrocess. Yet it is possible to simulate the 
process on a computer. Because of a large number ofj 
logical decisions involved in the flow process, digital 
simulation is preferred. Corrputer simulation has tho 
following additional advantages. 

(i) It may lead to a better understanding of the 
nature of the interactions involved including acceleration, 
retardation, overtaking etc, ( Drew, 19 68), 

(ii) It is useful to validate the component 

model s. 

(iii) Data can be generated and used to study the 
transient and steady state characteristics of the system 
( Naylor et al * , 19 66; Mize and Cox; 19 68). 

Hence it is proposed to use digital computer 


sirrulation for the analysis of mixed vehicular traffic 



5. 


COMPUTER SIMULATION 


5.1 G en er al 

Digital computer simulation may be defined, as, 

"a numerical technique for conducting experiments with 
certain types of mathematical models which describe the 
behaviour of a corrplex system on a digital computer over 
extended periods of time" (Naylor et al., 1966) . simulation 
serves different functions in system analysis and design. 
They include aids in problem formulation; insight into 
the sensitivity of design to wide ranges of parameters; 
guidance in predicting system performance, and validation 
of system model (Reitmann, 1971) . 

5,1.1 Steps in simulation 

There are generally a number of steps in simulation 
of an engineering system. Though the details may vary 
from time to time, a number of steps are common to several 
simulation studies. Fig. 5.1 is the representation of 
the various steps for simulation (Naylor et al. 1966 ; Mize 
and Cox, 19 68) . 

(I) Formulation of the problem: The problem is 


formulated in analytic terms. It includes the identification 




FIG, 5.1 FLOW CHART FOR SIMULATION 
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of the components of the system and its environment. 
Constraints on tho system are also to be identified, 

(II) Collection and analysis of data; Historical 
and experimental data concerning socioeconomic, engineering 
and other factors are collected and analysed. 

(Ill) Formulation of the system model; Based on 
experience, available Ten owl edge and information, a system 
model is formulated to represent the corrponents of the 
system, their parameters and behaviour and their relation- 
ships; the nature of the inputs to the system, throughputs, 
and outputs. 

(IV) Formulation of a computer programme: The 
formulation of a computer programme for the purpose of 
conducting simulation ejqperiments with a model of a complex 
system requires special considerations to the following 
three activities; (i) computer programme; (ii) data input 
and starting conditions; and (iii) data generation. 

(V) Estimation of parameters: Parameters of the 
system are estimated using historical and/or experimental 
data, 

(VI) Validation; The assumed mathematical model 
for tho system is a simple approximation to a more corrplex 
reality. It is necessary to validate the assumed system 
model in order that the mathematical model is a reasonable 
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approximation to reality* Validation involves the definition 
of reality and a reasonable approximation to it. Generally 
validation is done on the basis of corrparison between the 
actual recorded output and the output from simulation 
model for corresponding inputs. In cas^ the model is 
considered to be unsatisfactory, it is necessary to modify 
the model suitably in order that estimated output from the 
model agrees with the observed record, 

(VII) Design of sinulation experiments: Simulation 
may have to be done for analysis and optimisation of 
complex systems under different conditions of system inputs 
and outputs, and values of design variables, if any. The 
values of these variables are generally to be sampled. 

Hence it may be necessary to select factor levels and 
combination of factors and the order of experimentation . 

It is also necessary to ensure that sampling errors arc 
within bounds. 

(VIII) system simulation: The system is simulated 
under different conditions as per the experimental design. 
When random processes are involved, it may be necessary to 
generate and use data concerning such variables. 

(IX) Analysis of simulated data: Results from 
simulation are analysed to determine the characteristics 
of various corrponents of the process, their significance 
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and interdependence* a near optimal design may be obtained 
in appropriate cases* 

5.2 Model Formulation 

5.2.1 System Description 

The details of the system involving various components, 
variables, relationships and assumptions involved in the 
mathematical modelling of system simulation are discussed 
briefly. 


System Components; A graphical representation of the 
roadway section including the various system components 
is shown in Fig, 5.2, The components includes 



h 


-LENGTH = 2250 m 



PIG, 5.2 ROADWAY SECTION AND VARECUS COMPONENTS 
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(i) Length of roadway section - 2,5 km. 

(ii) VJidth - 2 lane , 6.0 metre wide cement concrete 
pavement with 1.5 metre wide shoulders on each side, 

(iii) Lateral obst motions - no natural or man made 
obstructions are present to restrict free flow. 

(iv) Horizontal alignment - the section is a 
straight reach with no sight distance restrictions duo to 
bends etc. 

(v) Vortical alignment - the longitudinal gradient 
is quite flat and does not affect the free flow. 

(vi) Presence of intersections and other regulatory 
and control measures. The section under study vras free 
from all such restrictions. 

Assurrptions: The following general assumptions have been made 

in model formulations* 

( 1 ) The roadway section was considered to be 
sufficiently away from traffic control points, e.g . , 
intersections at both ends. Vehicles were assumed to 
approach the section in a random fashion. The section 
selected for field studies generally satisfied this require- 
ment, 

(ii) Vehicles approach the section from both 
directions at their free speeds. The free speed of a certain 
category of vehicle is probabilistic and perhaps also a 
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function of time and space. However, the free speed distri- 
bution parameters we re assumed to be constants v/ith no 
variations in time and space. 

System Variables: Variables are those attributes of the 
system that take on different values under different 
conditions, or in different system states. They can be 
classified into several ways as exogeneous, endogeneous, 
status and random variables, etc. An exogeneous ( or input ) 
variable ( Naylor et al., 19 71; Mize and Cox , 19 68 ) is 

one that is quantifiable and which affects, but is itself 
unaffected by the system. Status variable is one that 
represents the state of an important characteristic of the 
system, which is generally a throughput variable. A variable 
having a value determined by other variables in the system, 
is an endogeneous variable ( Naylor et al,, 1971) and 
represents in a way the output of the system. The following 
variables may be identified in the traffic system under 
study; 

Exogeneous Variables: They include: 

(i) Interarrival time gaps A3 (I) and AS3l( J) of the 
vehicles approaching from either direction. The vehicles 
are to be identified and processed simultaneously from two 
directions. So separate names we re given to variables for 
either direction with the letter 1 before brackets identifying 
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the characteristics of vehicles in lane 2. Subscript I 
pertains to the serial number of vehicle in lane l and J to 
that in lane 2. Interarrival time gaps are used to compute 
the arrival times AT(I) and AT1(J); 

(ii) Hourly volume of traffic NX and NXl and the 
number of categories of vehicles in the mix NCAT, NCATi 
alongwith their proportions PROB(K) and PROBI(k) for 
K = 1 , NCAT (NCATI) ; 

(iii) Category of the vehicle entering the section 
from either direction ICATtI) and ICATi(J) . ICAT and ICATi 
are defined by codes 1-6 as follows: 

ICAT(i) or ICATl(J) = 1 for cars; 2 for trucks; 

3 for bullock carts; 4 for tongas; 5 for scooters; and 
6 for bicycles; and 

(iv) Free speed of arriving vehicles, V3?(l) and 
VF1(J). Free speeds follow a censored normal distribution 
for each category of vehicle* 

Status Variables: They include: 

(i) position of vehicles P(I) andpi(J) along with 
their running speeds V(I) and V1(J) at any instant; 

(ii) Available headways D(I) and D1(J) for each 
vehicle at any instant; 

(iii) Overtaking time OT(l) and OTl(J)in case a 
vehicle is overtaking at any instant. Otherwise OT(i) and 



114 


OTl(j) arc equal to zero; 

(iv) Traffic stream logic for the vehicle at any 
instant i.G,,, NVT(I) andNVTl(J) . Different stream logics 
arc indicated by codes; and 

(v) Number of vehicles in the section i.e., density 
of each category IQE(K) and IQE1(K) for K = 1, NCAT (NCaTI) 
at any instant. 

Endogeneous Variables: They include: 

(i) Actual departure time DEPT(I) and DEPTl(J) of 
the vehicles from the system. These are used to compute 
delay times for each vehicle DELT(i) and DELTKj) ; 

(ii) Number of vehicles of each category going out 
of the section NDPT(k) and NDPTi(K) 7 

(iii) Operating speed of the departed vehicles VEF(i) 
and VEFi(j) 7 

(iv) Number of vehicles delayed NDEL (K) and 
NDELl(K) in each category; and 

(v) Average operating speed AWE(k) and AWEi(k) for 
each category. 

Random Variables: The following three variables are 
probabilistic in nature ( Sec* 2,3 ) : 

(i) Arrivals afe random and interarrival time gaps 
AP(T) and AGl(J) follow a shifted exponential distribution; 
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(ii) Category of the arriving vehicles ICAT(I) and 
I CAT 1 ( J) . The probability that arriving vehicle belongs 
to a specific category is considered to be proportional to 
the proportion of the vehicles of the same category in the 
overall traffic composition; and 

(iii) Free speed of vehicles VF(I) and VF1(J). The 
free speeds of a particular category follow a censored 
normal distribution. 

System Relationships: Some of the general relationships 
between the variables of the mathematical model ate given 
below. Relationships that are more specific to the sirrulation 
analysis are given subsequently in Sec. 5,4. Equations are 
given for the two lanes separately only where necessary. 
Otherwise relationship for lane 1 is given and a similar 
relationship holds good for lane 2, 

(i) Available headway for a vehicle at any instant, 
i.e., D(l) and D1(J) . The vehicles from either direction 
are arranged in the order of their positions L and LI 
and KK(L) = I and KKl(Ll) = j. The available headways can 
be calculated as follows: 

D(I) = P ( KK ( L - 1 ) ) - P(I) (5.1) 

and D1(J) s pl(j) - PI ( KK1 ( LI - 1 )) (5.2) 

(iii) Eixpected departure time of the vehicles* i.e., 
EDEPT(i) and EDEPTl(j) , if they continue to move unimpeded 
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at their free speeds throughout the length of the section t 
i.e., there is absolutely no delay. 

LNGT 


EDEPT(I) = AT (I) + 


IX 


(5.3) 


where at(I) = arrival time of i th vehicle in lane 1? 

LNGT i= length of the section; and 

IX = distance covered by Ith vehicle in one second 
at free speed VF(I) in kmph 

= 0.28 VF(I) 

(iii) Delay time of a vehicle DELT(I) and DELTl(j). 
DELT(I) = DEPT(I) - EDEPT(I) (5.4) 

(iv) Operating speed of the vehicle in kmph 

LN3T 


VEF(I) = 


(5.5) 


0.28 ( DEPT (I) _ AT(D) 

(v) Average opeiating speed of each category of 

vehicle AVVE(k) and AVVEi(k) and their standard dfei-i-v-a-ti-on s / 


SDVEF(K) and SDVEFI(k) . 

1 


AWE(K) = 


NDPT(K) 


ndpt(k) 


1=1 


VEF(I) 


(5.6) 


SDVEF(K) 


E ( VEF( I) ) 2 - NEPT(K) (AWE(K) ) 2 


NDPT(K) - 1 


(5.7) 


5,2.2 Scanning Techniques 


The process is simulated along time and scanned or 
observed at intervals. The computer programme mast be 



117 


written in such a way that it moves the model over simulated 
time, causing events to occur m proper order and with a 
proper time interval between successive occurrences of 
events. 

The real components function simultaneously whereas the 
components are simulated sequentially. There are dependencies 
among different parts of the real system. It is important 
that a time flow mechanism be embedded and the process 
scanned so that simulated performances of the system compo- 
nents are synchronised in time^ This can be accomplished 
in two ways ( Blake and Gordon, 1964) ; (i) periodic scanning 

or uniform increment method; and (ii) event oriented scanning 
or variable increment method. 

In periodic scan method, the duration of the 
simulation run is divided into a number of equal time 
intervals, that is, scan and update the entire system ( like 
move vehicles, change speeds etc.) once per unit time 
interval. In the event scan method, when a specific event 
has ocoured, the status of a system is updated and among 
the possible events that can occur, a scanning is done to 
determine when the next event occurs and the system is 
again updated. The procedure is repeated. Hence event 
scanning essentially deals with the scanning of the process 
at nonuniform time intervals depending upon the occurrence 
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of the next ^vent; while periodic scanning is done at 
uniform time intervals. Event scan technique is faster and 
more accurate due to less variance of estimates/ hut it 
usual ly retires greater programming complexity . periodic 
scanning is more straight forward and easy to programme. 

The traffic flow process requires scanning both 
along time and space in order to determine the position and 
other characteristics of vehicles at different times. 
Continuous scanning in time and space is not possible in 
digital sirailation and so time and spacing were discretised. 
The roadway space is divided into one metre long one lane 
wide blocks. The vehicles are moved into the system based 
on their own characteristics and aJ so their interactions 
with other vehicles in the system. At different timos, a 
vehicle may be accelerating/ running at free speed, overtaking 
or decelerating depending on traffic conditions. With so 
many vehicles in the system having varying stream logics 
(Sec. 5.3 ) at different times, it is very difficult to 
identify the immediate next event. In such a case where 
vehicular characteristics like speed etc., are time dependent, 
periodic scanning is easy to formal ate and was adopted in 
this study. scanning was done uniformly at one second 
intervals based on the vehicular characteristics at that 
instant, and vehicles were stepped up one metre long 
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one lane wide blocks. 

New vehicles are brought into the system with 
arrival gaps following a certain probability distribution. 
Arrival times of all the vehicles likely to enter during 
next specified sirrulation period from either directions are 
generated and stored, Sirrulation starts with periodic 
scanning and at each time, it is tested whether a vehicle 
is likely to enter from any direction or not. If a vehicle 
enters, it is subsequently moved and scanned periodically. 

For overtaking vehicles, the overtaking time and the position 
of vehicle on corrpleting the operation are cortputed. Hence 
for the overtaking vehicle, scanning is not needed for this 
duration. But on completion of overtaking operation, the 
vehicle is brought back into the system for periodic scanning. 

5,2,3 Computer Representation of the Simulation Model 

Digital simulation deals with studying the process 
as it takes place through the system. This is acconplishod 
by representing the corrponents or their characteristics and 
the variables in the coirputor and manipulating them as the 
process takes place along time. Mixed traffic flow on the 
two lane highway involves the representation and manipulation 
of the vehicles over the two lanes. There are two general 
procedures for such representation ( Gerlough , 19 56, 58; 
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Wohl and Martin, 19 67) t namely, (i) Memorandum representation 
and (ii) Physical representation. 

In memorandum method, the characteristics of the 
vehicles like category, speed and position etc. are indicated 
by using one or more coded words, generally one word per 
characteristic. In each time interval, the data are updated 
and the imformation concerning the characteristics is changed, 
if and when necessary. 

In physical representation, the memory locations 
are associated with blocks in physical space. The information 
in the memory location characterises the information concer- 
ning the presence or absence of the vehicle in the given 
block and their characteristics. The data in the memory 
locations are scanned to identify the characteristics of 
the traffic and are updated in every time interval. 

There are 4500 one metre long one lane wide blocks 
to represent the space of the roadway section. Further more, 
there are large number of characteristics of an individual 
vehicle with wide variations in speeds and other character- 
istics. Because of all these factors, memory requirement is 
very large for physical representation. In addition, data 
manipulation is also very complicated. Hence in this study. 


it was decided to adopt memorandum representation of the 
simulation process. The vehicles were identified individually; 
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their characteristics were stored separately; and were 
updated and analysed in each time interval. 

5.3 Traffic stream Logic 

5.3.1 Introduction 

The logic that the vehicles follow as they move 
through the system is referred to as traffic stream logic 
(Drew, 19 68) . The logic for any vehicle depends upon its 
position and speed relative to other vehicles moving in 
both directions and as the process is dynamic, it will 
change from time to time. The stream logic was classified into 
the following three categories for this study: 

(i) Plow logic for unimpeded vehicles, i. e», vehicles 
that continue to move at their free speeds; 

(ii) Flow logic for overtaking vehicles that decide 
to overtake; and 

(iii) Flow logic for restrained vehicles i.e«for 
vehicles that desire to overtake; but are restrained due ’ 
to opposite traffic stream. 

5.3.2 Spacings for Vehicles 

A vehicle has to keep a certain minimum spacing 
with the vehicle ahead of it so as to avoid collision if 
the forward vehicle suddenly decelerates or stops. Spacing 
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is defined in this study as the distance measured from 
centre to centre of successive vehicles, The minirrum 
spacing is generally a function of driver reaction time, 
speed and length of vehicle, coefficient o£ friction, etc, 

A certain minimum spacing is also needed when a vehicle 
overtakes another vehicle and reenters its own lane. The 
spacing will govern the position of the vehicle after 
overtaking and hence the overtaking time. 

When the vehicle ahead is moving faster than the 
rear vehicle, headway availsbl e to the rear vehicle 
continues to increase. However, if the rear vehicle is 
moving faster, then it will either try to keep the minimum 
spacing or else overtake. The decision to overtake is 
assumed to be taken when the available headway is equal to 
or less than the specified rnaxiirum spacing for overtaking. 

Based on empirical field observations and available specifi- 
cations, these spacings are to be specified. 

Minirrum spacing for Vehicles; The minirrum spacing for vehicles 
is a function of speed and length of vehicles. The 
following equation ( IRC, 19 50 ) was adopted in this study, 

S(I) =0.2 V(I) + h 0 (5.8) 

where s(l) = minirrum spacing; V(I) = speed of vehicle I; and 
L 0 = constant which varies with the category of vehicles 
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and may be a function of length of vehicle I and of that 
ahead of it. 

Spacing for Overtaking Vehicles: A vehicle desiring to 
overtake needs more spacing to permit the driver to mane uvre 
the vehicle to che other lane. Maximum spacing for over- 
taking is a function of the speed of vehicle and also the 
drivers reaction time for taking the decision. This 
spacing is thus defined as the distance that can be covered 
by th'- 1 vehicle during the headway time of T q seconds/ viz., 

SD(I) = v(i). To 

= 0.28 V(l) .T 0 { 5.9 ) 

where V(I) is the speed of the overtaking vehicle in krnph 
and T 0 is the headway time in seconds which is a parameter. 

For slow moving vehicles, the maximum spacing for 
overtaking SD(I) , calculated on the basis of above equation 
may be smaller than the minimum spacing S(I) , in which 
case SD(I) was limited to S(I) . This was the case when a 
tonga or bicycle overtakes a bullock cart. 

5.3.3 Flow Logic for Unimpeded Vehicles 

Let D(I) be the available spacing between the 
vehicle I and that ahead of it at any instant ( Fig. 5.2 aru ^ 
Eqs, 5.1 and 5,2. Tho available spacings at an instant can be 
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computed by arranging all the vehicles moving in a lane 
in the order of their positions. When the free speed of a 
vehicle is less than that of the vehicle ahead, the rear 
vehicle can continue to move unimpeded irrespective of the 
available spacing which is increasing with time. When 
the available spacing D(l) for a vehicle at an instant 
is more than the maximum spacing needed for overtaking i.e,, 
D(I) > SD( I) 7 then rolso this vehicle can continue to move 
unimpeded at its free speed in its own lane. Otherwise the 
speed of the vehicle will be affected by traffic conditions. 

5,3,4 Overtaking Logic 

It may be noted that L gives the ordered position 
of the vehicle in a lane, KK(L) gives the sequence number 
of the vehicle based on arrival times, i» e, ,KI<(L) = I and 
KK1 ( Ll ) = J. 

When the speed of a vehicle is more than that of 
the vehicle ahead, i.e., V(KK ( L ) ) > V (KK ( L - 1) ) and 
the available spacing D (KK(L) ) is less than or equal to 
the maximum spacing for overtaking, ie., D{KK(l) SD(KK (L) ) ; 
then the vehicle has a tendency to overtake r 

(i) When the speed difference between two successive 
vehicles is less than 16 krrph, i.e,, 

V(KK(L-1) ) < V (KI<(L) ) < V (KK(L-l)) + 16 , 
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the roar vehicle can overtake and the following assumptions 
concerning the overtaking operations are made, 

Hie rear vehicle slows down to the speed of the 
vehicle ahead and follows it in its own lane for a time 
Tr , when the driver looks in the opposite stream for a gap 
sufficient enough for overtaking. If overtaking is possible, 
the vehicle accelerates, goes to the wrong lane and moves 
back to its own lane after overtaking. The rate of 
acceleration varies with speed and Indian Roads Congress 
(IRC, i9 5o) specifications for motor vehicles given in 
Table 5,1 are used. Fig, 5. 3 shows the details of the 
overtaking operation. 

table 5.1 maximjm overtaking accelerations at different speeds 


Speed 

knph 

Maximum Overtaking 

Accel oration 

Icnph/ sec 

25 

5.00 

3o 

4.80 

40 

4.45 

5o 

4.00 

65 

3.28 

80 

2.56 

100 

1.92 
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Let and represent the positions of two 
successive vehicles at any time with[V(KK(L- 1) ) 3 V(I) < 
V(KK(l- 1) +16] and the available headway D{I) < SD(I) , 

then the vehicle I ( = KK(l)) reduces its speed to that 
of vehicle KK(L~1) and the vehicles move from A-^ to ^ and 
and to B 2 respectively during the looking for gap 

time T r , if overtaking is possible, the vehicle I starts 
accelerating from A 2 / goes to the other lane and comes 
back to its own lane at A 3 , During the same period (T^) 
the vehicle KK(L-l) moves from B 2 to B 3 . The overtaking 
vehicle occupies the wrong lane for a distance d 2 given by 


d 2 = d + D(I) + S (KK(L-l)) (5.10) 

where d = distance travelled by the vehicle KK(L-l) during 
time ^=[0.28 V(KK(L-1)). ; D(I) = the available 

headway for overtaking vehicle I; and S(KK(L-1)) = minirrum 

spacing for vehicle being overtaken =[o.2V (KK(l-I)) + L 0 ]. 

r 

All the distances are measured longitudinally and the 
lateral movement of vehicles to the wrong lane is ignored , 
The overtaking vehicle(I) is travelling at speed 
V(KK(l~ 1 ) ) from A-^ to A 2 and accelerates from a 2 to A 3 . 


The distance d 2 can be represented as 


d 2 =0.28 V (KK(L-l) ).T 1 + 


1 


a T 


2 

1 


2 


(5.11) 
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where a is the acceleration in kmph/sec. 
From Eqs. 5,lo and 5,11 


and 


1 2 

— a T x = D(I) + S (KK(L-D) 


2(D(I) + S (KK(L-l) )) 


a 


(5.12) 


The total overtaking time needed by the vehicle to 
move from A-^ to A 3 is 

OT(i) = T r + T l (5.13) 

and thr- fineal position (A 3 ) of the overtaking vehicle can 
be given by 


PFIN(I) =P(I) + D(I) +0.28 V (KK(L-l)) 0T(l) 

+ S (KK(L-D) (5.14) 

(ii) when the speed difference between two 
successive vehicles is more than 16 knph ( as in case of 
motor vehicle following a slow moving vehicle) , overtaking 
'vehicle is assumed to travel at its normal free speed. 
During the looking for gap time (T R ) , th<~ vehicle generally 
decelerates and subsequently accelerates to its normal 
speed and hence covers a smaller distance. It is assumed 
that the distance lost in acceleration and deceleration 
is equal to [o.28 T R (V(I) - V (KK(l-I) )] , Hence tho 

relative distance covered by vehicle during overtaking 
cp oration is given by; 
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A = D(I) + S (KK(L-l)) (5.15) 

This distance is travelled at a relative speed of 
W = v(I) - V (KK(L-l)) during time 
T x = ,\/0.28 VV ' (5. 16 ) 

The total overtaking time is given by 0T( I) = T R + T 1# 

The overtaking vehicle occupies the wrong lane 
while moving from a 2 to A 3 for time Tj_. Overtaking is 
possible on] y if there is no conflict between the overtaking 
vehicle and the traffic stream coming from opposite direction. 
All the vehicles of the opposite stream are moved in their 
own lane for T^ seconds to check whether any vehicle lies 
between A 2 to A 3 • After the check/ vehicles of opposite 
stream are moved back to their original positions. If 
there is no liklihood of any conflict, the vehicle 1 moves 
to the wrong lane and comesback after OT(I) seconds to 
position PFIN(I) . 

Sometimes when a fast moving vehicle finds a 
platoon of slow moving vehicles moving ahead of it, and 
the available gap between any two vehicles of the platoon 
is not sufficient for it to reenter, then it desires to overtake 
the platoon in one operation. Fig, 5.4 shows the overtaking 
of a platoon of vehicles. 

Let the vehicle I ( = KK<L )) follow the vehicles 
KK (L-l) ,KK(L-2) , KK(L~3), etc. When vehicle 1 is at 
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and V(I)> V (KK(l~ 1)) / it desires to overtake KK(L-l) . 

The time to overtake KK(L-l) is ; 

OTMD(I) = T r + T x (5.17) 

where T-^ is the time during which vehicle KK(L-l) moves 
from B 2 to B 3 t 1 . g. / PMXD (KK(L-l)) and vehicle KK(l- 2) 
from C 2 to C 3 , i.e., PMID(kK(l- 2 ) ) , The available gap 
for (KK(L-l) ) is; 

D(KK(L-1)) = PMID(KK(L~2) ) - PMID(KK(L-1) ) (5.18) 

If D(KK(L-1) ) is less than (s(KK(l- 1)) + SD(I)) , then the 

gap is not enough for vehicle I to enter between KK(L-l) and 
KK(L— 2) . If V(I)> V (KK(L“2) ) , the vehicle I desires 

to overtake KK(L-2) also. When vehicle KK(L-l) is at B 3 , 
the position of vehicle I at A 3 is given by; 

PMID(I) =[P(I) + 0. 28 V(KK(L-1) ) *0TMD(I) + D(I) 

+ S(KK(L-1) ) ] (5.19) 

1- 

and D( I) = D(KK(L-1) ) - S(KK(L-1)) (5.20). 

The time to overtake KK(L-2) is during which 
vehicle KK(l- 2) and KK(l- 3) move from C 3 to C 4 and D 3 to 
respectively. When KK(l- 2) is at and the available 
gap D(I<K(L-2) ) > pD(I) + S (KK(l- 2))], the vehicle I can 
move between KK(l- 2) and KK(l- 3) to position A 4 . Otherwise 
I will overtake KK(L-3) also. Let T 3 be the time taken by 
vehicle 1 to movt from A 3 to A 4 .Vehicle I enters its own lane at 



position a. 4 given by ; 

PFIN(I) = PMED(I) + D(I) +0.28 V(KK(l-2) ) T 2 

+ S(KK(L-2) ) (5,21) 

after a total overtaking time of 

OT ( I) = OTMD(I) + 0?2 (5.22) 

The above procedure can be adopted for calculating 
overtaxing time and final position after passing a number 
of platooned vehicles. 

5,3.5 Slow Logic for Constrained Vehicles 

Overtaking is not possible When the vehicle desiring 
to overtake is likely to come in conflict with the opposite 
traffic stream. For the restrained vehicle to move in i ts 
own lane, a minimum spacing must be maintained. It can be 
done only if the vehicle reduces its speed to that of 
vehicle ahead, i.e, V(I) s V(KK(L-1)) . The vehicle I continues 
to move at this speed till it is able to overtake and 
maintain its own free speed. such constraints delay the 
vehicles thereby reducing the operating speeds. 

5,4 simulation Logic 

Simulation starts at time T = o second. The vehicles 
entering the section from each direction are identified by 
their sequence numbers and data concerning each of those 
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vehicles including arrival time, category and free speed 
etc. are stored by the memorandum process, initially the 
lanes are assumed to be empty and at each unit time 
interval ( 1 second) , the vehicles are moved according to 
their stream logic defined in Sec, 5.3. To accomplish this, 
the following procedure was adopted, 

XL1 the vehicles moving in a lane, including new 
arrivals and excluding overtaking vehicles at that instant, 
are arranged in the order of their positions. In one of the 
lanes, the arranged vehicles are scanned from the vehicle 
at the head of the stream. The available headway and the 
speed difference between successive vehicles are calculated 
for all other vehicles in the stream. The stream logic is 
determined for each vehicle. For any vehicle desiring to 
overtake, it is checked as to vhether opposite traffic 
comes in conflict with the overtaking operation or not. This 
is done by moving the opposing vehicles on an adhoc basis 
for the duration of overtaking operation at normal speeds. 

If there is no conflict, the vehicle can perform overtaking 
and be at tine predetermined position at the end of its 
overtaking times . It is hence not included in the consideration 
of stream logic for the duration of overtaking operation. 

It is considered again from the moment it reenters its own 
lane, in case of conflict, overtaking is not possible and 
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the vehicle is restrained to move in its own lane at the 
speed of the vehicle ahead. Scanning for a unit interval 
is complete when the last vehicle of the stream has been 
considered. Scanning is similarly done for the other lane. 
The vehicles in both the lanes are moved ahead for one 
second interval at the appropriate running speeds determined 
from stream logic. 

The above procedure is repeated for each time 
interval. Vehicles going out of the section are taken 
out of the system and their characteristics like departure 
time, delay time and operating speed,etc. # are computed and 
stored. Simulation is stopped at an appropriate time and 
group characteristics are computed for each category of 
vehicles . 

5.5 Dara for Simulation 
5,5,1 Weed 

Data concerning the several characteristics and the 
inputs of process are required for computer simulation. In 
case adequate and representative data are available, they 
may be read into the computer. Otherwise, it may be 
necessary to generate one or more sets of representative 
data in the computer for subsequent use in simulation. In 
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such a case, the problem is one of random selection from 
the probability distribution, of the exogeneous stochastic 
variable. 

5,5.2 Pseudorandom Number Generator 

Uniformly distributed random numbers play an 
inportant role in the generation of random variables drawn 
from other probability distributions. The numbers generated 
by computer subroutines are called pseudorandom numbers, as 
they are generated by a completely deterministic relationship 
and yet their statistical properties agree with those of an 
idealised chance device (Naylor ct al.,1971) . 

An ideal pseudorandom number generator should 
yield sequence of numbers that are uniformly distributed; 
statistically independent; reproducible ; and nonrepeating 
for any desired length. Further such a generator should 
bo capable of generating random numbers at high speed and 
require minimum amount of computer memory, Oongruentiol 
methods (Rotenberg, I9 60 ; Naylor et al , , 19 66) are generally 
used for generating pseudorandom numbers in a digital 
conputer. 

Multiplicative congruential method was used for 
generation and it is given by ; 



and 


(5.23) 
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+1 = 7^ n i ( mod m ) i 
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r i+l 


m 


5 3 5 

where /\= 5 and modulus m = 2 -1. starting number n Q 

was taken to be an odd number. By changing the initial 
value n Q , 16 different sequences of uniformly distributed 
random numbers were generated. The generated sequences 
were found to be un correlated and uniformly distributed. 


5,5.3 Random Numbers with Specified Distribution 


Generation of random variables of different dis- 
tributions can be derived from uniformly distributed 
random variables. Let the cfinulative probability function 
F(30 , of a stochastic process be 

x 

F(x) = Prob [x ^ X 1 « / f(t) dt (5.24) 

— oo 

1* 

where F(x) is defined over the range o F(X) ^ 1 and 
f(t) represents the probability density function of random 
variable X at X = t, 

As F( X ) has a uniform distribution over (0-1) range, 
uniformly distributed random numbers r in the (0-1) range 
are generated and set F(X) = r . For any particular value 
of r ( say r Q ) , which is generated, it is possible to find 
the value of X ( in this case 5^ ) corresponding to r Q by the 
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inverse function (Naylor et al, 1966), 

X o = F _1 (r Q ) (5.25) 

- 1 

whore F ( r Q ) is the inverse transformation of r on the 
unit interval into the domain of x. 

5.5.4 Data Generation for Input Variables 

The following three variables are probabilistic 
whose data are to be generated: 

(i) arrival time or intcrarrival time gap between 
consecutive vehi cl ss entering from either direction of the 
road section; 

(ii) category of the arriving vehicles like car, 
truck, bullock cart, bicycle etc.; and 

(iii) free speed of the arriving vehicles. 

Vehicular Arrival Times: The interarrival time gaps between 

the vehicles have been found to follow a shifted exponential 

distribution (Sec. 2,3 ) in which the probability P of a 

gap being loss than t seconds is given by; 

( t -l,o) 

P (h <t ) s [l - e (T/V-l.o) ] (5.26) 

or t = 1.0 - ( T/V - 1.0 ) Log 0 ( 1 - P ) (5.27) 

where ( 1-p ) is a random fraction whose value lies within 
(0-1) range. By Monte Carlo method if a uniformly 
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distributed random number of (0-1) range is generated and 
substituted for ( 1-P) in Eq. 5.27, the interarrival time 
gaps can be computed ( Von Neumann, 1951; Gerlough, 19 64) . 

These time gaps follow the shifted exponential distribution. 
From the sequence of arrival gap lengths t^ , the arrival 
time of vehicles ATj_ = ATj__ ^ + tj_ were computed. Separate 
sequence of arrival times were computed for the two lanes. 

Category of Arriving Vehicles: Vehicles have been classified 
into six different categories. When a new vehicle enters 
the section from either direction, it is necessary to know 
its category so that other characteristics like free speed etc. 
may be generated accordingly. It is possible to specify for 
any category of vehicle a subrange corresponding to composition 
in the overall range of (0-1) for F(x ). Any arbitrary 
allocation of the subranges is possible, for example, the 
Various subranges corresponding to a composition are shown 
in Table 5.2. The subranges for any desired composition 
arc stored. A uniformly distributed random number ( r) is 
generated and its value compared with various subranges. 

When r lies in a particular subrange, the arriving vehicle 
belongs to the category of that subrange. 
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TABLE 5.2 COMPOSITION AND CHARACTERISTICS OF MIXED 

vehicular traffic 


K 

Catecjo ry 
of 

Vchi cle 

Assumed 
Traffic 
Composi- 
tion in 

P or cent 

Sib range 

Mean Free 

Speed 

MVF(k) 

Standard 
Deviation of 
Free Speed 
SDVF(K) 

i 

Cars 

20 

o ,< o.2o 

49 .95 

0.05 

2 

Tru cks 

15 

0 . 2a£ r«C 0.35 

44. 30 

6.80 

3 

Tongas 

5 

0.3 5£r^ 0.40 

12.30 

1.50 

4 

Bullock 

Carts 

10 

0 . 4Q<: r x < o . 5o 

5.20 

1. 10 

5 

Scooters 

10 

0.50<r,< 0.60 

39 .80 

7.60 

6 

Bicycles 

40 

o.eo^r^ 1.0 

U.90 

1.S0 


Free Speed of Vehicles: The free speeds of a particular 

category of vehicle follow a censored normal distribution 

with the range ( ix - 3 sd ) to ( x + 3 sd) . Ihe central 

limit theorem can be used to derive a generator (RN) for 

normally distributed random variables by talcing the sum of 

twelve uniform variables ( Tocher, 1963), i.e., 

12 

^ = ( 2 r j_ - 6.0 ) . Knowing the normally distributed 

i=l 

random variable PN , the free speed of vehicle FS of category 
K was determined as follows: 


FS = MVF(K) + m (SDVF(K)) 


(5.28) 



140 


v/hcre MVF(K) and SDVF(k) are the mean and standard deviation 
of free speeds for category K. 

Uniformly distributed random numbers are needed 
foj. generating input data. In all, sixteen sots of random 
numbers were generated using the same subroutine by changing 
the initial value n Q , viz., 

(i) two sets for arrival gap distribution of either 
di recti on ; 

(ii) two sets for determining the category of arriving 
vehicles, i . e, , one for each lane; and 

(iii) twelve sets for free speed distribution of 
different categories of vehicles ( six for each lane) . 

5.6 Measures of Effectiveness 

Goals and objectives of a study should be e;<pressed 
in explicit terms so that definite measures can be used in 
evaluating the system behaviour. Commonly such measures 
are termed as criteria, figures of merit or measures of 
effectiveness (wohl and Martin, 1967) . These have three 
fundamental characteristics, viz., (i) they should be 
quantitative, (ii) they should measure the effectiveness of 
the entire system; and (iii) they should be efficient 
statistically, i . e, , have a small variance. 
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Quantitative measures of effectiveness that were 
used to validate the simulation model and to estimate the 
flow characteristics are: (i) travel time and hence 
operating speed for each category of vehicle in the mixed 
flow? (ii) proportion of delayed vehicles and delay times 
of different categories? and (iii) density of the section. 

5,7 Initial Conditions 

A process may be steady or transient. A transient 
process has to be studied from its actual initial conditions 
for the desired time period. But for a steady state process, 
it is necessary to ensure that the characteristics of the 
process have reached a steady state. Simulation of the 
steady state process is carried out for a sufficiently long 
period in order to estimate the characteristics of the process. 
It is necessary to select the starting conditions that make 
the transient period short, starting conditions can be set 
by any of the following methods ( Mize and Cox, 19 68 ) : 

(i) Choose starting conditions that approximate 
steady state conditions and insert explicit/ these starting 
conditions in the model; 

(ii) Use a long period simulation so that the data 
from the transient period is insignificant relative to the 
steady state data; 



(lii) Introduce an initial nonrecording period to 
c,,_t the process into steady state conditions. This is done 
by running the process till setady state conditions are 
achieved, then clearing all statistical accumulations ( tut 
leaving the state of system as it is ) and continuing the 
run for a sufficient period. 

Simulation runs are made for two cases, viz., 

(i) for validation of the formulated mathematical model; 
and (iii) for estimation of the characteristics of the 
process under steady state for varying volume levels and 
traffic composition . 

For mixed vehicular traffic at different composition 
the steady state conditions are not known and hence it is 
not possible to adopt the first method of starting with 
known steady state conditions. Second method involves errors 
in estimation of parameters as transient effects may become 
insignificant only for extremely long times. Hence the 
following procedure based on the third method was adopted 
in the study. 

Simulation is started with an empty system. The 
vehicles start arriving from both directions. First few 
vehicles are not significantly affected by the traffic 
stream and move at their free speeds. Constraints start 



143 


building up when vehicles arriving from either direction 
occupy the entire length of the roadway section. Steady 
state conditions are said to exist v/hen rate of arrival 
equals the rate of departure. These conditions are likely 
to occur only when some vehicles have passed through the 
entire section. The initial time period is affected by 
the length of the roadway section. Trial runs were made 
at various volume levels to obtain the steady state conditions. 
It was observed that time period for initialisation is 
also governed by the volume level? being less for high 
volume of traffic. This appears to be logical as with 
higher volume levels/ arrival gaps are less and vehicles 
being closely spaced/ approach the steady state conditions 
soon. Trial runs showed that for a 2.25 km. long section/ 
the maximum time period for initialisation was around 15 minutes 
for a traffic volume of 5o vehicles per hour per lane,whereas 
it was only 5 minutes for volume level of 700 and more. 

Based on the trial runs, initial periods selected for 
simulation experiments are given in Table 5,3. For some 
volume levels and compositions, jamming occurs and hence 
there arc no steady state conditions. In this case also, 
the time porrods for initialisation given in Table 5,3 were 
adopted. 
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table 5.3 initial time periods for obtaining steady 
state conditions 


Traffic Volume 
VPH in Either 

Di r-: cti on 

Initial Time Period in seconds 

1.0 km. Section 

2.25 km. Section 

5o - 100 

600 

900 

150 - 300 

500 

600 

3 5o - 600 

3oo 

400 

650 - 10O0 

200 

300 


5.8 stopping of simulation Runs 

Simulation experiments on a model can be conducted 
either at a particular point in time or over extended periods 
of time, in the former case, simulation is said to be static in 
which sampling takes place across the ensemble, in the 
latter case, the sirrulation is said to be dynamic, i ,e. , time 
series simulation in which sampling takes place over time 
(Naylor et al , 19 71) . A static simulation is achieved by 

replicating a given simulation run by changing the sequences 
of pseudorandom numbers used to generate stochastic variables. 
Dynamic sirrulation results by simply extending the length 
of a run over time. The main problem is how long the sample 
record must be ( Number of replications in static simulation 
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and time period for dynamic sinulation) in order to achieve 
a given level of statistical precision, A number of stopping 
rules are generally adopted and they are available from 
references (Buchhofer, 19 54; Winer, 1962; Gilman, 1968) . 

The sample size may be estimated by (i) finite sampling 
techniques where the sample size is initially determined for 
a given level of accuracy, or (ii) sequential sanpling techniques 
where the sample is tested at frequent intervals to determine 
the error of estimate and to stop sanpling when tho errors are 
within appropriate limits. A sequential procedure was adopted 
in this study. 

The main flow characteristics in this study are the 
delay time and operating speed of vehicles. At the end of 
initialisation, the system is in a steady state and system 
characteristics reach their steady state values. 

They will remain unchanged over further periods of simulation. 

The output characteristics are evaluated only after the 
initial time period. The average value of the characteristics 
were calculated for successive intervals of 5 minute duration. 
Operating speed of cars was used as the dominant characteristic. 

Sinulation was stopped when the average operating speed of 
cars for two consecutive intervals differ by not more than 
0.1 knph, Sinulation time was a function of traffic volume 
and it varied from around 3o minutes for high volume levels 
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of around 800 to more than an hour for volume levels 
of the order of loO„ 

5.9 Formal ati on of Computer Programme 

5.9.1 General 

Flow charts for computer simulation of mixed 
vehicular traffic are shown in Figs. 5.5 and 5.6. The 
computer programme may be written in either a general purpose 
language like FORTRAN or a special sirrulation language like 
GPSS ( Schriber, 19 7o), SUBSCRIPT ( Kiviat, 19 68), DYNAMO 
( Pugh, 1963) or SIMJLATE ( Hold, 1967 ). The special 
simulation languages have been written to facilitate the 
programming of certain types of systems. For example 
SIMJLATE was designed for simulating large size economic 
systems whereas GPSS/3 60 and SIMSRIPT are particularly suited 
for queuing problems. The principal advantage of these 
languages is that when found suitable, they require less 
programming effort and time than for FORTRAN programming. 

5.9.2 Computer programming 

Mixed vehicular traffic flow is quite complex and 
simulation languages available were not considered to be 
suitable for sirrulation. Furthermore, memory requirement 
was also a constraint. Hence FORTRAN language was used for 
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computer programming. Sinulation consists essentially of 
the following: 

(i) The volume level and composition of traffic is 

road; 

(ii) Initially the system is assumed to be empty? 

(iii) Vehicular arrival times and other related 
characteristics are generated and stored for every 5 minute 
interval s; 

(iv) Scanning is done at 1 second intervals. The 
vehicles arc arranged in the order of their positions/ their 
stream logic is determined and they are moved through appro- 
priate distances? 

(v) The process is run for the initial time period 
in order that a steady state may be reached? 

(vi) The process is simulated for 5 minute intervals/ 
at the end of which the output is processed and printed; and 

(vii) Simulation is stopped when the parameter estimates 
are within specified error limits. 

The programme consists of six subroutines and a 
main programme to link them. They include; 

(i) SdBl and SJB2 - to decide the stream logic for 
vehicles in the stream/ 

(ii) RA and NA - to generate pseudorandom numbers 
of uniform and normal distributions respectively? 
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(iii) MEETU - to separate oho output variables like 
average delay, operating speed, free speed, density etc. for 
six different categories of vehicles; and 

(iv) AVSD - to compute the mean and standard deviation 
of various output variables for each category of vehicle. 

The programme involves large number of variables and 
the maximim dimension of any variable representing vehicular 
characteristics can be 100 only. In simulation, a larger 
number of vehicles need to bo processed. A chain was set up 
to avoid memory overflow. As soon as the sequence number 
of arriving vehicles reached 100 , the departing vehicles were 
taken out of the system and sequence number of vehicles alongwith 
their characteristics were changed. New vehicles were generated 
with the sequence number not exceeding 100 at any time. By 
this chain system, the dimension was restricted to 100 without 
affecting the number of vehicles to be processed, 

5.10 Validation of Model 

5.10.1 General 

Formulation of a mathematical model involves the 
definition of a structure and a set of parameters. Validation 
involves the verification that the structure is correct and the 
parameter estimates are reasonable. Validation may also 
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involve parameter estimation. Generally the validation 
of the model is based on available information and an 
agreement betv/oen the throughput or output of the model with 
earlier observations. The better the agreement, better is 
the validity of the model. The model can be used to simulate 
future conditions only after it has been validated and found 
to be a reasonable representation of the actual system under 
study ( Naylor and Finger, 19 67) . 

5.10.2 Parameter Estimation and Validation 

A simulation model and the constituent equations have 
been described earlier in Secs. 5.2 and 5.3. It involves 
assumptions of specific relations for minimum spacing S(I), 
and maximum spacing for overtaking vehicles SD(I) in terms of 
parameters L 0 and T Q ( Eqs. 5.8 and 5.9 ) , and the reaction 
time of looking for gap T R , which is also a parameter. 
Initially the values of l q , t q and T R wore assumed on the 
basis of coirparable I.R. C. ( 19 5o) Specifications for 

overtaking sight distances. Validation was not satisfactory 
and so it was proposed to treat them as parameters and 
estimate their values by simulation. 

From records of field studies, data from six periods 
of one hour each with different compositions and volume 
levels, were selected for parameter estimation and validation. 
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For each hour, arrival time, category and free speed of 
vehicles were generated on the basis of fitted distributions. 
Using estimated values of parameters, the process was sirrulated. 
There are several measures of effectiveness. It was proposed 
to use for validation the frequency distribution of delay 
times for cars and trucks which are most affected by interaction. 
Simulation was done until these parameters could be estimated 
wichin the specified level of error. They were conparcd with 
field data for goodness of fit by 3C - test. When the fit 
is not satisfactory, the parameter values L 0 / T 0 end t r 
were modified till the simulated results were corrqp arable to 
field data. In this study the model was validated after four 
trials. The fitted parameters were as follows: 

(i) Minimum spacing for vehicles; 

S( I) =0.2 V(I) + L q 

where L 0 =6,5 metres for heavy motor vehicles, bullock 

carts and tongas 

= 5.5 metres for passenger cars ; and 
= 4,o metres for scooters and bicycles 

(ii) Headway gap (T Q ) of 2 seconds ,i . e. , sd(I)= 0.56 V(l). 
Overtaking is likely to take place only when the avail- 
able headway equals or is less than 2 seconds. 

(iii) Reaction time (T R ) of looking for gap in the 
opposite stream equals 2 seconds. 
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For the fitted model, validation also indicated a 
general agreement with frequency distribution of operationg 
speeds of different categories of vehicles and the proportion 
of delayed vehicles. The model is considered to be a 
satisfactory representation of the system and so it can be 
used for further computer simulation experiments with mixed 
vehicular traffic. 
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6. SIMJLATIOi.nI ANALYSIS OF MIXED VEHICULAR TRAFFIC 

6.1 Design of Simulation Experiments 
6 . 1,1 Introduction 

The behaviour of a system is affected by a number of 
factors. Simulation essentially deals with the evaluation 
of system response to changes in the factors affecting it. 

The aim of any experimental investigation might be to explore 
and describe the response surface over regions of interest 
in the factor space. This needs observations of the response 
at various factor levels. Some of the criteria for good 
experimental design are as follows ( Cochran and Cox, 19 57 ) : 

(l) The model and its underlying assumptions should 
b e app rop ri at e ; 

(ii) The analysis resulting from design should provide 
unambiguous information on the primary objective of the 
experiment; 

(iii) The design should provide maxclmum information 
v;ith respect to the major objectives and adequate information 
with respect to all objectives of the experiment ; and 

(iv) The design must be feasible within the working 
conditions that excist for the experimenter. 



There are a set of experimental designs that 
provide not only economy in the required number of experi- 
mental trials, but also additional qualities like minimum 
variance estimates, desirable confounding patterns and 
case of computation. The common methods of experimental 
design are full factorial, fractional factorial, rotatable , 
and response surface designs , ( Davies, I960 ; Bonini, 

19 63 ? Hcrzberg and Cox , 19 69 ) . 

6.1.2 Sampling for Parameters 

Interaction between the vehicles affects the basic 
flow charqictcri sti cs of mixed traffic at various volume 
levels. Important parameters include the following: 

(i) volume of the mixed traffic in either direction; 

(ii) its composition, i.e., proportion of the 
different categories of vehicles in the mix ; and 

(iii) the free speed distributions of the different 
categories of vehicles. 

In the present study, the free speed distribution 
parameters of different categories of vehicles were kept 
constant. Only traffic composition and volume levels 
were selected for exploring the response surface. Traffic 
volumes and compositions were also kept same in cither 
direction. 
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The operating speed of the slow moving vehicles 
does not depend generally on the traffic composition and 
volume levels as these vehicles can move etc their free 
speeds irrespective of other traffic factors. On the 

other hand/ the operating speed of motor vehicles is very 
much dependent upon the traffic factors and is below their 
free speeds due to restrained operations. Operating speed 
defines the level of service for any vehicle and interaction 
affects the level of service for cars which have the maximum 
speed among all the six categories. The experiment was 
designed for various combinations of different vehicles at 
varying proportions. For each combination, the response was 
estimated for increasing volume levels till at a particular 
volume, capacity is exceeded, resulting in platooning of 
vehicles. 

6.1.3 Ejqperimental Design for Homogeneous Traffic 

Simulation ejqperiments v;ere designed for homogeneous 
traffic consisting of either passenger cars or trucks only. 
These two vehicles move faster amongst all the vehicles 
and their operating speeds are likely to be affected under 
different sets of traffic conditions. The free speeds of 
vehicles of each category are probabilistic, and there is 
considerable interaction between vehicles of the same type 
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even in homogeneous flow. The volume level for each 
simulation run was increased till jamming occurs. The 
volume levels of passenger cars and trucks used for simulation 
are shown in Table 6.1. 

6.1.4 Eiqpcrimental Design for Mixed Traffic 

Experimental design for mixed traffic requires 
the selection of different combinations of vehicles. Each 
of the combinations is to be ejqperimented at various volume 
levels. Factorial designs require a very large number of 
design points to map the entire response surface. This was 
not possible due to constraints on computer time. Based on 
different traffic compositions/ a few selected design points 
were used to estimate the interactions between different 
combinations of vehicles. 

The interaction between passenger cars and any 
other category of vehicle was evaluated for all the five 
possible combinations, viz., each case included passenger 
cars and one of the remaining five categories with different 
proportions of vehicles in the mix. The proportions selected 
for design were 30 and 50 percent of trucks ; 25 and 5o percent 
of scooters; 25 , 5o and 75 percent of bicycles; 10 , 25 and 
5o percent of bullock carts; and lO and 25 percent of tongas. 
The volume level for each set was increased till it exceeded 
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TABLE 6.1 FACTOR COMBINATIONS FOR SIMULATION OF MIXED 
VEHICULAR TRAFFIC 


Combinations 


Homogeneous Cars 

Trucks 


Two Vehi cl e Cars 

Combination s 

70.0 
5o .0 

Cars 

90.0 

75.0 

Cars 

90.0 

75.0 

50.0 

Cars 

75.0 

50.0 

Cars 

75.0 

50.0 

25.0 


Traffic Composition 
in percent 


Tru cks 

30.0 

50. 0 

Tongas 

10.0 

25.0 

Bullock Carts 

10.0 

25.0 

50.0 

Scooters 

25.0 

50.0 

Bicycles 

25.0 

50.0 

75.0 


Volume Level (V) in 
Ei th er Di rection- VP H 


50/ 100, 150, 200, 300, 
400, 500, 600, 650, 700, 
800 


50, 100, 150, 200, 300, 
400, 500, 600 


50, 100, 200, 300, 400, 
500, 600 

50, 100, 200, 300, 400, 
500 


10O, 200, 300, 400, 500 
lOO, 200, 300, 400, 450 


100, 200, 300,3 50 
100, 200, 300 
lOO, 200 


100, 200, 250, 300, 400, 
600, 800, 850 
lOO, 200, 300, 400, 600, 
800 , 900, 1000 


lOO, 200, 300, 400,600, 

aoo 

10 o, 200, 300, 400, 600, 
800, 900 

lOO, 200, 300, 400, 600, 

800 / 900, 1000 


Contd. . . . 



table 6.1 ccmtd 


Combinations 

Traffic Composition 

Volume Level (V) in 


in percent 

Either Di recti on-VpH 

Three 

Cars 

Trucks 

Tongas 


Vehicl e 





Combinations 

63.0 

27*0 

10.0 

100, 200,300 


45.0 

45.0 

10.0 

100, 200, 300 


52. 5 

22.5 

25.0 

10 0, 200,300 


37.5 

37. 5 

25.0 

100, 200, 300 


Cars 

Trucks 

Bullock Carts 



63.0 

27.0 

10.0 

100, 200, 300 


45,0 

45,0 

10.0 

100, 200 


52.5 

22.5 

25.0 

100, 200 


37.5 

37.5 

25.0 

100, 200 


35.0 

15.0 

50.0 

loo 


Cars 

Tru cks 

Scooters 



52.5 

22. 5 

25,0 

100 , 200, 300, 400, 600 


37.5 

37,5 

25.0 

100, 200, 300, 400, 6oO 


35.0 

15,0 

50.0 

100 , 200, 300 , 400, 600, 80o 


25.0 

25.0 

50.0 

100, 200, 300, 400, 6oO,QOO 


Cars 

Trucks 

Bicycles 



52.5 

22.5 ' 

25,0 

10 0, 200, 300, 400, 600, 8 00 


37.5 

37,5 

25,0 

100, 200,300, 400, 600 


35.0 

15.0 

5o.O 

100, 200, 300, 400, 600,800 


25.0 

25.0 

50.0 

100, 200,300, 400, 600 


17.5 

7.5 

75.0 

100, 200, 300, 400, 600, 800 


12,5 

12.5 

75.0 

100, 200, 300, 400, 600, 800 

Six Vehicle 

Cars= 

17. 5, Trucks=7*5, 

100, 200,300, 400 

Combinations 

Tongas=6.o , 

Bullock 



Carts 

=4.Q, Scooters=l5«0/ 



Bicycles=5o . 

0 



Cars= 

15,0/ Trucks=l5,0/ 

100, 200,300, 400 


Tongas=6,0, 

Bui lock 



Carts 

-4.0, scooters^lO.O, 



Bicycl es=5o . 

0 




162 


capacity. The various sanple points aro also given in 
Tabl e 6.1. 

The interaction between three categories o£ vehicles 
was evaluated for four combi nations. These include cars, 
trucks and one of the remaining four categories of vehicles, 
viz., bullock carts, tongas, scooters and bicycles. Mixed 
traffic flow affects the characteristics of fast vehicles 
like cars and trucks and these two vehicles also interact 
with each other to a considerable extent. The car - truck 
ratio was maintained as 7o/3o and 5o/5G as in the case of 
two vehicle interaction studies. The combined effect duo to 
rhe three vehicles was evaluated for different proportions 
of the third vehicle in the mix for the two car - truck 
ratios in each of the four combinations. The proportions 
of the third vehicle in the total volume were 25 and 5o 
percent of scooters ; 25 , 50 and 75 percent of bicycles ; 
and 10 and 25 percent for both tongas and bullock carts. 

They were simulated at various volume levels given in 
Table 6.1, 

The interaction between four and five categories 
of vehicles in the mix were not investigated. Yet the 
interaction between all the six categories of vehicles in the 
traffic mix was estimated at the two observed peak traffic 
compositions at different volume levels ( Table 6, l ) . 
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In all about 180 experimental runs for different 
compositions and volume levels of various combinations of 
vehicles were made requiring more than 100 hours of processing 
time on IBM 7044/1401 system at I.I.T. , Kanpur, 

6,2 Simulation Characteristics 

Simulation for each set of data was started on an 
empty system till steady state conditions were achieved at 
the end of the initial time period. Simulation was continued 
further and output variables were determined at the end of 
every 5 minute interval till the parameter estimates were 
within specified error limits. The specific characteristics 
obtained from simulation runs include: 

(i) proportions of each category of vehicles delayed 
and their delay times; 

(ii) Operating speeds of different vehicles ( VEF ) 
and also of different categories, ( AWE }; 

(iii) Density of the section or the number of vehicles 
of each category per unit length ; and 

(iv) Rate of departure of ouch category of vehicles 


from the section 
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6.3 Homogeneous Traffic Flow 

In homogeneous traffic , the interaction is between 
vehicles of the same category. Simulation of homogeneous 
traffic consisting of only cars or trucks was performed 
at various volume levels. 

6.3.1 Homogeneous Car Traffic 

Simulation results of homogeneous car traffic are 
shown in Figs. 6.1 to 6.5, At low volume levels, say upto 
150 vehicles per hour ( VPH ) for either direction, the 
available headways are more and very few manoeuvring 
operations are restrained. This results in very samll 
proportion of the vehicles being delayed ( Fig. 6.1 ) ; 
the average delay for all the vehicles at low volumes is very 
small ( Fig. 6 . 2 )there are very few vehicles in the section, i . e. , 
density is low ( Fig, 6.3 ) ; and operating speed of cars 

is very near their free speeds ( Fig. 6,4 ) . As volume 
increases, the available spacings get reduced and in turn 
a large number of overtaking vehicles are restrained. The 
increased volume thus increases the delay and density 
thereby considerably reducing the operating speed. 

Simulation results ( Figs. 6.1 to 6 . 5 ) indicate the 
following four distinct characteristics: 
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(i) Upto a volume level of l5o VPH in either 
direction, the rates of increase of delayed vehicles and 
average delay time are very small. Density also builds 
up slowly and operating speed - volume relationship is 
nonlinear with operating speed very near the free speed. 

(ii) Beyond volume level of 150 VPH , there is a 
sudden change in the proportion of delayed vehicles, density 
increases sharply, delay times increase, and the operating 
speed also decreases at a rapid rate. Between volume 
levels of 15o to 6o0 , there is a regularity in the relation- 
ships between the flow characteristics. Operating speed 
and the density change linearly whereas the proportion of 
the delayed vehicles and the delay times increase non- 
1 in early with volume. 

(iii) Between volume levels of 600 and 650 , the 
proportion of the delayed vehicles, delay times and the 
density increase rapidly. Hie operating speed also 
decreases at a faster rate, 

(iv) At a volume level of 700 VPH , traffic flow 
is always in the transient state and steady state is never 
reached. The characteristics of flow change with time ? 

j 

are indicated by marked open circles in Figs. 6.1 , 6.3 
and 6,4. These indicate that density of vehicles is 
building up in the road because of jamming with a consequent 
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increase in clelny time and decrease in operating speed and 
hence the volume. 

Based on sirrul ation results# the following mathema- 
tical relationships have been established by regression 
an al y si s » 

'The proportion of the delayed vehicles is given 
by 

, , 2 

P d = 0.135 ( V-100 ) +0,0001 ( V-100 ) 

for 150 ^ V 6o0 ( 6.1 ) 

where = percentage of delayed vehicles and V = traffic 
volume in either direction ( VPH ) . 

The average delay per vehicle increases with volume 
level and is represented as 

, , 2 

d = 0.0045 V+ O.00002( V- 100 ) 

for 150 ^ V 600 ( 6.2 ) 

where d = average delay in seconds. 

The volume density relationship shown m -Fig, 6,3 
indicates an abrupt increase in density between volume 
levels of 150 to 2Co # though the rate of increase in 
density up to volume level of 150 is same as that between 
200 to 600 VPH. The relationships can be represented as: 

D. = 0.054 v for V *< 150 ( 6.3 ) 

t > 
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and D t = 13 + 0.054 ( V-200 ) 

= 2.2 + 0.054 V for 2o °\^ v \< 600 ( 6,4 ) 

where D-)- is the deasity in terms of vehicles per km. length 
o f road s ecti on . 

Beyond volume level of 600 , the density increases 
rapidly and at 700 VPH , the density continues to increase 
with time while the outgoing volume rate decreases indicating 
jamming. Jamming seems to occur when density exceeds around 
40 vehicles per lane per km. of roadway section. 

The speed volume relationship ( Fig, 6.4 ) indicates 
that up to a volume level of 150 VPH , the operating speed 
decreases nonlincarly and is very near the free speed. It 
can be represented as 

0S HT = 50 - 0,0014 V - 0.000019 v 2 

for V ^ 150 < 6.5 ) 

where OS^-p = operating speed ( kmph ) of homogeneous car 
traffic having mean free speed of 5o kmph. 

Between traffic volume of l5o to 65o , the relation-* 
ship can bo given by 

0^ = 49.36 - 0.0101 ( V - 150 ) 

= 50.875 - 0.0101 V 

for 150^ Vv< 650 


( 6.6 ) 
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For a volume of 700 VPH, the operating speed 
decreases to 43.28 kirtph after 15 minutes of simulation and 
to 39 knph after 30 minutes. The outgoing traffic rate 
is around 650 VPH at the end of 3o minutes. With an input 
volume of 800 VPH the jamming starts in just 5 minutes of 
simulation and operating speed decreases to 39 krrph in 
10 minutes. It is generally seen that jamming starts when 
operating speed of cars is around 43.30 kmph. 

The speed density relationship can be derived from 
Eqs» 6.3 to 6.6 as follows: 


and 


D-^ i^t 

08^ = 50 - 0.0014 V ( ) - 0,0000 19 ( 

0.054 0.054 

( 6.7 ) 


2 


for D t 8 

D t _ 2.2 

OSum = 50.875 - 0.0101 ( ) C 6.8 ) 

0.054 

14 >< D <*< 40 


Fig 6.5 shows that the speed density relationship 
is nonlinear at lower density ( less than 8 ) and is linear 
for densities between 14 to 4o vehicles per lane per km. 
This agrees with Eqs. 6.7 and 6.8, 


Capacity and Level of Service: The concept of level of 
service can be defined best in terms of measures familiar 
to drivers. H.R.B, ( 1965 ) specifies the level of service 
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by; (i) operating speed.; and (ii) ratio of service volume 
to capacity. 

In order that a highway provides an acceptable 
level of service, it is necessary that service volume is 
less than capacity. Service volume is thus the maximum 
volume that can be carried at any selected level of service. 

The speed volume and volume density relationships 
are not same at various volume levels of traffic flow. 

These may be classified into the following levels of 
service : 

(i) Free flow; Upto a volume level of 15o VPH, 
the operating speed is very near the free speed, density 
is low , and only less than 5 percent of vehicles arc delayed. 
Thus upto this volume level , there is little or no 
restriction to the manoeuvrability due to other vehicles 
in the stream. 

(ii) Stable flows This zone lies between 150 to 
650 VPH. Here manoeuvrability begins to be restricted 
by other vehicles due to interaction and the operating 
speed decreases linearly with volume. Yet steady flow is 
achieved and maintained. 

(iii) Unstable flows This zone lies between volume 
levels of above 650 to less than 700 , There is a significant 
change in operating speeds and about 5o percent of the 



vehicles are delayed. Volume speed relationship is 
nonlinear and this indicates that flow is approaching 
capacity. 

(iv) At volume level of around 700 or more, the 
outgoing volume ( or service volume ) is less than incoming 
volume and this continues to decrease with time thereby 
increasing density and reducing operating speed. Speed and 
service volume may both drop to zero under jarrming conditions 
The capacity of the 2 lane 2 way highway can thus he 
said to be between 65o to 700 VPH. However at a volume 
level of 700 , the outgoing volume rate reduces to 65o only 
after half an hour. Considering this to be an acceptable 
level of service, the capacity can be considered to be 
700 vehicles per hour in either direction or a total of 
1400. This may be compared with the capacity of 2000 
specified by HRB ( 1965 ) , The difference may bo attributed 

to lower free speed of passenger ears in India and perhaps 
a differing stream logic. Various levels of service and 
t'neir characteristics are presented m Fig, 6.4 and Table 6.2 

6.3.2 Homogeneous Truck Traffic 

The general characteristics of homogeneous truck 
traffic are similar to those of homogeneous car traffic. 

Speed volume relationship for homogeneous truck obtained 
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table 6.2 characteristics of homogeneous car traffic at 
various levels of service 


Level 

of 

Service 

Arriving 

Volume 

Maximum 

Service 

volume 

(MSV) 

MSV 

Capacity 

Operating Speed 
knph 

Free 

FI ow 

0-150 

150 

4 0.214 

50-0.0014 V-0 .000019 V 2 

Stabl e 

Flew 

150— 650 

650 

^ 0.857 

50 .8 7 5-0.0 101 V 

Un sta- 
bl e 

Flow 

650-700 

^<700 

4 1.00 

Variable 

Forced 

Flow 

> 700 

Widely ^ 1,0 

Variable • 

Capacity 

Approaches 

Zero 

May drop to nearly 
zero 
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from simulation is shown in Fig. 6.6. When traffic volume 
is low, i.e., upto 100 VPH , the available headways arc 
largo, there is hardly any constraint to manoeuvrabili ty 
of vehicles and the operating speed of truck is very near 
the free speed. The speed decreases nonlinearly with 
volume and can be represented by; 

2 

0S T rucks = 44,3 “ 0.0009 V - 0.000043 V 

for V N < 100 ( 6.9 ) 

Between volume levels of loo to 450 , the operating speed 
deceases due to smaller spacings and more interaction at 
higher volume levels. The operating speed of trucks in 
the zone of stable flow is given by; 

0S T rucks = 43 * 78 “ 0,0143 V ( 6. 10 ) 

for 100 < V < 450 

The flow characteristics are unstable at volume 
levels between 450 to 500 VPH. At volume level of around 
5oO and beyond, the operating speed continues to decrease 
with time alongwith service volume. They may both drop to 
zero when traffic gets completely jammed. The capacity 
of the 2 lane 2 way highway may be considered to be around 
475 trucks per hour for either direction or a total of 9 50. 

The various levels of service and their characteristics 

k> 

are summarised in Table 6.3 
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TA3LE 6.3 CHARACTERISTICS OF HOMOGENEOUS TRUCK TRAFFIC AT 

various levels of service 


Level 

of 

Arriving 

Volume 

Maxirrum 
Servi ce 

MSV 

Operating Speed 
ktrph 

Service 

Volume 

(MSV) 

Capacity 

Free 

Flow 

0 - 100 

100 

\< 0.211 

44. 3-0, 0C09V-O. 000 043 V' 

Stable 

Flow 

100 - 450 

450 

^0.947 

43,78-0.0143 V 

U n st a- 
bie 

Pi ow 

450 - ( < 50 0) 

475 

N < l.o 

Variable 

Forced 

Flow 

> 500 

Wi del y 
Variabl e 
Capaci ty 
Approaches 
zero 

<< 1.0 

May drop to nearly 
zero 
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6.4 interaction Between Two Categories of Vehicles 
6.4.1 Introduction 

VJhen traffic is composed of two categories of 
vehicles, i,e,, passenger cars and those of another category, 
the flow characteristics like operating speed, density, 
delay etc. are affected by the following interactions i 

(i) between cars themselves; 

(ii) between cars and vehicles of the other category; 

and 

(iii) between vehicles of the other category. 

Thu manoeuvrability of cars is restrained by vehicles 
that move slower than cars. As the speed difference increases, 
the number of desired passing manoeuvres increase and they 
are further restrained by vehicles of the opposite stream. 

If a slov/ moving vehicle comes in the zone of conflict of 
an overtaking vehicle, it takes more time to cross the 
zone thereby restraining the overtaking vehicle and causing 
more delay. The interaction between the two categories 
of vehicles is likely to be affected by the composition 
of these vehicles in the mix. The interaction between slow 
moving vehicles is small as there is less of overtaking 
among them. The interaction between cars themselves was 
discussed in Subsec. 6.3,1, The maximum delay is likely 
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to be for cars clue to their high speed. The operating 
speeds of cars at varying volume levels and proportions 
of the two vehicles in the mix were estimated from simulation. 

6,4.2 .Speed Volume Relationships for Different Combinations 

Cars and Trucks: The free speed distributions of cars and 
trucks overlap over a sufficiently large extent as the 
difference between their mean free speeds is only 5.7 krnph. 

A number of trucks may thus have more speed than some of 
the cars. Trucks being longer need a larger minimum spacing. 

So car and truck traffic cannot be categorised as homogeneous. 

Fig. 6,7 shows the operating speed of cars at 
various volume levels and for different car truck compositions. 
The operating speed decreases nonlineariy with volume and 
may bo represented as folloy/s: 

0S MT2 = 57 * 7 V" 0 * 038 for 50 % <V^ 300 { 6. 11 ) 

for ear - truck ratio of 70/30 
0S MT2 ” 59,8 v’" 0 * 0 ^ 4 for 50^ V ^ 300 ( 6.12 ) 

for car - truck ratio of 5o/5o 
0Spn?2 = operating speed of cars in the mixed traffic 
V = total traffic volume . 

There is a considerable reduction in operating 
speeds even at low volume levels as interaction builds up 
quickly. Furthermore, for the some volume of total traffic. 


and 

where 

and 



1 ! 



•67 SPEED VOLUME RELATIONSHIP FOR CAR TRUCK 
COMBINATIONS 
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operating speed decreases with increased proportion of 
trucks in the mix as cars arc restrained by a larger number 
of trucks in the stream. The operating speed of cars is 
thus a function of volume level and traffic composition. 

When there are 30 percent of trucks m the mix, the operating 
speed continues to decrease with time ar a volume level 
of 600. Even after 10 minutes of simulation, service 
volume falls to 56o and decreases further alonguith operating 
speed. So for this composition, the road can carry per lane a 
maximum traffic volume of 5oO. with 5o percent of trucks, 
jamming starts occuring at a volume level of 4o0 after 
25 minutes. Thus capacity is only 400 VPH for this 
composition . 

Gars and Tongas: There is a wide difference in the free speed 
°f cars and tongas and there is very little possibility 
of any delay to tongas. The main interactions are between 
cars themselves and between cars and tongas. Fig. 6,8 
shows that for lo percent of tongas, there is a linear 
relationship between the operating speed of cars and total 
traffic volume ( V ) up to 400 VPH, Beyond volume level 
of 4oO , operating speed decreases nonlinearly and jamming 
starts occuring at a volume level of 5oO after 3o minutes. 

When there are 25 percent tongas in the mix, the operatirtg 



If!/, 



S-6'8 SPEED VOLUME RELATIONSHIP FOR CAR 
TONGA COMBINATIONS 



speed decreases nonlinearly with volume and jajnming starts 
oc curing at a volume of 450 VPH. The road can thus carry 
500 and 45o VPH respectively with tonga proportions of 
10 and 25 percent. 

The speed volume relationships for the stable 
zone can be represented as: 

0S MT2 = 50 - °.°129 V for V x < 400 ( 6. 13 ) 

for 10 percent tongas 

0S MT2 = 61 ' 1 V~°* 069 for 400 ( 6. 14 ) 

for 25 percent tongas 

Cars and Bullock Carts j Fig. 6,9 shows the speed volume 
relationships for varying coirposi tions of bullock carts 
in the mix. Hillock cart, being the slowest vehicle causes 
maximum delay to the fast vehicles like cars. In case a 
car is unable to overtake, it has to reduce its speed to 
that of bullock cart, thereby considerably reducing its 
operating speed. The operating speed decreases exponentially 
with traffic volume. Even with only lo percent bullock 
carts in the mix, the roadway can carry only 300 VPH and 
at a volume level of 3 5o , jamming builds up immediately. 

As the proportion of bullock carts increases, the operating 
speed and capacity are considerably reduced. 



■+ ARRIVING VOLUME 

■► + SIMULATION TIME 

NOTE: PERCENTAGES ARE IN TERMS 
OF TOTAL VOLUME 
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Cars and Scooters: The mean free speed of scooters and 
motorcycles is less than that of cars by 10.3 kmph and there is 
a very small overlap of their free speed distributions. 

Thus very fati scooters may overtake cars. scoot er, is 
a two wheeler ,it occupies very little pavement width and 
its simulation logic ( sec. 5.3 ) is different from that 
of four wheelers. A car may overtake a scooter without 
moving into the wrong lane. However, if a scooter overtakes 
a car, in travels in the wrong lane during the overtaking 
operation. Fig. 6,lo shows that v;ith 25 percent scooters 
in the mix, the operating speed is very close to the free 
speed upto a volume of 2oO VPH, The operating speed of 
cars is given by; 

0s MT 2 = 50 - 0.0008\r - 0,000012 V 2 ( 6. 15 ) 

f o r V 200 

Beyond a volume level of 200 , the operating 
speed reduces linearly upto a level of 8 00, at which 
jamming starts building up after more than half an hour. 

The relationship is given by; 

0S MT2 = 5l * 14 - °-°°9 v C 6.]6 ) 

for 200^< 800 

Vflien there are 50 percent scooters in the mix, 
free flow occurs upto a volume level of 250 and operating 
speed changes linearly from 250 to 900 , Jamming occurs 
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immediately at volume 1 -svel of 1000, The operating 
speed of cars is given by 

0S MT2 “ 50 " 0.0008 V - 0.000007 V 2 ( 6.17 ) 

for V ^ 250 

and °^MT2 = 51 ‘ 50 “ 0.00856 V ( 6, 18 ) 

for 250^ 900 

The capacity is 85o and 900 VPH respectively with 
25 and 50 percent of scooters in the mix. 

Fig 6.io show s that operating speed of cars at 
a parti ail ar volume level is more with higher percentage 
of scooters in the mix. This is in contrast to that for 
the combination of cars with trucks/ tongas or bullock 
carts ( Fig. 6.7 to 6.9 ) / in which operating speed of 
cars decreases with increased proportion of the second vehicle. 
With increased proportion of scooters in the total volume, 
there are a smaller number of cars and so there is less 
interaction between cars and scooters and between cars them- 
selves. In a car truck combination/ a higher proportion 
of trucks may reduce the interaction amongst cars; but 
there is considerable increase in the car - truck interaction 
and between trucks themselves. Consequently the operating 
speed of cars is considerably reduced. 
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Car and Bicycles: The bicycles are two wheelers and have 
simulation logic similar to scooters. The mean free speed 
of bicycles is rruch less than that of cars and there is 
no overlap of their free speed distributions. The main 
interaction is thus between the cars themselves. Pig, 6.11 
shows that when there are 25 percent bicycles in the mix, 
the cars move under free flow conditions up to a volume of 
200 VPH. The free flow is up to 25 q and 300 VPH respectively 
with 50 and 75 percent bicycles in the mixed flow. The 
operating speeds of cars in the free flow are given by the 
following relationships; 

°Smt2 =* 50 " 0.0016 V - 0,000009 V 2 ( 6. 19 ) 

for V^< 200 

for 25 percent bicycles 

OS W2 = 50 - 0.0013 V - 0,000006 V 2 < 6.20 ) 

for V ,< 250 

for 5o percent bicycles; and 

0S MT2 * 50 " 0.0010 V - 0.000005 V 2 ( 6.21 ) 

for V ^ 300 
for 75 percent bicycles. 

At higher volume levels the flow is stable and there 
is a linear reduction in operating speed with the volume 
level upto 800 VPH, The operating speed of cars in the 
stable zone is given by 




0S MT2 = 51. li - 0.00905 V for 200 ^ V ^ 8 00 

( 6.22 ) 

for 25 percent bicycles 

0S MT2 = 51,39 - °.°0836 V for 250^ V s < 8 50 

( 6.23 ) 


for 5o percent bicycles ; and 

OS MT2 = 51,71 - °*°0 82 V for 3o0 ^ V ^ 900 

( 6.24 ) 

for 75 percent bicycles. 

Fig. 6.H shows that operating speed of cars at 
a particular volume 3 evel is more with higher proportion 
of bicycles in the mix. Maximum service volumes that 
approach capacity are 800 , 8 5o and 10OO VPH respectively 
at 25 , 5o and 75 percent of bicycles in the mix, 

The above results demonstrate clearly that intera- 
ction is a function of volume level and traffic composition 
Mathematical equations derived may be valid only for specifi 
conditions and corrposi ti on s. Interaction between cars them- 
selves have been estimated in Sabsec. 6.3.1. When there are 
two vehicles, viz., cars and vehicles of another category, 
there is a need to estimate the interaction between: 

(i) cars and vehicles of the other category like 
trucks, tongas, bullock carts, scooters and bicycles; and 
(ii) among the second category of vehicl 6s them- 


selves 
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These interactions reduce the operating speed from 
that for homogeneous car traffic. The change in operating 
speed is a function of the volume of cars and the proportion 
of the second category expressed as percentage of car volume . 
Hence Fig. 6.12 shows the operating speed volume relation- 
ship for cars for different compositions of the second 
category of vehicles. These may be compared with Figs, 6.7 
to 6.11 where the relationships are in terms of total volume 
of traffic. The effect of interactions is more clearly 
evident in Fig, 6,12, Equations were fitted to define the 
speed volume relationships of cars for different compositions 
of two vehicle combinations and are given in Table 6,4. it 
may be noted that in Fig. 6. 12, the plotted points are 
derived from simulation and plotted curves correspond to 
fitted relationships ( Table 6,4 ) . 

6.4.3 Maximum Service Volumes for Different Combinations 

For a given composition of two vehicle combination# 
there is a maximum possible service volume which defines 
the roadway capacity. Five curves each for a two vehicle 
combination are shown in Fig. 6.13 and these define the 
capacity of a two lane two way highway. Car scooter and 
car bicycle curves are very close to each other at the 
Various compositions. Maximum service volume is a little 
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6.4 SPEED VOLUME RELATIONSHIPS OF CARS FOR TWO 
VEHIOjE COMBINATIONS 


P roportion 
of second 
Vehicle as 
Percentage 
of Car Volume 

Operating Speed of 
krrph 

Cars ( 

■ ° S MT2 * 


2 



- 

50-0. 0014 V-j-0. 000019 V x 

for 

Vj^iso 


50.875-0.0101 V l 

for 

150%£ 


V = Car Volume 



' 42.9 

57.6 Vf 0 ' 0405 

for 

35 ^ M 230 

100.0 

57.6 Vl -0.0525 

for 

25 N < V^150 

11.1 

50-0,01433 V a 

for 

V 1 ^ 360 

33.3 

56.7 v^ 0 ' 0427 

for 

V 2 ^300 

11.1 

-0.047 

58. 2 V x 

for 

V l v< 270 

33,3 

58.2 v^O.065 

for 

V 1 < 75 

33.3 

51. 14-0.012 V-j^ 

for 

150 ^ V ^ 600 

100.0 

51. 50-0,0171 V x 

for 

125 ^V-^400 

33.3 

51.13-0.01207 V ± 

for 

150 ^ V^SOO 

100,0 

51. 39-0.0167 V 1 

for 

125 < v l < 400 

300.0 

51. 7l-0. 0328 V 1 

for 

75^ V-^2 25 
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COMPOSITIONS of two vehicle combinations 


more for car scooter combination than car bicycle 
combination at the same composition. This indicates a 
larger interaction between cars and bicycles than between 
cars and scooters. For vehicles of wide base and for the 
same proportion, maximum service volume increases as the 
relative speed between the car and the second vehicle decreases, 
viz., from car bullock carts combination through car - 
tonga combination to car truck combination. 

6.4.4 a Multiplicative Model for Speed Volume Relationship 
of Mixed Traffic Flew 

The results of simulation analyses indicate that 
for a given volume of cars, the operating speed of cars 
decreases with increased proportion of the second category 
of vehicles. The interaction factor, IF , represents the 
ratio of operating speed of cars in mixed traffic ( ^ 

to the operating speed of cars in homogenecus traffic 
( 0s HT ) , i.e., 

IF = ° SmT2 ( 6.25 ) 

OS HT 

The operating speed of cars in mixed flow can be 


represented as; 
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where = volume of cars in the mixed flow; P 2 = proportion 
of the second vehicle expressed as a percentage of ? 
and Cat =the second category of vehicles in mixed flow. 

The interaction factor can also be r&p resented by 

IP « f 2 ( v x f P 2 , CAT ) ( 6.27 ) 

From Eq. 6.25 and 6.27 

0S M?2 = 0S HT ( V 1 ) • f 2 ( V l' P 2 ' CAT } < 6 ' 28 ) 

Hence the operating speed of cars in mixed traffic 
is a product of the operating speed of cars alone and the 
interaction factor, IF , which depends upon the volume of 
cars ( 7 a nd the type ( CAT ) and proportion ( P 2 ) of 

vehicles of second category. This is referred to as 
a mil tipi icative model for speed volume relationship in 
the mixed traffic flow. It may be noted that the multi- 
plicative model is designed to agree with relationships for 
homogeneous and mixed traffic flows. 

From the results of simulation analysis given 
in Table 6,4 and Fig. 6.12, the foil wing relationships were 
developed for different combinations of vehicles. 

Let p, represent the percentage of vehicle of 
category ( . ) in the two vehicle combination. Let 
subscript T , KK / BC , SC and BY represent respectively 
trucks, tongas, bullock carts, scooters and bicycles. 
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Then for; 

(i) Car - truck combination 


V, 


IF = 1.0 - 0.0012 P T + (0.0023+0,00112 P T ) ( 10) Uo3 + 0,4666P t ) 

( 6.29 ) 

(ii) Car - tonga combination: 


V- 


XF = 1,0 -(0.020 + 0.036 P KK ) [ 1 + 10 


(293 + 13.8 6 P 1GC ) 


] 


( 6.30 ) 


(iii) Car - bullock cart combination: 


V, 


IF = 1.0 - 0.009 P BC + 0.0062 P Bc j0.O) ( 230 + 26.1 P BC ) ] 

{ 6.31 ) 


(iv) Car - scooter combination: 


IF = 1.007 - (0.0007+ 0.000043 Pg^ (lo/ 7 °° ~ 2 ' 8 PsC> t ] 

( 6.32 ) 


(v) Car - bicyc] e combination: 

V 1 

P BY (486-0. 97 5P by ) 

XF = 1.0 - [ - 0,0015 + 0,000126 ( ) ( lO) j 

1 + 0.01 p by 

( 6.33 ) 

Though simulation analysis performed in this study 
is fairly extensive it had to be limited to only a few 
proportions for each combination because of limitations of 
computer time. Yet it indicates clearly the nonlinearity 
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o Jr the interactions and furthermore there is a general 
consistency in the results. Hence Eqs. 6.29 to 6,33 may 
be considered satisfactory even though further refinement 
or parameters and equations for IF is possible. 

Using the Eqs. 6.29 to 6,33 , the interaction factor 
IF for different combinations and for various proportions 
were derived and they are shown in Figs. 6. 14 to 6.18. They 
may be used as design charts for traffic having characteristics 
similar to those considered in this study. In the absence 
of additional in formation, it may be presumed that , IF 
is a function of V-^ , P 2 and Gat only, and is independent 
of the free speed distribution of vehicles. In such a case, 
the speed volume relationship of cars alone for any other 
free speed distribution of cars can be established and 
used with the design charts ( Figs. 6.14 to 6.18 ) to deter- 
mine the speed volume relationship for mixed traffic flow. 

6,5 Interaction Between Three Categories of Vehicles 

6,5.1 Introduction 

Simulation was carried out for combinations of three 
categories of vehicles and was limited to car, truck and 
one of the remaining four categories, viz., tonga, bullock 
cart, scooter and bicycle. In such combinations, the flow 



1 100 200 300 

CAR VOLUME IN VFH 

'14 INTERACTION FACTOR FOR TRUCKS ON 
CARS 



>•16 INTERACTION FACTOR FOR BULLOCK 
, : CARTS ON CARS 
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characteristics are affected by the following*. 

(1«) Interaction between cats themselves and trucks 
themselves? 

(ii) Interaction between cars and trucks? 

(iii) Interaction of vehicles of third category 
among themselves? and 

(iv) Interaction between vehicles of third category 
with cars and trucks. 

6.5.2 Speed Volume Relationships for Different Combinations 

The flow characteristics were estimated by simulation 
for each combination at different compositions and volume 
levels C Table 6, l ) . Figs 6, 19 to 6. 22 show the operating 
speed of cars as a function of total traffic volume for 
different compositions of various three vehicle combinations. 
The results indicate that in case bullock carts or tongas 
are present alongwith cars and trucks , the operating speed 
of cars decreases with volume due to smaller headways and 
large number of constraints, furthermore/ for the same 
volume level the operating speed decreases with increased 
proportion of bullock carts or tongas, indicating a signi- 
ficant interaction between slow and fast moving vehicles. 

The operating speed also decreases with increased proportions 
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of trucks in the mix. 

When the third category of vehicle is scooters or 
bicycles, the operating speed is less with a smaller 
proportion of scooters or bicycles. This is because two 
wheelers like scooters or bicycles occupy very less width 
and interaction with fast vehicles is not very significant. 
The results agree with those of two vehicle combinations. 

The speed volume relationships are generally exponential 
and they are not reported in this study. 

6,5.3 Feasibility and Levels of Service for Three Vehicle 
Combi nations 

The operating speed of cars varies with volume of 
traffic, the combination of vehicles and their coirposition . 
In this study, it is proposed to define the level of service 
in terms of meag operating speed of cars. Hence for homoge- 
neous traffic, the level of service worsens as the volume 
of cars increases, and in mixed traffic, for the same volume 
of cars , the level of service worsens with the proportion 
of other categories of vehicles. Different combinations 
and compositions of vehicles may represent the same lav cl 
of service and so they may be connected by isolines or 
isosurfaces of Same level of service. 
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The level of service may be considered as tho 
response of the system to different combinations of vehicles 
and for three vehicle combination, it can be plotted in 
an isometric graph. The plots for three vehicle combinations 
with bicycles and bullock carts as the third vehicle , art- 
shown respectively in Figs. 6.23 and 6,24. It may be 
noted that the feasible combinations of traffic are bounded 
by the throe planes and an upper concave surface defining 
the capacity of a two lane two way highway. The surface 
defining capacities is bounded by its intersections with tho 
bounding planes, which in turn define tho capacity curves 
for two vehicle combinations. 

The level surfaces for car truck bicycle combination 
for operating speeds of 48 , 46 and 44 kmph are shown in 
Fig. 6.23. It may be noted that level surfaces are also 
concave from above and may be bounded by the bounding 
planes or the capacity surfaces. For car truck bullock 
cart combination, only tho level surface for 46 knph and 
level curves of 48 and 44 knph are shown in Fig. 6.24 to 
avoid confusion. They are similar to those of bicycles 
in Fig. 6.23 , but the level surfaces , curves and capacity 
surfaces for bullock carts are much steeper than in case 
of bicycles. 
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6.5.4 Conditional Multiplicative Model 

Spaed volume relationships indicate that the 
general characteristics of three vehicle combinations are similar 
to those of two vehicle combinations, for which a rrultipli- 
cative model has been developed. The operating speed of 
cars for different combinations depends upon the traffic 
composition and volume level of each category. Based on 
Eq. 6.26 , it can be expressed as ; 

OS HT3 = f ( V x , P 2 - p 3 - cat ) ( 6.34 ) 

vjhere = car volume in the total mix; and P 2 and P 3 are 
proportions ( expressed as percentage of car volume ) of 
trucks and vehicles of third category ( CAT ) respectively. 

Assuming that the interactions of trucks and the 
third category of vehicleson cars arc independent of each 
other, the multiplicative model was generalised as follov/ss 

0^MT3 = 0 S HT< V 1 ) • ( V x , P 2 >« lF 3 * V 1,P3' CAT ) 

( 6.35 ) 

where XF 2 represents the interaction factor for trucks 
on cars and IF 3 , the interaction factor for third category 
of vehicles on cars, They are estimated from Eqs, 6.29 to 

6,33, 

The operating speed of cars for different three 
vehicle combinations an'd at different compositions and volume 
obtained from simulation and were compared 


1 ovel s. 


wire 
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with the values calculated by Eq. 6.35. The results were 
not consistent indicating that the interaction between 
trucks and third category of vehicles is also significant. 
Hence the multiplicative model ( Eq. 6.35 ) was modified 
as a conditional multiplicative model given by 

0 S MT3 = OS HT ( V x ). IF 2 ( V L , P 2 ) . IP 3 [ < V 1 + V 2 ), P 3 /CAT] 

C 6.36 ) 

I 

where P 3 is tho proportion of the third category of vehicle 
m terms of car truck volume ( + v 2 ) t and 13 

the conditional interaction factor for the third category 
of vehicles on the fast moving vehiclC3,cars and trucks. 

Assuming that IF 3 has the same parameters and 
forms as in Eqs. 6,29 to 6.33 , but with V-^ replaced by 
car truck volume ( ) » the operating speed of car& 

in mixed traffic were calculated using Eq, 6,36 for different 
three vehicle combinations at different proportions and 
volumes. A sample calculation is shown in Tablu 6.5. They 
were conrpared with results from simulation analysis. 

(Table 6.6 and Figs. 6.25 to 6,28 ) . The errors in the 
operating speed of cars has an average absolute value of 
around 0.2 percent with a maximum of around 0.5 percent* 

There is thus a very good agreement between simulation 
results and the fitted model confirming the validity of 
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TABLE 6.5 SAMPLE CALCULATIONS FOR OPERATING SPEED OP CARS 
IN MIXED VEHICULAR TRAFFIC 

(i) Traffic composition : scooters a 25 % ; 

Car truck ratio = 70/3 O; 

Hence Cars = 52.5 % ; Trucks = 22,5 % 
and Scooters = 25, o A 


(ii) Traffic Volume (V) = 200 

Car Volume (Vi) =105 ; 
Truck Volume(Vo) = 45 ; 

Scooter Volumefv 3 ) =50 


(iii) Operating speed of cars for homogeneous traffic ( OS^t, ) 
at volume of 105 VPH ( Eq. 6.5 ) = 49.64 knph 

(iv) percentage P 2 of trucks in terms of car volume 
= 30/70 = 42.9 % 


(v) interaction factor, IP 2 for trucks on cars ( Eq. 6.29 ) 
= 0.9593 


(vi) Percentage P 3 of scooters in terms of car truck volume 
= 25/75 =33.3 % 


(vii) Interaction factor, IF 3 y for scooters on car truck 
combination { Eq. 6,32 ) = 0.9969 


(viii) Operating speed of cars ( osj*^ ) in mixed flow 
= OSj-jt * IF2 * IP3 
= 47.43 

(ix) From simulation results, °S^t3 = 47.52 


( x) Error = 47,43 - 47,52 = -0.09 
= 0, 189 % 



217 


TABLE 6.6 COMPARISON OF SIMULATION RESULTS WITH THOSE 
COMPUTED FROM MULTIPLICATIVE MODEL FOR THREE 
VEKECLE COMBINATIONS (EQ. 6.36 ) 


Combination Composition Volume Operating speed of 
in percent Cars (OS Jir ^) kmph 

From From 

Model Simulation 


Error 


Cars, Trucks 

Cars 

=63.0 

100 

47. 27 

47,33 

-0.06 

and Tongas 

Tru cks 

=27,0 

200 

45.85 

46.01 

-0. 16 


Tongas 

=10 . o_ 

300 

44.45 

44,56 

-0.11 


Cars 

=52. 5~ 

100 

47. 23 

47,05 

0.18 


Trucks 

Tongas 

=22.5 

=25.0 

200 

45.13 

45.27 

-0. 14 


Cars 

Tru cks 
Tongas 

=45. o" 
=45.0 
=10.0 

100 

46. 20 

46.25 

-0.05 


Cars 

Trucks 

Tongas 

=37.51 
=37.5 
=25. Oj 

100 

l 

46.3 2 

46,18 

0.14 

Cars, Trucks 

Cars 

=63.0 

100 

45.75 

45.94 

-0.19 

and Bullock 
Carts 

Trucks 

Bullock 

=27.0 

cartsj 

200 

=10 

44.04 

44.15 

-0. 11 


Cars 
Trucks 
Bud lock 
Carts 

=45.0"* 
=45.0 
=10 . 0_ 

100 

44.71 

44.50 

0, 21 

Cars, Trucks 

Cars 

=52.5 

100 

48.52 

48.46 

0.06 

and Scooters 

Trucks 

= 22.5 

200 

47.43 

47.52 

-0.09 


Scooters =25.0 

300 

400 

600 

46.7 2 
46.03 
44.75 

46.98 

46.10 

44.80 

-0. 26 
-0.07 
-0,05 


Cars 

=35.0 

100 

48.90 

49.06 

- 0.16 


Trucks 

=15.0 

200 

47.81 

48.02 

-0,19 


Scooters =50 . 0 

300 

400 

600 

47.07 

46.46 

45.14 

47.21 

46.51 

45,09 

-0.14 

-0.05 

-0,05 

Contd. , . 
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table 6,6 coutd... 


Combination 

Corrpoci tion Volume 

in percent 

Operating Speed of 
Cars (OSy^) kmph 

From From 

Mo del Simj 1 ation 

Error 

Cars 

, Trucks 

Cars 

=37 , 5 

100 

47.58 

47. 33 

0. 25 

and 

Scooters 

Trucks 

=37.5 

200 

45.76 

45.69 

0.07 



Scooters 

=25.0 

300 

44,61 

44.57 

0.04 





40 0 

43.82 

43,70 

0.12 



Cars 

=25. 0 1 

100 

48. 24 

48. 25 

-0.01 



Tru cks 

=25.0 

200 

46 , 61 

46.44 

0. 17 



Scooters 

=50. 0_ 

400 

44.48 

44.65 

-0.17 

Cars 

, Tracks 

Cars 

=52 . 5 

ICO 

48. 51 

40.65 

-0.14 

and 


Trucks 

=22,5 

200 

47.41 

47 . 40 

o.oi 

Bicycles 

Bi cycles 

=25.0 

300 

46.68 

46,69 

-0.01 




400 

45.96 

45.93 

0.03 





600 

44.51 

44.41 

0 . io 



Cars 

=35 . 0 

100 

48.87 

48.82 

0.05 



Trucks 

=15.0 

200 

47.77 

47.75 

0.02 



Bicycles 

=50.0 

300 

47.02 

46.96 

0.06 




40 0 

46.40 

46. 22 

0. 18 





600 

44.99 

45.11 

-0. 12 



Cars 

=17 , 5 -1 

100 

49 . 41 

49 , 48 

-0.07 



Tru cks 

= 7.5 

200 

40.48 

48 . 66 

-0. 18 



Bicycles 

=75. 0_ 

300 

47.68 

47.85 

-0. 17 




4-00 

46.93 

47.08 

-O. 15 





600 

45.32 

45. 28 

0.04 



Cars 

=37.5" 

100 

47 . 57 

47,48 

0.09 



Trucks 

=37.5 

200 

45.74 

45.73 

o.oi 



Bicycles 

=25. 0_ 

' 300 

44.58 

44 , 49 

0.09 




400 

43.75 

43.77 

-0.02 



Cars 

=25.0 

100 

48. 21 

'18 . 22 

-0.01 



Trucks 

=25.0 

200 

46.58 

46.49 

0,09 



Bicycles 

=50. 0_ 

300 

45.36 

45, 23 

—0.13 




400 

44.41 

44,48 

-0.07 



Cars 

=12.5 

200 

47.83 

47.71 

0.12 



Tru cks 

=12.5 

300 

46,75 

46.8 2 

0.07 



Bicycles 

=75.0 

400 

J 600 

45.76 

43.73 

45.73 

43.89 

0 . 03 
- 0.16 
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the conditional nrultiplicative model ( Eq. 6,36 ) , 

6,6 Interaction Between six Different Categories of Vehicles 

On the basis of earlier results, a generalised 
conditional rru 1 tipi i cat ive model was formulated to define 
the operating speed of cars in mixed traffic having several 
categories of vehicles ( NCAT ), namely ; 

NCjff 

0 %iTtJCAT = 0S HT ( V x ) , IF 2 ( V x , P 2 ) Tf IF ± [ ( V +V, ), 

i =3 -L ^ 

pJ_, CA^ ] ( 6.37 ) 

where is the percentage of vehicles of i th category- 
in terms of car truck volume . 

It may he noted that homogeneous traffic, two or 
three vehicle combinations considered earlier are special 
cases of this model. 

Simulation analysis was carried out for six vehicle com- 
bi nations only for two specific compositions and for differenc 
volume levels. The results were conpared with the operating 
speeds predicted, by the conditional rnilriplicative model 
( Table 6.7 ). The absolute errors had an average value of 
0.3 percent with a maxinum of around 0.6 percent. Even 
though they are slightly larger than for three vehicle 
combinations ( Table 6.6 ) , these results are considered 
to be quite satisfactory. 
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Table 6,7 comparison of simulation results with those 

COMPUTED FROM MULTIPLICATIVE MODEL ( E Q* 6. 37) FOR 
COMBINATION OF six vehicles. 


Corrposition 
in percent 


Interaction Factors 
L evel s 

at Various 

Volu me 




Volume Levels 

100 

200 

300 

Cars 

=17.5 





Trucks 

= 7.5 

if 2 

0.9925 

0.9870 

0.9684 

Tongas 

= 6,0 

if 3 

0.9906 

0.9821 

0.9743 

Bullock Carts = 4,o 

if 4 

0.9469 

0.9393 

0.9323 

Scooters 

=15.0 

I f 5 

0.9970 

0.9966 

0.9962 

Bicycles 

=50,0 

if 6 

0.9916 

0.9097 

0.9875 

OSht 



49.97 

49.93 

49 .87 

°SmT3 

from Model 

46.34 

44.84 

43.80 

°Smt3 

from Sinulation 

46.17 

45.10 

44,07 

% Error 



0.37 

“0.58 

-0,61 

Cars 

=15.0 





Tru cks 

=15.0 

if 2 

0.9868 

0.9646 

0.9469 

Tongas 

= 6,0 

if 3 

0.9895 

0.9802 

0.9719 

Bullock Carts = 4.0 

if 4 

0.9531 

0.9449 

0,9376 

Scooters 

= 10.0 

IF 5 

0.9983 

0.9980 

0.9977 

Bicycles 

=50 

if 6 

0,9921 

0.9921 

0,9879 

OSht 



49.97 

49.94 

49.89 

0 S MT3 

from Model 


46.05 

44.09 

43.63 

OS MT3 

from sinulation 

46. 21 

44. 16 

43,61 

% Error 



-0.35 

-0.15 

0,05 
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Ibo results of this section indicate that separa- 
bility of the interacting factors ( XF i ) in terms of a 
conditional multiplicative model is valid. So the genera- 
lised conditional multiplicative model can be used for 
estimating the operating speed of cars in mixed traffic 
flow for different combinations, compositions and volume 
levels, and perhaps even for different free speed distribu- 
tions of vehicles. 

6.7 Characterisation of Mixed Traffic Flow 
6.7.1 Introduction 

The characterisation of mixed traffic flow is quite 
complex due to wide variations of characteristics of vehicles 
in countries liXe India, The characterisation of mixed 
vehicular traffic is normally done in terms of passenger 
car equivalents ( PCEs ) of different, vehicles. The PCES 
are defined as, "Nurriber of passenger cars displaced in the 
traffic flow by another vehicle under the prevailing roadway 
and traffic conditions'*. In this study, the mixed flow is 
characterised in terms of PCEs which express the relative 
effects of each type of vehicle on the operating speed of 
passenger cars, Uie speed volume relationships derived in 
Secs, 6.4 to 6,6 were used for characterisation of traffic 


flow 
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The operating speed of cars in the mixed flow is 
a nonlinear function of traffic volume and its composition. 

The equivalence between mixed traffic and homogeneous traffic 
is in terms of the level of service. Hence pCUS equivalent 
to any given mixed traffic is represented by the volume 
of homogeneous car traffic which has the same operating 
speed. Knowing the composition of mixed traffic, the 
PCES for different categories of vehicles can be derived 
from the PCUS for mixed traffic, 

6.7.2 Characterisation of Two Vdiicle Combination 

PCU S are defined explicitly in terms of by 

Eqs.6,5 and 6.6, viz,, 

OSyjp = 50 - 0.0ol4( PCUS) - 0.000019 ( pCUS ) 2 

for PCU S ,< 150 ( 6.38 ) 

and 

OSj^ = 50,875 - 0.0101 ( PCUS) 

for 150 x < PCU S^ 650 ( 6.39 ) 

Given 0SMT2 ' the corresponding PCUS(2) may be established 
from these equations, say, 

PCUS( 2 ) = \ ( OSj^ )]= f” 1 [OS HT (V 1 ).IP(V 1 ,P 2 ,CAT)] 

C 6.40 ) 
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where PCU S ( 2 ) is the equivalent car volume for a two 
vehicle combination. The PCU S of various car truck combi- 
nations at different volume levels and compositions are 
shown in Fig, 6,29. 

The PCEs of the second category of vehicles in a 
mix can be estimated as follows: 


p C U S ( 2 ) - 

P C E S ( 2 ) 

( P 2 V /1O0 ) 


( 6.41 ) 


where p 2 is the percentage of second category of vehicle 
in terms of car volume V 1# The PCEs for different categories 
of vehicles in 2 vehicle combinations are shown m Figs. 6.30 


to 6.34 at different volume levels and proportions. It is 
assumed that j aiming occurs vhen PCU S equal TOO VPH, 

The results indicate that for tru cks, hullock carts 


and tongas, the PCES reduce with increased car volume. 
Furthermore, they decrease as their proportion in the mix 
increases. This is because the relative interaction effect 
of these vehicles are smaller at higher volumes ( Fig, 6. 12 ) 
of cars and themselves. On the other hand for scooters and 
bicycles, the PCES increase with car volume and vith 
higher proportion of scooters and bicycles in the mix. There 
seems to be a tendency for the PCEs to approach 1 at higher 
volume of cars and higher proportion of second category o 1 
vehicles, viz., for wide base vehicles like trucks, bullock 
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carts and tongas, which have a PCES larger than 1 , the PCEs 
decreases and xn case of two wheelers like scooters and 
bicycles, which have a PCES of less than 1 , the PCES 
increases with car volume and proportion of the secon i 
vehicle. The effect of the proportion of vehicles is 
rruch smaller in case of PCES of trucks than in case of 
other vehicles. 

The range of PCEs for different categories of 
vehicles are given in Table 6 , 8 . The constant values 
suggested by CRRI and t40T are valid utmost over a small 
range of combinations and these ranges are also indicated 
in Table 6 . 8 . It is clear from Pigs, 6.30 to 6.34 and 
Table 6.8 that PCES are variable functions of composition 
and volume level of mixed traffic. 

6.7,3 Characterisation of Three Vehicle Combination 

Vtfhen there are three categories of vehicles in 
the mix, the operating speed of cars for car truck 
combinations ( ) an ^ three vehicle combinations ^ os t-pr 3 ^ 

can be derived respectively by Eq, 6.28 and 6,36, The 
PCus for mixed flow for car truck combinations PCUS ( 2 ) 
and for three vehicle combinations PCUS ( 3 ) can be 
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T/^BLE 


Category 

of 

Vehicle 


Trucks 

Tongas 

Bullock 

Carts 

Scooters 

B1 cycl gs 


.8 RANGE OF PCES FOR DIFFERENT CATEGORIES OF 
vehicles at varying traffic COMPOSITIONS 


CRRI 


Passenger Car Equivalents (PCES) 


MOT 


From sinulation /nalysis 


proportion Car Volume PCES 

as percentage Range in Range 

of Car Volume VpH 


2.7 

2.6 


10 . 7 

0.2 

0 . 4 


3,0 


6.0 

6.0 


1.0 


1.0 


20 - 100 

50 - ISO 

9.0 - 3.5 

150 

25 - 75 

8.0 - 5.0 

10 - 20 

50 — 3 50 

20.0 - 4.7 

20 - 30 

50 - 3Q0 

13.0 - 4.1 

30 - 50 

50 - 250 

3.7 - 3.5 

10 - 20 

75 - 300 

33.0 - 12.0 


25 _ 

50 

50 

- 

500 

0.2 

- 

0.4 

50 - 

10 o 

50 


400 

0.2 

- 

0.5 

100 - 

150 

50 

- 

300 

0 . ’ 

— 

0.37 


25 

50 


500 

0.2 



0,45 

25 - 

50 

50 


450 

O.n 

- 

0.60 

50 - 

100 

53 

— 

400 

0,2 

- 

0.60 

10 o - 

150 

50 

— 

350 

0.2 

- 

0,60 

150 - 

200 

50 

— 

300 

0.2 

- 

0.60 

200 - 

300 

50 

- 

200 

0.2 

- 

0.70 
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determined from the respective operating speeds of cars 
by the equations 

PCU S ( 2 ) = f' 1 ( , ( 6 .„ 2 ) 

and PCUS (3) = r 1 ( 0S MT3 ) ( 6.43 ) 


Tiic PCES for trucks can be estimated from Eq,6,4l. 

The PCES of the third vehicle can be determined as follows 


PGES (3 ) 


PGUS (3) - PCUS ( 2) 

___ 


{ 6,44 ) 


where volume of third vehicle 


V 3 =P 3 (V 1+ ^> (6.«, 

The PCEs of the third category of vehicles for 
different percentages of car truck volume are indicated 
in Figs. 6,35 to 6,38 and these may be used as design charts 
when appropriate. The variations of PCES^ are similar to those 
of the same category for two vehicle combinations. 


6,7,4 Characterisation of Multivehicle Combinations 

When the number of categories of vehicles in the 
mix exceed three, they can be considered in a specified 
order and PCES for each category of vehiclescan bu calculated 
as follows; 

PCUS U) - PODS (i . l) 


PCES (i) = 


( 6,46 ) 




w 


TRUCK VOLUME 20% CAR VOLUM 

TRUCK VOLUME 60% CAR VOLUM 

TRU CK VOLUME 1O07oCAP volume 

NOTE : PERCENTAGES ARE IN TERMS 

qp car +t ruc^ volume 


\ 


\ 

\ j, \ 

\ 

\ 




' -jqc 150 

VOLUME of OARS r TRUCKS 


, G pcE5 OF BULLOCK CARTS IN CAR TRUCK AND 
BULLOCK CART COMBINATIONS 
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where volume of i th category of vehicle V,- =? p ] (V. + ^ 2 \ 

1 11 100 ‘ 

It should be noted that PCEs estimated depend 
upon the sequential order in which vehicles ate considered. 
It is suggested that car truck combination be considered 
first/ the vehicles of wide base considered next and finally 
two wheelers, furthermore/ in each subgroup, the vehicles 
with higher average speeds be considered ahead of vehicles 
with lower speed, viz., in the order trucks, tongas, bullock 
carts, scooters and bicycles. 

The results of this study indicate that the 
nonlinear interaction between different categories of 
vehicles may be represented by a generalised conditional 
mil tiplicative model for operating speed of cars which 
defines the level of service. This can be used to define 
the e<gi±valent volume of cars (PCUS) for the same level 
of service and in turn to determine the variable PCEs of 
different categories of vehicles in mixed traffic flow. 



7. SUMMARY, CONCLUSIONS AND SUGGESTIONS FOR FUTURE STUDY 


7.1 Summary 

Characterisation and estimation of traffic are 
two important aspects in the planning and design of a 
highway transportation system. In developing countries 
like India, the traffic is of mixed nature cons^ing of 
both slow and fast moving vehicles between which there 
are wide variations in speed. Generally constant values 
of PCES of different categories of vehicles are specified 
by appropriate authorities. 

This study considers the nonlinear interaction 
between the different vehicles in mixed traffic flew situa- 
tions and attempts to characterise the variations of PCES, 

i 

in terms of traffic composition and volume. Traffic flow is 
stochastic in nature, the free speeds are probabilistic and 
the process is affected by logical decisions concerning 
acceleration, retardation, overtaking etc. Analytic 
techniques of traffic flew theory, are suitable for the 
analysis of si np 1 e homogeneous vehicular traffic and it 
does not seem feasible to use exact analytic procedures 
for the analysis of complex mixed traffic flow. Qonputer 
simulation is adopted in this study to analyse the mixed 
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traEfic flow and infer the PCES under different volumes 
and compositions of traffic. 

Traffic flow on a highway is a stochastic process 
with seasonal, weekly and daily variations. Using historical 
monthly and daily traffic data of seven years on five 
main highways approaching Kanpur, stochastic models have 
been developed for the monthly and daily traffic data for 
three categories of goods carriers. Trends, annual, six 
monthly and weekly cycles may be present in the data series 
and the residual series can be represented by a first 
order AR model. Alternatively, monthly series can be 
represented toy a multiplicative seasonal ARIHA model of 
order ( 1, 0, 0 ) x ( 0, l, l ) ^ Stochastic models can 

also be used for making forecasts about traffic volumes. 

Detailed data were needed for simulation and so 
field studies were conducted for ten days round the clock 
on a 2.25 km stretch, of Grand Trunk Road. Prom the obser- 
vations, peak hourly volumes; variations of volumes within 
the day of six vehicular types; and interarrival time gap 
distributions were estimated. The free speeds, delay time 
and operating speeds of different categories of vehicles 
were also estimated at varying volume levels and traffic 
compositions. 
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A conputer simulation model of mixed traffic flow 
on a two lane highway was formulated. Vehicles move from 
either direction of the roadway section and scanning is done 
at one second intervals. Any vehicle in the roadway moves 
at its free speed where it had enough headway no travel 
unimpeded. In case it cannot move unimpeded and it has a 
higher speed than the vehicle ahead of it, then it may try 
to overtake. Overtaking is possible only if there is no 
conflict with the opposing traffic stream during overtaking 
operation. Otherwise the vehicle is forced to reduce its 
speed to that of the vehicle ahead until sometime later 
overtaking is possible. The above logic is applied to each 
of the vehicles in the section and those entering the section 
for each time interval. As a vehicle leaves the section, 
time of leaving is noted and its characteristics like 
travel time, delay time and overall running speed etc., 
are determined. 

Initially historical data of arrival times, traffic 
composition and system parameters were used to identify the 
system and validate the system model. The roadway is 
divided into one metre long sections and vehicles are 
moved from either direction at one second intervals. 

Vehicles move at their free speed when there is enough 
headway to travel unirrpeded for two seconds. The change 
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of speed of vehicles is assumed to he instantaneous and a 
looking for gap time of two seconds is provided for 
overtaking operations. The overtaking vehicle accelerates 
only if the speed difference with the overtaken vehicle 
is less than 16 kmph. Otherwise overtaking takes place at 
normal speed. For vehicles being overtaken/ minimum spacing 
related to speed and length of vehicle, is specified. The 
results from the simulated model were consistent with the 
observed characteristics and the model is thus validated. 
Simulation is hence found to be a versatile tool for under- 
standing of the flow process and identifying the system 
model , 

Traffic on a highway is highly nonstationary with 
hourly trends, and cycles, and persistence. However, for 
design considerations a stationary peak rate corresponding 
to a specified level of risk may be assumed. So this 
study deals with simulation of stationary traffic flow. 

It may be noted that when appropriate, more complicated 
non station ary processes can be simulated to derive conpara- 
ble results. Generated data were used in simulation of 
stationary mixed traffic flow at different volumes and 
compositions. When generated data are used in conputer 
simulation, it is necessary to initialise the system. In 
this study initially the stretch of the road was assumed to 
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be empty. The process was ran until a steady state was 
reached. This initial time period varied generally with 
the volume level, being more for lower traffic volumes. 
Simulation was carried out further for a sufficiently long 
time for estimation of the characteristics of the process. 
Simulation was performed for homogeneous, tvo, 
three and six vehicle combinations for specific conpositions 
and volume levels. The characteristics observed include 
(i) proportion of delayed vehicles; (ii) average delay time; 
(iii) operating speed of different categories of vehicles; 
and (iv) density of section. 

The study indicates that, a nonlinear interaction 
between different categories of vehicles in mixed traffic 
flow can bo analysed by simulation; the capacity and 
operating characteristics of highways can be established; 
mathematical models can be formulated for representing the 
interactions; and PCUS for mixed traffic and PCES for each 
category can be derived as a function of traffic volume and 
composition. 
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7.2 Conclusions 

On the basis of this study, the following 
conclusions can be drawn; 

(i) Traffic flow on a highway is a stochastic 
process with seasonal/ weekly and daily variations. Trends 
and annual, six monthly and weekly cycles may be present in 
the data series and the residual series can be represented 
by a first order AR model. Alternatively, the monthly 
series can be represented by a multiplicative seasonal ARIMA 
model of order < 1, 0, 0 ) x ( o, 1, 1 ) 2 . stochastic 
models can be used for making forecasts about traffic 
volumes. 

(.ii) A simulation model which includes unitrpcded, 
overtaking and restrained stream logics is formulated 
Using field data. Parameters of the model are estimated 
and the simulation model is validated, simulation leads 
to a better understanding of the corrponents and interactions 
of the cortplex process , It facilitates mathemati ceil 
modelling of the system and the process. 

(iii) Traffic flow can be classified into four 
distinct levels of service, viz., free flow, stable flow, 
unstable flow and forced flow, simulation helps m identifi- 
cation of these levels for homogeneous and mixed traffic 
and also for estimation of highway capacity. 
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(iv) The level of service can be specified 

in terms o E operating speed of passenger cars in mixed 
traffic flow. Surfaces of equal levels of service can also 
be derived for mixed traffic flow, 

(v) From results of simulation,, it is possible to 
derive regression equations relating volume, operating 
speed and/or density. Relationships for homogeneous and 
two vehicle combinations are given in this study, 

(vi) Simulation analysis of homogeneous car traffic 
and of two vehicle combinations indicate that interaction 
is a function of volume level and traffic composition. 

The operating speed of cars in mixed flew ( 0s^2 ) is a 
product of operating speed of cars alone { ) and the 

interaction factor, IF, v;hich depends upon the volume of 
cars ( Vjj ) , and the type ( CAT ) and preportion of 

vehicles of second category. Mathematical relationships 
for IF' s have also been derived, 

(vii) For three or more vehicle combinations 
a generalised conditional multiplicative model with 
separable interacting factors/ IP 1 s, of the following type 
has been formulated. 


ncat 


OShthw = «wr < v ! >• ir 2 « V p 2 > • TT Mi<V v 2 ' p i' CM 'i> 

i=3 


where - operating speed of cars in mixed flov; 
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having N CAT categories of vehicles / P2 = track volume v 2 
expressed as a percentage of car volume V-j. ; = proportion 

of ith category of vehicle ( CAT^ ) expressed as percentage 
of car truck volume ( Vj + V 2 ) / and IF^ = interaction 
factor as derived for two vehicle combinations but in terms 

r 

of ( v 1 + v 2 ) and P i , 

The above model has been validated from simulation 
results for three and six vehicle combinations for a number 
of compositions and volume levels. 

(viii) The generalised model can be used for estimating 
the operating speed of cars in mixed traffic flow and also 
for determining equivalent passenger car volume ( PCU s) , 

Tliis can in turn be used to derive PCES for different 
categories of vehicles. Design charts for PCEs for two and 
three vehicle combinations have also been drawn. 

(ix) For tracks/ tongas and bullock carts/ the 
PCES decrease with increased car volume and with increased 
preportion of these vehicles in the mix, on the other hand 
for scooters and bicycles. PCES increase with car volume 
and with higher preportion of scooters and bicycles in 
the mix. Thera seems to be a tendency for the PCEs to 
approach unity at higher volume of cars and higher proportion 
of second category of vehicles. 
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(x) The study hence clearly demonstrates the 
nonlinear interaction in mixed traffic flow; that these 
interactions can he defined in terms of the generalised 
model? and that PCU s and PCEs can be derived from such a 
modal , 

(xi) It seems possible to combine stochastic 
modelling and simulation results for the design of highways. 
For cxanple, from stochastic models of different categories 
of mixed traffic, the A DT of each category for a given 

level of risk may be determined. From the composition of 
vehicles at critical penod(s) and the derived ADT< s, the 
volume of each category of vehicle in the critical period(c) 
may be estimated. From the PCES for different categories 
of vehicles for given volume and composition, the PC US of 
mixed traffic at critical penod(s) can be determined, Ihe 
highways can then be appropriately designed. 

7,3 Suggestions for Future study 

Based on the results of this study, the follov/ing 
suggestions are made for future work in this area: 

( 1 ) This study deals with the simulation of mixed 
traffic flow on a two lane highway and is limited to specific 
conpositions, free speed distributions/ straight section 
with no lateral obstructions and equal directional 
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di stribu tioh of traffic. While the results of the study 
indicate the nonlinear interactions, it seems necessary 
to extend the study to different corrposi tions, different 
free speed distributions of vehicles, varying geometries 
like grades, curves, lateral obstructions etc,, varying 
directional distribution of traffic, breakdown of vehicles 
etc., and perhaps for intersections and intercity traffic. 
The results of the study including the mathematical models 
may need modifications because of these factors, 

(ii) It seems necessary to develop a problem 
oriented simulation Language for the simulation of mixed 
traffic flow in muitilane highways, Scuh a programme 
will be helpful in simulating complex mixed vehicular 
traffic on highways with different number of lanes and 
other geometries. 

(iii) Stochastic modelling of some categories of 
traffic has h^on indicated in this study. It seems desirable 
to develop multivariate stochastic models which relate 
temporal Variations of mixed traffic and preserve corre- 
lations among tho different categories of vehicles. Using 
such models, data can be generated and critical periods 
may be identified. Using detailed generated data for 
these critical periods, the mixed traffic flow can be 
sirrulated for a particular design of highway system* The 
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risk associated with the design may be evaluated and this 
can bo used in an iterative design of the system. 

(iv) The capacity of highways is generally given 
in terms of homogeneous traffic. Using simulation, it is 
possible to derive the capacity and level surfaces for 
different combinations of mixed traffic. For combinations 
of more than throe categories of vehicles, a graphical 
representation is not possible and it may be necessary to 
either prepare sets of graphs or combine the vehicles to 
atmost three groups of vehicle combinations* The variability 
of capacity and level surfaces as a function of the free 
speed distributions, geometries etc., also need further 
in vestig ation . 
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